Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,872,359 T→C 9.4% K374R (AAA→AGA)  gyrB ← DNA gyrase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,872,3590TC 9.4% 593.3 / 46.5 202K374R (AAA→AGA) gyrBDNA gyrase, subunit B
Reads supporting (aligned to +/- strand):  ref base T (68/114);  new base C (11/8);  total (79/123)
Fisher's exact test for biased strand distribution p-value = 8.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

Only 100 of 204 total aligned reads displayed.

ACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACT  >  CP009273/3872212‑3872502
                                                                                                                                                   |                                                                                                                                               
 ccaGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAAtttt                                                                                                                                                >  1:469833/1‑148 (MQ=255)
    ggTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGcc                                                                                                                                             >  1:2653201/1‑148 (MQ=255)
       aCAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAAc                                                                                                                                          <  1:1887265/148‑1 (MQ=255)
        cAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAAcc                                                                                                                                         <  2:3196950/148‑1 (MQ=255)
        cAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAAcc                                                                                                                                         >  2:1452508/1‑148 (MQ=255)
           ttCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACg                                                                                                                                      >  2:3171885/1‑148 (MQ=255)
           ttCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACg                                                                                                                                      >  2:1090741/1‑148 (MQ=255)
            tCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGa                                                                                                                                     <  2:2055198/148‑1 (MQ=255)
              ggAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGAtt                                                                                                                                   <  1:1288904/148‑1 (MQ=255)
               gAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGAttt                                                                                                                                  <  1:2000055/148‑1 (MQ=255)
                aaaGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGAtttt                                                                                                                                 <  2:1726942/148‑1 (MQ=255)
                   gcgcCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTCTGCCAACCACGATTTTcgc                                                                                                                              >  1:2897602/1‑148 (MQ=255)
                    cgcCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCGCCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTcgcg                                                                                                                             >  1:1853640/1‑148 (MQ=255)
                       cGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTCTGCCAACCACGATTTTCGCGTCg                                                                                                                          <  1:665316/148‑1 (MQ=255)
                        ggATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCgg                                                                                                                         <  2:3259527/148‑1 (MQ=255)
                          aTCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGtt                                                                                                                       <  2:1331323/148‑1 (MQ=255)
                           tCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTg                                                                                                                      <  2:1863561/148‑1 (MQ=255)
                           tCGCGGTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTg                                                                                                                      <  2:3166942/148‑1 (MQ=255)
                            cgcgTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTCTGCCAACCACGATTTTCGCGTCGGTTgg                                                                                                                     >  1:2193948/1‑148 (MQ=255)
                            cgcgTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTCTGCCAACCACGATTTTCGCGTCGGTTgg                                                                                                                     >  1:2247687/1‑148 (MQ=255)
                             gcgTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTTCCAACCACGATTTTCGCGTCGGTTggg                                                                                                                    <  2:681549/148‑1 (MQ=255)
                             gcgGTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTTCCAACCACGATTTTCGCGTCGGTTggg                                                                                                                    <  2:680813/148‑1 (MQ=255)
                                 tCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGtttt                                                                                                                >  1:1235438/1‑148 (MQ=255)
                                  ccTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTc                                                                                                               <  1:2323259/148‑1 (MQ=255)
                                    tGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCca                                                                                                             <  2:2442374/148‑1 (MQ=255)
                                    tGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCca                                                                                                             >  2:2771488/1‑148 (MQ=255)
                                     ggCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCcag                                                                                                            <  2:2756764/148‑1 (MQ=255)
                                     ggCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCcag                                                                                                            >  2:1938980/1‑148 (MQ=255)
                                      gCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCcagc                                                                                                           >  1:1710998/1‑148 (MQ=255)
                                      gCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCcagc                                                                                                           >  1:539902/1‑148 (MQ=255)
                                      gCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCcagc                                                                                                           >  2:1253874/1‑148 (MQ=255)
                                      gCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTCTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCcagc                                                                                                           >  1:2824997/1‑148 (MQ=255)
                                      gCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTCTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCcagc                                                                                                           >  1:2825100/1‑148 (MQ=255)
                                       cAGTCTGCCAGTGTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCcagca                                                                                                          <  2:3260930/148‑1 (MQ=255)
                                         gTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAgg                                                                                                        >  2:770459/1‑148 (MQ=255)
                                          tCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGt                                                                                                       >  1:1714481/1‑148 (MQ=255)
                                             gCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTAtt                                                                                                    <  1:2720925/148‑1 (MQ=255)
                                             gCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTAtt                                                                                                    <  2:2950062/148‑1 (MQ=255)
                                              ccAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTc                                                                                                   >  1:494445/1‑148 (MQ=255)
                                              ccAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTCTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTc                                                                                                   <  2:2404004/148‑1 (MQ=255)
                                                aGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTg                                                                                                 >  1:2744023/1‑148 (MQ=255)
                                                aGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTg                                                                                                 >  2:45036/1‑148 (MQ=255)
                                                aGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTg                                                                                                 >  2:46127/1‑148 (MQ=255)
                                                  tttGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGcc                                                                                               >  1:3176660/1‑148 (MQ=255)
                                                  tttGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGcc                                                                                               >  1:3032519/1‑148 (MQ=255)
                                                     gCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCcagc                                                                                            <  2:1392274/148‑1 (MQ=255)
                                                       ccGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCcagcag                                                                                          <  1:2537947/148‑1 (MQ=255)
                                                         ggCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAgtt                                                                                        <  1:2344224/148‑1 (MQ=255)
                                                            aGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAgttcgt                                                                                     <  2:1979387/148‑1 (MQ=255)
                                                             ggCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAgttcgtt                                                                                    >  2:1700838/1‑148 (MQ=255)
                                                             ggCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTCTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAgttcgtt                                                                                    >  2:2308680/1‑148 (MQ=255)
                                                             ggCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTCTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAgttcgtt                                                                                    >  2:1775205/1‑148 (MQ=255)
                                                              gCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAgttcgttc                                                                                   >  1:2132288/1‑148 (MQ=255)
                                                              gCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAgttcgttc                                                                                   >  1:1551288/1‑148 (MQ=255)
                                                               cccGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCa                                                                                  >  2:800080/1‑148 (MQ=255)
                                                                ccGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCAt                                                                                 <  2:747549/148‑1 (MQ=255)
                                                                 cGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATc                                                                                <  1:1005821/148‑1 (MQ=255)
                                                                   cTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATctg                                                                              <  2:2968001/148‑1 (MQ=255)
                                                                      aGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATctgctg                                                                           <  1:1111346/148‑1 (MQ=255)
                                                                      aGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATctgctg                                                                           <  1:769740/148‑1 (MQ=255)
                                                                        tCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGtt                                                                         <  2:1853629/148‑1 (MQ=255)
                                                                        tCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGtt                                                                         <  2:1853669/148‑1 (MQ=255)
                                                                         cGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTc                                                                        <  2:235292/148‑1 (MQ=255)
                                                                           aGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCaa                                                                      >  2:2193263/1‑148 (MQ=255)
                                                                               cACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAAccgc                                                                  <  1:1422264/148‑1 (MQ=255)
                                                                                aCCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAAccgcc                                                                 >  1:2433448/1‑148 (MQ=255)
                                                                                 ccTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAAccgccg                                                                >  2:1055658/1‑148 (MQ=255)
                                                                                  cTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGa                                                               >  2:579798/1‑148 (MQ=255)
                                                                                   tttACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGAt                                                              <  1:535001/148‑1 (MQ=255)
                                                                                       cggcggGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCa                                                          <  2:2269663/148‑1 (MQ=255)
                                                                                        ggcggGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTCTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCAc                                                         <  1:966073/148‑1 (MQ=255)
                                                                                        ggcggGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTCTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCAc                                                         <  1:1817369/148‑1 (MQ=255)
                                                                                              tCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCaga                                                   <  1:1476864/148‑1 (MQ=255)
                                                                                                 tttCACGCGCGCGACGCGCCGCTTCACTGGCACGGGGAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCagaaga                                                <  1:2200770/148‑1 (MQ=255)
                                                                                                  ttCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAgaagaa                                               <  1:808807/148‑1 (MQ=255)
                                                                                                   tCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGaaa                                              <  2:960771/148‑1 (MQ=255)
                                                                                                     aCGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAAcc                                            <  2:5349/148‑1 (MQ=255)
                                                                                                           gcgACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTg                                      <  1:511559/148‑1 (MQ=255)
                                                                                                           gcgACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTg                                      <  2:1421558/148‑1 (MQ=255)
                                                                                                              aCGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCANCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCt                                   <  1:1679083/148‑1 (MQ=255)
                                                                                                              aCGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCt                                   <  2:807404/148‑1 (MQ=255)
                                                                                                              aCGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCt                                   <  2:2653201/148‑1 (MQ=255)
                                                                                                               cgcgCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCtt                                  <  2:287710/148‑1 (MQ=255)
                                                                                                               cgcgCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCtt                                  >  2:1361565/1‑148 (MQ=255)
                                                                                                               cgcgCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTACAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTGtt                                  >  1:578910/1‑148 (MQ=255)
                                                                                                                   ccgcTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGt                              <  2:2748780/148‑1 (MQ=255)
                                                                                                                   ccgcTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTGTGTCTTTGGt                              >  1:1240877/1‑148 (MQ=255)
                                                                                                                      cTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTg                           >  2:3116309/1‑148 (MQ=255)
                                                                                                                             ggcacgggcaGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGgaggag                    <  2:1710998/148‑1 (MQ=255)
                                                                                                                             ggcacgggcaGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGgaggag                    <  2:539902/148‑1 (MQ=255)
                                                                                                                              gcacgggcaGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGgaga                   <  1:3137308/148‑1 (MQ=255)
                                                                                                                              gcacgggcaGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGgaga                   <  2:899283/148‑1 (MQ=255)
                                                                                                                                 cgggcaGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAAtt                >  1:1298101/1‑148 (MQ=255)
                                                                                                                                    gcagcaTCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCg             <  1:2610293/148‑1 (MQ=255)
                                                                                                                                       gcaTCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGt          <  2:2530691/148‑1 (MQ=255)
                                                                                                                                        caTCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTc         <  2:1607107/148‑1 (MQ=255)
                                                                                                                                           cGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCgg      >  1:3265048/1‑148 (MQ=255)
                                                                                                                                            gATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGc     >  1:2756604/1‑148 (MQ=255)
                                                                                                                                             aTAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCa    <  1:3171885/148‑1 (MQ=255)
                                                                                                                                              tAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTGTGGGAGGAGAATTTCGGGTCCGGCAc   >  2:2707477/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
ACCAGGTACAGTTCGGAAAGCGCCGGATCGCGTTCCTGGCAGTCTGCCAGTTTGCCCGGCAGGCCCGCTAAGTCGAGCGCACCTTTACGGCGGGTCATTTCACGCGCGCGACGCGCCGCTTCACGGGCACGGGCAGCATCGATAATTTTGCCAACCACGATTTTCGCGTCGGTTGGGTTTTCCAGCAGGTATTCTGCCAGCAGTTCGTTCATCTGCTGTTCAACCGCCGATTTCACCTCAGAAGAAACCAGTTTGTCTTTGGTCTGGGAGGAGAATTTCGGGTCCGGCACT  >  CP009273/3872212‑3872502

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: