Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,870,977 T→A 5.4% intergenic (‑152/+88) yidB ← / ← gyrB DUF937 family protein/DNA gyrase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,870,9770TA 5.4% 917.9 / 8.7 332intergenic (‑152/+88)yidB/gyrBDUF937 family protein/DNA gyrase, subunit B
Reads supporting (aligned to +/- strand):  ref base T (182/125);  new base A (7/17);  total (190/142)
Fisher's exact test for biased strand distribution p-value = 5.06e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 344 total aligned reads displayed.

TAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAA  >  CP009273/3870828‑3871126
                                                                                                                                                     |                                                                                                                                                     
tAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                                                                                       <  1:235844/150‑1 (MQ=255)
 aaTCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                                      >  2:3809841/1‑150 (MQ=255)
   tCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACa                                                                                                                                                    >  1:1870759/1‑150 (MQ=255)
      ccTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACaaaa                                                                                                                                                 <  2:2715594/150‑1 (MQ=255)
        tGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATc                                                                                                                                               <  1:125630/150‑1 (MQ=255)
           aaaGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATg                                                                                                                                            >  1:613834/1‑150 (MQ=255)
             aGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGaa                                                                                                                                          <  2:2520801/150‑1 (MQ=255)
                gTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAAt                                                                                                                                       >  1:4292692/1‑150 (MQ=255)
                gTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAAt                                                                                                                                       >  1:4590108/1‑150 (MQ=255)
                gTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAAt                                                                                                                                       <  1:3525402/150‑1 (MQ=255)
                gTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCCACAAAATCATGAAAAt                                                                                                                                       >  1:4325892/1‑150 (MQ=255)
                 ttACGTGTTGACCCAAAAATTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAAtt                                                                                                                                      <  1:2576218/150‑1 (MQ=255)
                 ttACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAAtt                                                                                                                                      >  2:4246465/1‑150 (MQ=255)
                 ttACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAAtt                                                                                                                                      >  1:3861388/1‑150 (MQ=255)
                 ttACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTGCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAAtg                                                                                                                                      >  2:4136738/1‑149 (MQ=255)
                  tACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAAtta                                                                                                                                     >  1:3664207/1‑149 (MQ=255)
                  tACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTc                                                                                                                                     >  1:1024119/1‑150 (MQ=255)
                          gACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACAtt                                                                                                                             <  1:3235320/150‑1 (MQ=255)
                          gACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACAtt                                                                                                                             >  1:2927615/1‑150 (MQ=255)
                          gACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACAtt                                                                                                                             >  1:2836691/1‑150 (MQ=255)
                           aCCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTg                                                                                                                            <  2:896668/150‑1 (MQ=255)
                             ccAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGGAAATTCAATACATTGCa                                                                                                                          <  1:867885/150‑1 (MQ=255)
                               aaaaaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAg                                                                                                                        <  2:4466822/150‑1 (MQ=255)
                               aaaaaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAg                                                                                                                        >  2:1142372/1‑150 (MQ=255)
                                 aaaGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAt                                                                                                                      <  1:4421256/150‑1 (MQ=255)
                                    gTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAtttt                                                                                                                   <  2:3069384/150‑1 (MQ=255)
                                           tttGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGc                                                                                                            >  2:1505449/1‑150 (MQ=255)
                                             tGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCt                                                                                                          >  1:5111285/1‑150 (MQ=255)
                                              gCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTg                                                                                                         <  1:3029468/150‑1 (MQ=255)
                                               cGGGATGATAATTGCGGATTGCAGAAATAACAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGa                                                                                                        <  1:4979436/150‑1 (MQ=255)
                                               cGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGa                                                                                                        <  1:3426573/150‑1 (MQ=255)
                                               cGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGa                                                                                                        <  1:222898/150‑1 (MQ=255)
                                                gggATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGAt                                                                                                       <  1:2284440/150‑1 (MQ=255)
                                                 ggATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATa                                                                                                      >  2:2772744/1‑150 (MQ=255)
                                                   atgatAATTGCGGATTGCAGAAATAAAATGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAg                                                                                                    >  2:2388810/1‑150 (MQ=255)
                                                     gatAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCg                                                                                                  <  1:1774520/150‑1 (MQ=255)
                                                      atAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGt                                                                                                 <  2:1191991/150‑1 (MQ=255)
                                                      atAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGt                                                                                                 <  2:438442/150‑1 (MQ=255)
                                                       tAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAATTTCAATACATTGCTATATTTTCGTAGGCCTGATAAGCGTa                                                                                                >  2:3133998/1‑150 (MQ=255)
                                                        aaTTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAg                                                                                               <  1:4055157/150‑1 (MQ=255)
                                                         aTTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAgc                                                                                              <  1:1194593/150‑1 (MQ=255)
                                                         aTTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAgc                                                                                              <  1:325529/150‑1 (MQ=255)
                                                              ggATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATc                                                                                         >  2:1794560/1‑150 (MQ=255)
                                                              ggATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCGTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAATATTTTCGTAGGCCTGATAATCGTAGCGCATc                                                                                         >  2:4973748/1‑150 (MQ=255)
                                                              ggATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCGTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAATATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATc                                                                                         >  2:3558066/1‑150 (MQ=255)
                                                               gATTGCAGAAATAAACAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCa                                                                                        <  1:4784217/150‑1 (MQ=255)
                                                                  tGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                                                     >  2:1070679/1‑150 (MQ=255)
                                                                    cAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAATATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAc                                                                                   <  2:762314/150‑1 (MQ=255)
                                                                    cAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAc                                                                                   >  1:2082415/1‑150 (MQ=255)
                                                                    cAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGAAGCGGCGAGAACGCCAAATCCGGCCTACAAAATAAAGAAAATTCAAAAAATTGAAACATTTTAGAAGTCATGATAATCGTATCCCATCAGGCaa                                                                                   >  1:2518135/1‑149 (MQ=255)
                                                                     aGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTTGGCCTGATAAGCGTAGCGCATCAGGCACg                                                                                  >  1:3741374/1‑150 (MQ=255)
                                                                       aaaTAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCt                                                                                >  1:3297669/1‑150 (MQ=255)
                                                                       aaaTAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCt                                                                                <  2:2998689/150‑1 (MQ=255)
                                                                         aTAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCCTCAGGCACGCTCg                                                                              >  1:4726490/1‑150 (MQ=255)
                                                                            aaaaGGGCGGTTATGATCTCCGCCCTTTTCTTATAGCTTGTTGCCGGATTGTGGGTGAACGCCTTATCCGGCATAAAAAATCTTTAAAATTAAATACATTGAAAGATTTTCGTAGGCCTGATAAGCTTAGCGCATCAGGCACGCTCGCAt                                                                           <  2:324077/150‑1 (MQ=21)
                                                                              aaGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACCTTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCAGCAGGCACGCTCGCATgg                                                                         >  1:738815/1‑150 (MQ=255)
                                                                               tggggggAGTTTATCCCCGCCCTTTTCTTATAGCTTCTTGCGGGTTGCGGCGTGACCGCCTTACCCGGCCTACAAAATCATCAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGt                                                                        <  2:2480226/149‑1 (MQ=21)
                                                                               aGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGt                                                                        >  2:2125981/1‑150 (MQ=255)
                                                                                      agaTGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAATCGTAGCGCATCAGGCACGCTCGCAAGGTTAGCGcc                                                                 >  1:3350843/1‑150 (MQ=255)
                                                                                        atgatCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATACGGCCTACAAAATCATGAAAATTCAATACATTGCAATATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCAt                                                               >  1:4540140/1‑150 (MQ=255)
                                                                                            tctcCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGGGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGAGTGTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTaaa                                                           >  2:5078124/1‑150 (MQ=255)
                                                                                             ctcCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAat                                                          >  1:3956159/1‑150 (MQ=255)
                                                                                             ctcCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCATATCCGGCCAACAAAAAAATGAAAATTCAATAAAATTCAAGATTTTAGTTGGCCTGATAATCGGAGCGCAAAAGGCACGCACGCCAGGTTAGCGCCCATAAat                                                          >  1:2697314/1‑150 (MQ=18)
                                                                                               ccgcatttttATTATAGCTTCTTGCCGGTTGGGGCGTCAACGCATTACCCGGAAAAAAAAAAAATTAAAATTCAATACATTGAATTTTTTTCGTAGGCTTTATATGCGTAGCTAATCTGGCACGCTCGCATGGTTAGCGCCATTAAatat                                                        <  1:516790/144‑1 (MQ=255)
                                                                                                cGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATc                                                       >  2:1575831/1‑150 (MQ=255)
                                                                                                 gCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAATCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCg                                                      >  2:3297889/1‑150 (MQ=255)
                                                                                                   ccTTTTCTTATAGCTTCTTGCCGGTTGCGGCGTGAACGCCTTACCCGGCCTACAAAATCATGAAAATTAAATACATTGAAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGat                                                    <  2:3350322/150‑1 (MQ=255)
                                                                                                   ccTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAATATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGat                                                    <  2:787735/150‑1 (MQ=255)
                                                                                                   ccTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGat                                                    >  1:677933/1‑150 (MQ=255)
                                                                                                   ccTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCATACAAAATCATGAAAATTCAGTACATTGCAATATTTTCTTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGat                                                    <  1:2095365/150‑1 (MQ=255)
                                                                                                   ccTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATACGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTAGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGat                                                    >  1:4206883/1‑150 (MQ=255)
                                                                                                      tttCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATAtt                                                 <  2:5048084/150‑1 (MQ=255)
                                                                                                      tttCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATAtt                                                 >  2:1352068/1‑150 (MQ=255)
                                                                                                       ttCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTc                                                <  2:3934277/150‑1 (MQ=255)
                                                                                                        tCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTcg                                               >  1:1551982/1‑150 (MQ=255)
                                                                                                           tataGCTTCTGGCCGGATGCGGCGTGAACGCCTTATCCGGCAAACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTcgccg                                            <  1:846093/150‑1 (MQ=255)
                                                                                                            ataGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAAGATCGATATTcgccgc                                           >  2:2226708/1‑150 (MQ=255)
                                                                                                            ataGCTTCTTGCCGGATGCGGCGTGAACGCCTTATACGGCCTACAAAATCATGAAAATTCAATACAATGCAAGATTTTCGTAGGCCAGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTcgccgc                                           >  1:1939194/1‑150 (MQ=255)
                                                                                                            agaGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGAAAATTCAATACATTGCATGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTcgccgc                                           <  1:1745785/148‑1 (MQ=255)
                                                                                                             taGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCt                                          >  2:1735466/1‑150 (MQ=255)
                                                                                                              cgCTTCGGGCCGGATGCGGCGTGAACGCCTTAGCCGGCCTACCAAATCATCAAAATGAAATACATTCCAATATTTTCGTTGGCCTTATATGCGTAGCGAATCTGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCtt                                         <  1:4136738/149‑1 (MQ=25)
                                                                                                                cttcttGCCGGATGCGGCGAGAACGCCTTATCCGGCCTACAAAATCCAGAACATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTc                                       <  2:5288327/150‑1 (MQ=255)
                                                                                                                 ttctgGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCTGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCa                                      <  2:5143285/150‑1 (MQ=255)
                                                                                                                   cttGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAgg                                    >  1:3554875/1‑150 (MQ=255)
                                                                                                                       ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGt                                >  2:2014738/1‑150 (MQ=255)
                                                                                                                       ccGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGt                                >  1:968192/1‑150 (MQ=255)
                                                                                                                         ggATGCGGCGTGAACGCCTTATCCGGCCTACAAAATAATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTtc                              >  2:4736163/1‑150 (MQ=255)
                                                                                                                                   tGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATaa                    >  2:3493637/1‑150 (MQ=255)
                                                                                                                                      aCGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACg                 >  1:4003041/1‑150 (MQ=255)
                                                                                                                                        gCCTTATCCGGCCTACAAAATCATTAAAATTCAATACATTGCAAGATTTTCGTAGGCCTTATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCa               <  2:86381/150‑1 (MQ=255)
                                                                                                                                         ccTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCAc              >  2:3392046/1‑150 (MQ=255)
                                                                                                                                         ccTTATCCGGCATAAAAAATAATGAAAATTCAATACATTGCACTATTTTCGTTGGCCTGATAAGCTTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATAGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCAc              <  2:4602412/150‑1 (MQ=255)
                                                                                                                                           ttATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAATATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACgg            >  1:1405271/1‑150 (MQ=255)
                                                                                                                                           ttATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACgg            >  1:3114841/1‑150 (MQ=255)
                                                                                                                                                 ggCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGt      <  2:176647/150‑1 (MQ=255)
                                                                                                                                                    cTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCa   >  1:1726017/1‑150 (MQ=255)
                                                                                                                                                    cTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCa   >  1:2198910/1‑150 (MQ=255)
                                                                                                                                                    cTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCa   >  2:1925488/1‑150 (MQ=255)
                                                                                                                                                     tACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCaa  >  2:3663682/1‑150 (MQ=255)
                                                                                                                                                     tACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCaa  >  2:1428408/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAATCACCTGTAAAGCGTTACGTGTTGACCCAAAAAGTATAGATTTGCGGGATGATAATTGCGGATTGCAGAAATAAAAAGGGCGGAGATGATCTCCGCCCTTTTCTTATAGCTTCTTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAATCATGAAAATTCAATACATTGCAAGATTTTCGTAGGCCTGATAAGCGTAGCGCATCAGGCACGCTCGCATGGTTAGCGCCATTAAATATCGATATTCGCCGCTTTCAGGGCGTTCTCTTCAATAAACGCACGGCGCGGTTCAA  >  CP009273/3870828‑3871126

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