Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,175,657 G→T 10.0% E191* (GAA→TAA)  cobB → deacetylase of acs and cheY, chemotaxis regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,175,6570GT10.0% 630.4 / 32.0 190E191* (GAA→TAA) cobBdeacetylase of acs and cheY, chemotaxis regulator
Reads supporting (aligned to +/- strand):  ref base G (83/88);  new base T (5/14);  total (88/102)
Fisher's exact test for biased strand distribution p-value = 8.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 190 total aligned reads displayed.

TTCCGGCACCCTTGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAGCGGGA  >  CP009273/1175511‑1175803
                                                                                                                                                  |                                                                                                                                                  
 tCCGGCACCCTTGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACTaa                                                                                                                                                  <  1:692639/148‑1 (MQ=255)
    ggCACCCTTGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCg                                                                                                                                               <  2:2267754/148‑1 (MQ=255)
     gCACCCTTGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGa                                                                                                                                              >  2:3263849/1‑148 (MQ=255)
     gCACCCTTGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGa                                                                                                                                              >  2:2391230/1‑148 (MQ=255)
      cACCCTTGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGaa                                                                                                                                             >  1:389536/1‑148 (MQ=255)
       aCCCTTGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGaaa                                                                                                                                            >  2:1491746/1‑148 (MQ=255)
           ttGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGc                                                                                                                                        <  1:3030691/148‑1 (MQ=255)
           ttGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGc                                                                                                                                        >  1:2906526/1‑148 (MQ=255)
            tGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCa                                                                                                                                       >  2:2213478/1‑148 (MQ=255)
            tGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCa                                                                                                                                       >  1:3191092/1‑148 (MQ=255)
            tGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCa                                                                                                                                       >  2:423243/1‑148 (MQ=255)
             gcgcCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACTAAGCGAAACTGCAt                                                                                                                                      <  1:1269029/148‑1 (MQ=255)
             gcgcCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACTAAGCGAAACTGCAt                                                                                                                                      <  1:1269355/148‑1 (MQ=255)
             gcgcCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCAt                                                                                                                                      >  1:1347904/1‑148 (MQ=255)
                 ccGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAATCGAAACTGCATGgcg                                                                                                                                  >  2:2287883/1‑148 (MQ=255)
                        tAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCacac                                                                                                                           >  2:2606343/1‑148 (MQ=255)
                          gtgtGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCg                                                                                                                         <  2:2682754/148‑1 (MQ=255)
                             tgGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTgg                                                                                                                      <  1:3333870/148‑1 (MQ=255)
                               gTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGaa                                                                                                                    >  2:2168623/1‑148 (MQ=255)
                                tttGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAAc                                                                                                                   >  1:54299/1‑148 (MQ=255)
                                tttGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAAc                                                                                                                   >  1:2897998/1‑148 (MQ=255)
                                   ggCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGa                                                                                                                <  1:1022125/148‑1 (MQ=255)
                                    gCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGaa                                                                                                               >  2:1341221/1‑148 (MQ=255)
                                     cGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAAt                                                                                                              <  1:872364/148‑1 (MQ=255)
                                     cGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAAt                                                                                                              <  1:1765203/148‑1 (MQ=255)
                                      gAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATc                                                                                                             <  2:1564383/148‑1 (MQ=255)
                                      gAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATc                                                                                                             <  2:3326606/148‑1 (MQ=255)
                                       aaaTGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCt                                                                                                            >  1:2392287/1‑148 (MQ=255)
                                         aTGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTg                                                                                                          >  1:1002231/1‑148 (MQ=255)
                                          tGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGa                                                                                                         >  2:1388913/1‑148 (MQ=255)
                                          tGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGa                                                                                                         >  2:2922358/1‑148 (MQ=255)
                                           gCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACTAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGaa                                                                                                        <  2:3313557/148‑1 (MQ=255)
                                              aCTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCg                                                                                                     <  1:3201754/148‑1 (MQ=255)
                                              aCTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCg                                                                                                     >  2:2433467/1‑148 (MQ=255)
                                               cTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGa                                                                                                    >  2:1633037/1‑148 (MQ=255)
                                                 cGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGt                                                                                                  <  1:1908163/148‑1 (MQ=255)
                                                 cGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGt                                                                                                  <  1:2792922/148‑1 (MQ=255)
                                                  ggCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTc                                                                                                 <  2:3407978/148‑1 (MQ=255)
                                                   gCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCa                                                                                                <  2:2740485/148‑1 (MQ=255)
                                                      tGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGACTCTTGAACCGAGTCAGGt                                                                                             >  2:2459967/1‑148 (MQ=255)
                                                      tGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGt                                                                                             >  2:2459939/1‑148 (MQ=255)
                                                       ggATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGtt                                                                                            <  2:2272900/148‑1 (MQ=255)
                                                         aTGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTgg                                                                                          >  1:2182756/1‑148 (MQ=255)
                                                            aaaTTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACTAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTaa                                                                                       <  2:1061596/148‑1 (MQ=255)
                                                              aTTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATg                                                                                     >  2:2345922/1‑148 (MQ=255)
                                                                ttATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGaa                                                                                   >  1:3012008/1‑148 (MQ=255)
                                                                 tataTGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAAt                                                                                  <  2:2645695/148‑1 (MQ=255)
                                                                 tataTGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAAt                                                                                  <  1:1664657/148‑1 (MQ=255)
                                                                    atGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTg                                                                               >  2:2545536/1‑148 (MQ=255)
                                                                     tGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACTAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGc                                                                              >  2:1454310/1‑148 (MQ=255)
                                                                       gCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCg                                                                            <  2:798416/148‑1 (MQ=255)
                                                                             tCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAat                                                                      <  2:1235970/148‑1 (MQ=255)
                                                                             tCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAat                                                                      <  2:2086056/148‑1 (MQ=255)
                                                                             tCGATGGCCGACATTGTCATTGCCATTGGCACATCCGGGCATGATTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGACGAGAAat                                                                      <  2:2086098/148‑1 (MQ=255)
                                                                                 tGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTa                                                                  <  1:3161944/148‑1 (MQ=255)
                                                                                   gCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACTAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACg                                                                >  1:2370848/1‑148 (MQ=255)
                                                                                     cGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGc                                                              >  1:3286806/1‑148 (MQ=255)
                                                                                     cGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGc                                                              <  1:2449993/148‑1 (MQ=255)
                                                                                      gATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGcc                                                             <  2:3012008/148‑1 (MQ=255)
                                                                                       atatTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGccc                                                            >  2:3059513/1‑148 (MQ=255)
                                                                                       atatTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGccc                                                            >  2:1289020/1‑148 (MQ=255)
                                                                                       atatTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGccc                                                            >  2:1289548/1‑148 (MQ=255)
                                                                                        tatTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCg                                                           >  1:2288118/1‑148 (MQ=255)
                                                                                           tttCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCa                                                        >  2:200247/1‑148 (MQ=255)
                                                                                                 tGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACNCCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAg                                                  >  1:1805611/1‑148 (MQ=255)
                                                                                                 tGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAg                                                  >  2:3489651/1‑148 (MQ=255)
                                                                                                   ccATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAggt                                                >  2:3540334/1‑148 (MQ=255)
                                                                                                   ccATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAggt                                                >  2:663070/1‑148 (MQ=255)
                                                                                                     aTTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAggtgg                                              >  1:1069168/1‑148 (MQ=255)
                                                                                                     aTTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAggtgg                                              >  1:3258033/1‑148 (MQ=255)
                                                                                                      ttGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAggtggt                                             >  1:1268005/1‑148 (MQ=255)
                                                                                                      ttGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAggtggt                                             <  2:2759188/148‑1 (MQ=255)
                                                                                                      ttGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAggtggt                                             <  2:981269/148‑1 (MQ=255)
                                                                                                      ttGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAggtggt                                             >  1:967739/1‑148 (MQ=255)
                                                                                                       tGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGgtggtg                                            >  2:627024/1‑148 (MQ=255)
                                                                                                         gCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGcc                                          <  2:460536/148‑1 (MQ=255)
                                                                                                          cacaTCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGcat                                         <  2:565438/148‑3 (MQ=255)
                                                                                                             aTCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAg                                      >  2:724833/1‑148 (MQ=255)
                                                                                                              tCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGt                                     <  1:907603/148‑1 (MQ=255)
                                                                                                               ccGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGtt                                    >  2:1674290/1‑148 (MQ=255)
                                                                                                               ccGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGtt                                    <  2:3490256/148‑1 (MQ=255)
                                                                                                                   gCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTttgtt                                <  2:1016078/148‑1 (MQ=255)
                                                                                                                       gTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGaaa                            <  2:389536/148‑1 (MQ=255)
                                                                                                                        tttATCCGGCGGCTGGGTTTGTTCACTAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGaaaa                           <  1:1454310/148‑1 (MQ=255)
                                                                                                                        tttATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGaaaa                           >  1:2794913/1‑148 (MQ=255)
                                                                                                                          tATCCGGCGGCTGGGTTTGTTCACTAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGt                         <  1:2226210/148‑1 (MQ=255)
                                                                                                                                   gcTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAgg                >  2:1960937/1‑148 (MQ=255)
                                                                                                                                    cTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAggg               <  1:2345922/148‑1 (MQ=255)
                                                                                                                                    cTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAggg               <  2:103295/148‑1 (MQ=255)
                                                                                                                                    cTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAggg               <  1:1491746/148‑1 (MQ=255)
                                                                                                                                          ttGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTaaa         >  2:3522799/1‑148 (MQ=255)
                                                                                                                                          ttGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTaaa         <  1:3211012/148‑1 (MQ=255)
                                                                                                                                           tGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTaaaa        >  2:1697712/1‑148 (MQ=255)
                                                                                                                                           tGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTaaaa        <  2:907221/148‑1 (MQ=255)
                                                                                                                                            gTTCACTAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAg       <  1:2587984/148‑1 (MQ=255)
                                                                                                                                            gTTCACTAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAg       <  1:2584134/148‑1 (MQ=255)
                                                                                                                                            gTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAg       <  1:3099463/148‑1 (MQ=255)
                                                                                                                                                aCGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAGCggg   >  1:887894/1‑148 (MQ=255)
                                                                                                                                                aCGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAGCggg   >  2:1617787/1‑148 (MQ=255)
                                                                                                                                                 cGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAGCGGGa  >  2:2775427/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTCCGGCACCCTTGCGCCCGCACGTAGTGTGGTTTGGCGAAATGCCACTCGGCATGGATGAAATTTATATGGCGTTGTCGATGGCCGATATTTTCATTGCCATTGGCACATCCGGGCATGTTTATCCGGCGGCTGGGTTTGTTCACGAAGCGAAACTGCATGGCGCGCACACCGTGGAACTGAATCTTGAACCGAGTCAGGTTGGTAATGAATTTGCCGAGAAATATTACGGCCCGGCAAGCCAGGTGGTGCCTGAGTTTGTTGAAAAGTTGCTGAAGGGATTAAAAGCGGGA  >  CP009273/1175511‑1175803

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