Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi1,311,8300AG20.2% 332.1 / 79.8 158intergenic (+22/+18)tonB/yciAmembrane spanning protein in TonB‑ExbB‑ExbD transport complex/acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base A (54/72);  new base G (15/17);  total (69/89)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-03
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 158 total aligned reads displayed.

TGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAA  >  NC_000913_3_pae_pgi/1311682‑1311978
                                                                                                                                                    |                                                                                                                                                    
tGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGAc                                                                                                                                                     <  1:1535274/150‑1 (MQ=255)
   gAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCgg                                                                                                                                                  >  4:1018626/1‑150 (MQ=255)
         gTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCtt                                                                                                                                            <  2:1776915/150‑1 (MQ=255)
         gTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtt                                                                                                                                            <  2:319563/150‑1 (MQ=255)
          tttGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCttt                                                                                                                                           <  2:1272064/150‑1 (MQ=255)
          tttGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCttt                                                                                                                                           <  4:730399/150‑1 (MQ=255)
          tttGAGCGTGAGGTGAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtttt                                                                                                                                          <  2:994367/150‑1 (MQ=255)
           ttGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCCCCGAGATTCAGTAAGCAGAAGGTCAAAAGCCTCCGACCGGAGGCtttt                                                                                                                                          <  2:533606/150‑1 (MQ=255)
           ttGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCtttt                                                                                                                                          <  1:1265115/150‑1 (MQ=255)
           ttGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtttt                                                                                                                                          <  3:361384/150‑1 (MQ=255)
            tGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                                                                         <  1:1745334/150‑1 (MQ=255)
            ggAGCGTGAGGTGAAACATGCGATGCGCAGATGGCGTTATGAGGCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTg                                                                                                                                         <  1:782765/149‑1 (MQ=255)
             gAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGa                                                                                                                                        <  1:363488/150‑1 (MQ=255)
             gAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGa                                                                                                                                        <  1:937512/150‑1 (MQ=255)
                  tGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCGGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTGACtttc                                                                                                                                   <  3:575197/150‑4 (MQ=255)
                      gTGAAAAATGCGATGCGCAGATGGGGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAAGATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTc                                                                                                                               >  1:1448239/1‑150 (MQ=255)
                      gTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTc                                                                                                                               >  3:1059956/1‑150 (MQ=255)
                       tGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCtttt                                                                                                                                          >  1:261274/1‑138 (MQ=255)
                       tGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCtttt                                                                                                                                          <  2:261274/138‑1 (MQ=255)
                       tGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCa                                                                                                                              >  1:110338/1‑150 (MQ=255)
                        gAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGa                                                                                                                                        <  4:420304/139‑1 (MQ=255)
                         aaaaaTGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAAc                                                                                                                            >  2:857875/1‑150 (MQ=255)
                         aaaaaTGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAAc                                                                                                                            >  3:132398/1‑150 (MQ=255)
                          aaaaTGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACa                                                                                                                           <  2:675751/150‑1 (MQ=255)
                          aaaaTGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACa                                                                                                                           <  2:1779538/150‑1 (MQ=255)
                             atgcgatgcgCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGt                                                                                                                        <  3:485893/150‑1 (MQ=255)
                             atgcgatgcgCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGt                                                                                                                        <  3:485909/150‑1 (MQ=255)
                               gcgatgcgCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTTGTTAAGGTGGTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTaa                                                                                                                      <  1:1505710/150‑1 (MQ=255)
                                      gcAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGa                                                                                                               <  4:1059956/150‑1 (MQ=255)
                                                ttATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACt                                                                                                                                <  1:417776/123‑1 (MQ=255)
                                                ttATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACt                                                                                                                                >  2:417776/1‑123 (MQ=255)
                                                ttATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCtt                                                                                                     <  1:1483493/150‑1 (MQ=255)
                                                 tATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCCCCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTc                                                                                                    >  3:86553/1‑150 (MQ=255)
                                                   tGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAg                                                                                                  <  1:200312/150‑1 (MQ=255)
                                                      gCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                                                                         >  1:1615405/1‑108 (MQ=255)
                                                              aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGgcgc                                                                                                                 >  4:921883/1‑124 (MQ=255)
                                                               gCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGa                                                                                      >  2:1097974/1‑150 (MQ=255)
                                                                ccAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtttt                                                                                                                                          >  1:1195525/1‑97 (MQ=255)
                                                                ccAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCtttt                                                                                                                                          <  2:1195525/97‑1 (MQ=255)
                                                                  aGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACa                                                                                                                           >  1:693547/1‑110 (MQ=255)
                                                                  aGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACa                                                                                                                           <  2:693547/110‑1 (MQ=255)
                                                                  aGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACAt                                                                                   <  1:1000720/150‑1 (MQ=255)
                                                                   ggCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATa                                                                                  <  3:952303/150‑1 (MQ=255)
                                                                      aGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACtt                                                                               <  2:50421/150‑1 (MQ=255)
                                                                        tGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTaa                                                                             >  2:1546460/1‑150 (MQ=255)
                                                                          ggATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaat                                                                           <  4:132398/150‑1 (MQ=255)
                                                                           gATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaata                                                                          >  4:410448/1‑150 (MQ=255)
                                                                             ttGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGttt                                                                                                          <  1:691968/116‑1 (MQ=255)
                                                                              tgtgGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATg                                                                       <  1:1459660/150‑1 (MQ=255)
                                                                                  gtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCGACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTc                                                                   >  3:978825/1‑150 (MQ=255)
                                                                                  gtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTc                                                                   >  2:1488767/1‑150 (MQ=255)
                                                                                     aaTATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGttt                                                                                                          <  1:41252/108‑1 (MQ=255)
                                                                                     aaTATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGttt                                                                                                          >  2:41252/1‑108 (MQ=255)
                                                                                         tCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTc                                                                                                                               >  4:180113/1‑83 (MQ=255)
                                                                                         tCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTc                                                                                                                               <  3:180113/83‑1 (MQ=255)
                                                                                         tCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCt                                                            <  2:1448239/150‑1 (MQ=255)
                                                                                         tCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCt                                                            >  1:525920/1‑150 (MQ=255)
                                                                                          ccTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTAc                                                                                                                                 >  1:476691/1‑80 (MQ=255)
                                                                                           cTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCttt                                                          >  2:1449137/1‑150 (MQ=255)
                                                                                              tttAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCtt                                                                                                     >  1:366258/1‑104 (MQ=255)
                                                                                              tttAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCtt                                                                                                     <  2:366258/104‑1 (MQ=255)
                                                                                              tttAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTtata                                                       >  3:891751/1‑150 (MQ=255)
                                                                                                tAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGCCTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGc                                                     >  2:1556728/1‑150 (MQ=255)
                                                                                                 aaaaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCtttt                                                                                                                                          <  2:585288/64‑1 (MQ=37)
                                                                                                 aaaaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCtttt                                                                                                                                          >  1:585288/1‑64 (MQ=37)
                                                                                                   aaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGtt                                                  >  2:495494/1‑150 (MQ=255)
                                                                                                     ttAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACttt                                                                                                                                    >  2:223714/1‑64 (MQ=255)
                                                                                                     ttAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACttt                                                                                                                                    <  1:223714/66‑3 (MQ=255)
                                                                                                          ggCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAAt                                           <  2:525920/150‑1 (MQ=255)
                                                                                                          ggCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAAt                                           <  2:1729938/150‑1 (MQ=255)
                                                                                                            caccacCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTg                                         <  4:690879/150‑1 (MQ=255)
                                                                                                             accaccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTgg                                        <  1:319803/150‑1 (MQ=255)
                                                                                                               caccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGtt                                      <  3:970019/150‑1 (MQ=255)
                                                                                                               caccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGtt                                      <  3:79849/150‑1 (MQ=255)
                                                                                                                     aaTTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaataat                                                                        >  3:680813/1‑110 (MQ=255)
                                                                                                                     aaTTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACg                                <  4:86553/150‑1 (MQ=255)
                                                                                                                        tCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCttt                                                          <  2:1031634/121‑1 (MQ=255)
                                                                                                                        tCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCttt                                                          >  1:1031634/1‑121 (MQ=255)
                                                                                                                         cAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTAc                            <  4:283858/150‑1 (MQ=255)
                                                                                                                             aaGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACtttt                        >  4:240874/1‑150 (MQ=255)
                                                                                                                                cAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGcc                                               <  1:822009/124‑1 (MQ=255)
                                                                                                                                cAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGcc                                               >  2:822009/1‑124 (MQ=255)
                                                                                                                                 taataGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCa                    >  4:89484/5‑150 (MQ=255)
                                                                                                                                  gAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCAc                   <  3:741107/150‑1 (MQ=255)
                                                                                                                                  gAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCAc                   <  1:1097974/150‑1 (MQ=255)
                                                                                                                                   aaaGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTa                             >  2:1402776/1‑139 (MQ=255)
                                                                                                                                      gTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCcac               <  2:110338/150‑1 (MQ=255)
                                                                                                                                       tCAAAAGCCTCCGGTCGGAGGCTTTTGACTTTTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCcaca              <  2:475686/150‑1 (MQ=255)
                                                                                                                                       tCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCcaca              >  3:31488/1‑150 (MQ=255)
                                                                                                                                         aaaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACtt                                                                               >  1:1747310/1‑83 (MQ=255)
                                                                                                                                         aaaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACtt                                                                               <  2:1747005/83‑1 (MQ=255)
                                                                                                                                         aaaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCt                                                                  >  2:267855/1‑96 (MQ=255)
                                                                                                                                         aaaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACt            >  1:917408/1‑150 (MQ=255)
                                                                                                                                          aaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGt                                                                >  3:280130/1‑97 (MQ=255)
                                                                                                                                          aaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACtt           >  2:1497480/1‑150 (MQ=255)
                                                                                                                                          aaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACtt           >  4:307477/1‑150 (MQ=255)
                                                                                                                                              ccTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAat       <  2:1605128/150‑1 (MQ=255)
                                                                                                                                              ccTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAat       >  4:562230/1‑150 (MQ=255)
                                                                                                                                               cTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAata      <  3:689716/150‑1 (MQ=255)
                                                                                                                                                   gACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTaa  >  1:1397929/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGTTGCCCAATTGGTTCAGACGCTACTTTTTTCACCCACACTTCAATATTAA  >  NC_000913_3_pae_pgi/1311682‑1311978

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: