Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,122,248 A→T 46.3% 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,2480AT46.3% 69.7 / inf 272intergenic (+53/‑296)metL/metFBifunctional aspartokinase/homoserine dehydrogenase 2/5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base A (55/91);  new base T (55/71);  total (110/162)
Fisher's exact test for biased strand distribution p-value = 3.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-01

Only 100 of 276 total aligned reads displayed.

CGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTT  >  CP009273/4122099‑4122397
                                                                                                                                                     |                                                                                                                                                     
cgATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTa                                                                                                                                                       <  1:1596741/150‑1 (MQ=255)
  aTAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAAt                                                                                                                                                     >  2:3629141/1‑150 (MQ=255)
   tAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATg                                                                                                                                                    <  1:166801/150‑1 (MQ=255)
    aaCCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGa                                                                                                                                                   <  1:3244485/150‑1 (MQ=255)
        ctctGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAata                                                                                                                                               <  2:1038646/150‑1 (MQ=255)
        ctctGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAata                                                                                                                                               >  1:3955648/1‑150 (MQ=255)
         tctGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAatat                                                                                                                                              <  2:3280121/150‑1 (MQ=255)
          ctGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAatata                                                                                                                                             >  1:696725/1‑150 (MQ=255)
          ctGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAatata                                                                                                                                             >  1:164315/1‑150 (MQ=255)
              tGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGa                                                                                                                                         >  1:2333227/1‑150 (MQ=255)
                   cgcgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTAt                                                                                                                                    <  2:3827016/150‑1 (MQ=255)
                    gcgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATc                                                                                                                                   <  1:2523765/150‑1 (MQ=255)
                    gcgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATc                                                                                                                                   <  1:2249729/150‑1 (MQ=255)
                     cgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCa                                                                                                                                  >  2:3128817/1‑150 (MQ=255)
                     cgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCa                                                                                                                                  <  2:2340105/150‑1 (MQ=255)
                      ggACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTATAATGAAATATATTGAATTATCAt                                                                                                                                 >  1:1990970/1‑150 (MQ=255)
                       gACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATa                                                                                                                                <  1:1579394/150‑1 (MQ=255)
                        aCCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAg                                                                                                                               <  1:3619540/150‑1 (MQ=255)
                         ccTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAgg                                                                                                                              <  1:3339068/150‑1 (MQ=255)
                          cTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGa                                                                                                                             >  1:2434777/1‑150 (MQ=255)
                           tGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATTCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGAt                                                                                                                            >  1:2947661/1‑150 (MQ=255)
                           tGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGAt                                                                                                                            >  2:1928772/1‑150 (MQ=255)
                           tGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGAt                                                                                                                            <  2:3394646/150‑1 (MQ=255)
                             gcgcTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTa                                                                                                                          <  2:3295339/150‑1 (MQ=255)
                               gcTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAgg                                                                                                                        <  1:2266099/150‑1 (MQ=255)
                                 tGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGcc                                                                                                                      <  1:4077853/150‑1 (MQ=255)
                                 tGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGcc                                                                                                                      <  2:1560526/150‑1 (MQ=255)
                                 tGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGcc                                                                                                                      <  1:1890481/150‑1 (MQ=255)
                                    gcgcGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGa                                                                                                                   <  2:629085/150‑1 (MQ=255)
                                        gACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTaa                                                                                                               >  2:1440655/1‑150 (MQ=255)
                                         aCGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAg                                                                                                              <  2:275943/150‑1 (MQ=255)
                                            tCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGt                                                                                                           >  2:636943/1‑150 (MQ=255)
                                                   cgGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTacgac                                                                                                    <  2:1524643/150‑1 (MQ=255)
                                                   cgGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTacgac                                                                                                    <  1:2091033/150‑1 (MQ=255)
                                                    ggggCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTacgacg                                                                                                   >  1:3017687/1‑150 (MQ=255)
                                                         gATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATcc                                                                                              <  2:2515253/150‑1 (MQ=255)
                                                         gATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATcc                                                                                              >  2:1123172/1‑150 (MQ=255)
                                                             cAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCa                                                                                          >  2:2389614/1‑150 (MQ=255)
                                                              aGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATTAAAAATATTGAATAATAATATGATTAGGCCGGATTAATCGTAAACGACCAATCCGCCaa                                                                                         >  1:4309283/1‑150 (MQ=255)
                                                                   gATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGc                                                                                    <  2:52357/150‑1 (MQ=255)
                                                                    atatCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCa                                                                                   >  2:3853153/1‑150 (MQ=255)
                                                                      atCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAAt                                                                                 <  2:216626/150‑1 (MQ=255)
                                                                       tCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATa                                                                                <  1:1086389/150‑1 (MQ=255)
                                                                        cAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATaa                                                                               <  1:3869097/150‑1 (MQ=255)
                                                                        cAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATaa                                                                               <  2:3831141/150‑1 (MQ=255)
                                                                         aaCCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAg                                                                              <  1:1057601/150‑1 (MQ=255)
                                                                         aaCCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAg                                                                              >  1:3468367/1‑150 (MQ=255)
                                                                            cGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTAc                                                                           <  2:3319104/150‑1 (MQ=255)
                                                                              gCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACAt                                                                         <  2:1863648/150‑1 (MQ=255)
                                                                                tGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATgg                                                                       <  2:2821942/150‑1 (MQ=255)
                                                                                tGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATgg                                                                       >  2:2910417/1‑150 (MQ=255)
                                                                                tGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATgg                                                                       <  2:1645498/150‑1 (MQ=255)
                                                                                 ggCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGt                                                                      <  2:572238/150‑1 (MQ=255)
                                                                                   cacaGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTa                                                                    >  1:3744904/1‑150 (MQ=255)
                                                                                      aGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGtt                                                                 >  1:1381755/1‑150 (MQ=255)
                                                                                      aGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCACTAAGTACATGGTGAGtt                                                                 >  2:3934795/1‑150 (MQ=255)
                                                                                        ttgttgTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTta                                                               <  1:743472/150‑1 (MQ=255)
                                                                                            tgtAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATAtt                                                           >  2:2673642/1‑150 (MQ=255)
                                                                                              tAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTg                                                         >  1:2835260/1‑150 (MQ=255)
                                                                                               aaTTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGc                                                        >  1:2040298/1‑150 (MQ=255)
                                                                                                aTTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCa                                                       >  2:2950831/1‑150 (MQ=255)
                                                                                                aTTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCa                                                       >  1:3142157/1‑150 (MQ=255)
                                                                                                   tCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCACTAAGTACATGGTTAGTTTATATTTGCAGTc                                                    >  2:4192028/1‑150 (MQ=255)
                                                                                                   tCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTc                                                    >  2:3739280/1‑150 (MQ=255)
                                                                                                    cAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTcc                                                   >  2:3054414/1‑150 (MQ=255)
                                                                                                     aGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCg                                                  <  1:1112964/150‑1 (MQ=255)
                                                                                                      gAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCgg                                                 >  2:3563342/1‑150 (MQ=255)
                                                                                                        aaTTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtt                                               <  1:3206567/150‑1 (MQ=255)
                                                                                                        aaTTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtt                                               >  2:4266987/1‑150 (MQ=255)
                                                                                                             aataatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgctt                                          <  2:3929081/150‑1 (MQ=255)
                                                                                                              ataatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgcttt                                         <  1:3850425/150‑1 (MQ=255)
                                                                                                                aatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGAGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCAGGtttgctttgc                                       >  1:2982142/1‑150 (MQ=255)
                                                                                                                aatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgctttgc                                       >  1:3565676/1‑150 (MQ=255)
                                                                                                                 atGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCa                                      >  2:1775952/1‑150 (MQ=255)
                                                                                                                  tGCCCGGTACACATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCAt                                     <  2:4354446/150‑1 (MQ=255)
                                                                                                                   gCCCGGTACTCATGTTTTCGGGTTTATTGTTTCTAATGAAATATTTTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGTATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATa                                    <  2:424281/150‑1 (MQ=255)
                                                                                                                   gCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATa                                    >  2:1532245/1‑150 (MQ=255)
                                                                                                                   gCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATa                                    <  2:424008/150‑1 (MQ=255)
                                                                                                                    cccGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATAc                                   <  1:2360593/150‑1 (MQ=255)
                                                                                                                    cccGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATAc                                   <  1:1641170/150‑1 (MQ=255)
                                                                                                                      cGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCg                                 >  1:4242625/1‑150 (MQ=255)
                                                                                                                       ggTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATCTTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCgg                                >  2:3121143/1‑150 (MQ=255)
                                                                                                                          aCTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAtt                             >  1:2144864/1‑150 (MQ=255)
                                                                                                                            tCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAtttt                           <  1:458542/150‑1 (MQ=255)
                                                                                                                              aTGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttct                         >  1:2024207/1‑150 (MQ=255)
                                                                                                                              aTGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttct                         >  1:3947149/1‑150 (MQ=255)
                                                                                                                               tGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctt                        <  1:873003/150‑1 (MQ=255)
                                                                                                                               tGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctt                        <  1:1826755/150‑1 (MQ=255)
                                                                                                                                gTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttcttt                       >  1:962860/1‑150 (MQ=255)
                                                                                                                                gTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttcttt                       >  1:1297632/1‑150 (MQ=255)
                                                                                                                                     cGGGTTTATGGTTTCTAATGACATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATtttctttttctt                  <  2:1737062/150‑1 (MQ=255)
                                                                                                                                      gggTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTa                 <  1:553448/150‑1 (MQ=255)
                                                                                                                                       ggTTTATGGTTTCTAAAGAAATATAGTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTAc                <  2:722521/150‑1 (MQ=255)
                                                                                                                                        gTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTAcc               <  2:2355600/150‑1 (MQ=255)
                                                                                                                                             tGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCt          >  2:1307315/1‑150 (MQ=255)
                                                                                                                                             tGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCt          >  2:3499530/1‑150 (MQ=255)
                                                                                                                                              ggTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTg         >  2:3106511/1‑150 (MQ=255)
                                                                                                                                                 ttCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAg      >  2:2182589/1‑150 (MQ=255)
                                                                                                                                                  tCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGt     >  1:1018630/1‑150 (MQ=255)
                                                                                                                                                     aaTGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGtttt  >  1:581105/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGTTTT  >  CP009273/4122099‑4122397

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: