Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,122,248 A→T 90.7% intergenic (+53/‑296) metL → / → metF Bifunctional aspartokinase/homoserine dehydrogenase 2/5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,122,2480AT90.7% 962.9 / 72.0 301intergenic (+53/‑296)metL/metFBifunctional aspartokinase/homoserine dehydrogenase 2/5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base A (16/12);  new base T (127/145);  total (144/157)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 301 total aligned reads displayed.

ATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGT  >  CP009273/4122101‑4122394
                                                                                                                                                   |                                                                                                                                                  
aTAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTAt                                                                                                                                                  >  1:2571384/1‑150 (MQ=255)
 tAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATg                                                                                                                                                 >  2:212960/1‑150 (MQ=255)
 tAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATg                                                                                                                                                 >  1:2162145/1‑150 (MQ=255)
  aaCCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGa                                                                                                                                                >  2:3587409/1‑150 (MQ=255)
  aaCCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGa                                                                                                                                                >  2:2861507/1‑150 (MQ=255)
       tctGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGACatat                                                                                                                                           <  2:1945682/150‑1 (MQ=255)
         tGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAatatat                                                                                                                                         <  2:274606/150‑1 (MQ=255)
         tGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAatatat                                                                                                                                         >  1:2411168/1‑150 (MQ=255)
         tGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAatatat                                                                                                                                         <  1:996126/150‑1 (MQ=255)
          ggTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATAtt                                                                                                                                        <  2:1214615/150‑1 (MQ=255)
            tGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGa                                                                                                                                      >  1:3168183/1‑150 (MQ=255)
            tGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGa                                                                                                                                      <  1:251142/150‑1 (MQ=255)
              aTCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAAt                                                                                                                                    >  1:712038/1‑150 (MQ=255)
              aTCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAAt                                                                                                                                    <  1:458908/150‑1 (MQ=255)
               tCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAAtt                                                                                                                                   <  2:2642091/150‑1 (MQ=255)
                   cgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCa                                                                                                                               >  1:2604372/1‑150 (MQ=255)
                             gcTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAgg                                                                                                                     <  1:4120099/150‑1 (MQ=255)
                             gcTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAgg                                                                                                                     >  1:2116737/1‑150 (MQ=255)
                             gcTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAgg                                                                                                                     >  1:1466411/1‑150 (MQ=255)
                               tGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGcc                                                                                                                   >  1:512510/1‑150 (MQ=255)
                                gggCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCg                                                                                                                  <  1:4251016/150‑1 (MQ=255)
                                gggCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCAGGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCg                                                                                                                  >  2:1692315/1‑150 (MQ=255)
                                 ggCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCgg                                                                                                                 <  1:3278278/150‑1 (MQ=255)
                                   cgcgACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGAt                                                                                                               >  1:603703/1‑150 (MQ=255)
                                   cgcgACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGAt                                                                                                               <  2:4172176/150‑1 (MQ=255)
                                   cgcgACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGAt                                                                                                               >  1:394648/1‑150 (MQ=255)
                                     cgACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTa                                                                                                             >  1:2899978/1‑150 (MQ=255)
                                          tCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGt                                                                                                        <  1:1995386/150‑1 (MQ=255)
                                           cACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGtt                                                                                                       >  1:2950345/1‑150 (MQ=255)
                                           cACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGtt                                                                                                       >  1:139737/1‑150 (MQ=255)
                                           cACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGtt                                                                                                       <  1:3681053/150‑1 (MQ=255)
                                            aCCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCCCAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGGTTATGGTTTCTTATTAAAAATATTGAATTATCATAGGATTATGACGGATTAAGCGtgt                                                                                                      >  1:3485722/1‑148 (MQ=255)
                                            aCCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGttt                                                                                                      <  1:311141/150‑1 (MQ=255)
                                               gccgGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTacg                                                                                                   <  2:1750532/150‑1 (MQ=255)
                                                 cgGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTacgac                                                                                                 >  1:1007093/1‑150 (MQ=255)
                                                  tgggCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTacgacg                                                                                                >  2:229913/2‑150 (MQ=255)
                                                  tgggCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTacgacg                                                                                                >  2:2545870/2‑150 (MQ=255)
                                                  ggggCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTacgacg                                                                                                >  1:3930493/1‑150 (MQ=255)
                                                    ggCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGaa                                                                                              >  1:1514322/1‑150 (MQ=255)
                                                      cGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATc                                                                                            >  1:2835463/1‑150 (MQ=255)
                                                      agATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATc                                                                                            <  2:1513868/149‑1 (MQ=255)
                                                       gATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATcc                                                                                           >  2:3107690/1‑150 (MQ=255)
                                                       gATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATcc                                                                                           >  2:3107667/1‑150 (MQ=255)
                                                       gATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATcc                                                                                           >  1:4357606/1‑150 (MQ=255)
                                                        aTTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCg                                                                                          <  1:4418073/150‑1 (MQ=255)
                                                        aTTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCg                                                                                          <  2:4392946/150‑1 (MQ=255)
                                                           cAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCa                                                                                       >  1:4508294/1‑150 (MQ=255)
                                                               cGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaagaa                                                                                   >  1:4019951/1‑150 (MQ=255)
                                                                 gATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGc                                                                                 <  2:125479/150‑1 (MQ=255)
                                                                  atatCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCa                                                                                <  1:1282027/150‑1 (MQ=255)
                                                                   tatCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCaa                                                                               <  2:575754/150‑1 (MQ=255)
                                                                    atCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATCATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAAt                                                                              <  1:2743850/150‑1 (MQ=255)
                                                                      cAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATaa                                                                            <  2:4528751/150‑1 (MQ=255)
                                                                       aaCCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAg                                                                           <  2:2638101/150‑1 (MQ=255)
                                                                         ccGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTa                                                                         <  1:2316413/150‑1 (MQ=255)
                                                                             cTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATg                                                                     <  2:2109098/150‑1 (MQ=255)
                                                                             cTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATg                                                                     >  2:3758461/1‑150 (MQ=255)
                                                                             cTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATg                                                                     >  1:3867371/1‑150 (MQ=255)
                                                                              tGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATgg                                                                    <  2:759892/150‑1 (MQ=255)
                                                                              tGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATgg                                                                    >  1:1593544/1‑150 (MQ=255)
                                                                              tGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATgg                                                                    >  2:1443454/1‑150 (MQ=255)
                                                                                gCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGtt                                                                  <  1:3055342/150‑1 (MQ=255)
                                                                                      ttgttgTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTta                                                            <  2:2732368/150‑1 (MQ=255)
                                                                                       tgttgtAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTtat                                                           <  2:4420467/150‑1 (MQ=255)
                                                                                         ttgtAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTatat                                                         >  1:2259106/1‑150 (MQ=255)
                                                                                         ttgtAATTTCAGAAATTTAATAATGCACGGTACTCATGTTTTCGGGTTTATGGTGTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTatat                                                         <  2:1470438/150‑1 (MQ=255)
                                                                                                ttCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGt                                                  <  2:394645/150‑1 (MQ=255)
                                                                                                 tCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTc                                                 <  2:706391/150‑1 (MQ=255)
                                                                                                   aGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCg                                               <  2:3168183/150‑1 (MQ=255)
                                                                                                       aTTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGttt                                           <  2:110790/150‑1 (MQ=255)
                                                                                                        tttAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttg                                          >  1:3124931/1‑150 (MQ=255)
                                                                                                           aataatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgctt                                       <  2:2041240/150‑1 (MQ=255)
                                                                                                            ataatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgcttt                                      >  2:1961154/1‑150 (MQ=255)
                                                                                                             taatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgctttg                                     <  1:1210563/150‑1 (MQ=255)
                                                                                                             taatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgctttg                                     <  2:344255/150‑1 (MQ=255)
                                                                                                              aatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgctttgc                                    >  2:701910/1‑150 (MQ=255)
                                                                                                              aatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgctttgc                                    >  2:45149/1‑150 (MQ=255)
                                                                                                              aatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgctttgc                                    >  2:4313221/1‑150 (MQ=255)
                                                                                                              aatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgctttgc                                    >  2:3247411/1‑150 (MQ=255)
                                                                                                               atGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCa                                   >  1:1131266/1‑150 (MQ=255)
                                                                                                               atGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCa                                   <  1:3424817/150‑1 (MQ=255)
                                                                                                                  cccGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATAc                                <  2:3310489/150‑1 (MQ=255)
                                                                                                                   ccGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATAcc                               <  2:2411168/150‑1 (MQ=255)
                                                                                                                     ggTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCgg                             <  2:3485722/150‑1 (MQ=255)
                                                                                                                      gTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGa                            <  2:3608438/150‑1 (MQ=255)
                                                                                                                         ctactGTTTTCGGGTTTATTGTTTCTTATTAATTATATTGAATTATCATAGGATTAGGCAGGATTAAGCGTTTTCGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCTGTCCGGTTTGCTTTGCATACCGGAttt                         <  1:564730/146‑1 (MQ=255)
                                                                                                                         cTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttt                         <  1:4179727/150‑1 (MQ=255)
                                                                                                                          tCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAtttt                        <  1:833864/150‑1 (MQ=255)
                                                                                                                              gTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttcttt                    <  1:1873470/150‑1 (MQ=255)
                                                                                                                              gTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttcttt                    <  2:2899978/150‑1 (MQ=255)
                                                                                                                                tttCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttcttttt                  >  1:2867917/1‑150 (MQ=255)
                                                                                                                                 ttCGGGTTTATGGTTTCTTATGAACTATATTGAATTATCAGAGGATTAGGCCGGATTAAGCGATTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctttttc                 <  2:3965088/150‑1 (MQ=255)
                                                                                                                                 ttCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctttttc                 <  2:2036214/150‑1 (MQ=255)
                                                                                                                                   cGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATtttctttttctt               <  2:512510/150‑1 (MQ=255)
                                                                                                                                      gTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTAcc            <  1:1506659/150‑1 (MQ=255)
                                                                                                                                       tttATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCa           >  1:3794363/1‑150 (MQ=255)
                                                                                                                                         tattgattCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTTCAGTCCGGTTTGATTTGCATACCGGATTTTCTTTTTCTTACCATc         >  2:1675831/7‑150 (MQ=255)
                                                                                                                                            ggTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTg      <  2:449594/150‑1 (MQ=255)
                                                                                                                                            ggTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTg      >  2:538735/1‑150 (MQ=255)
                                                                                                                                                tCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGt  <  1:3338915/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGT  >  CP009273/4122101‑4122394

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: