Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,122,248 A→T 68.1% 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,2480AT68.1% 395.9 / 263.3 269intergenic (+53/‑296)metL/metFBifunctional aspartokinase/homoserine dehydrogenase 2/5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base A (41/45);  new base T (77/106);  total (118/151)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 271 total aligned reads displayed.

CGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGT  >  CP009273/4122099‑4122394
                                                                                                                                                     |                                                                                                                                                  
  aTAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAAt                                                                                                                                                  <  2:298631/150‑1 (MQ=255)
        ctctGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAata                                                                                                                                            <  1:4038604/150‑1 (MQ=255)
         tctGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAatat                                                                                                                                           <  2:3909081/150‑1 (MQ=255)
           tGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAatatat                                                                                                                                         <  2:2289328/150‑1 (MQ=255)
            ggTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATAtt                                                                                                                                        >  2:4036710/1‑150 (MQ=255)
              tGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGa                                                                                                                                      <  2:2571306/150‑1 (MQ=255)
                  ccGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTa                                                                                                                                  <  2:1014498/150‑1 (MQ=255)
                  ccGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTa                                                                                                                                  <  1:2696626/150‑1 (MQ=255)
                  ccGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTa                                                                                                                                  <  2:2461120/150‑1 (MQ=255)
                  ccGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTa                                                                                                                                  <  1:1826980/150‑1 (MQ=255)
                   cgcgGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTAt                                                                                                                                 <  2:3979894/150‑1 (MQ=255)
                          cTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGa                                                                                                                          >  1:3974085/1‑150 (MQ=255)
                           tGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGAt                                                                                                                         <  2:4166420/150‑1 (MQ=255)
                              cgcTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAg                                                                                                                      <  2:1143590/150‑1 (MQ=255)
                                cTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGc                                                                                                                    >  2:725667/1‑150 (MQ=255)
                                 tGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGcc                                                                                                                   >  1:1024189/1‑150 (MQ=255)
                                  gggCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCg                                                                                                                  >  1:2182279/1‑150 (MQ=255)
                                  gggCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCg                                                                                                                  <  1:407254/150‑1 (MQ=255)
                                     cgcgACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGAt                                                                                                               <  1:775085/150‑1 (MQ=255)
                                     cgcgACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGAt                                                                                                               <  1:530914/150‑1 (MQ=255)
                                     cgcgACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGAt                                                                                                               <  2:3947666/150‑1 (MQ=255)
                                           gTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCg                                                                                                         >  2:4256698/1‑150 (MQ=255)
                                           gTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCg                                                                                                         >  2:2820517/1‑150 (MQ=255)
                                                cgccgGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTTATTTCAGAAATTTAATAATTCCCGGTACTCATGTTTTCGGGTTTATGTTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTac                                                                                                    <  2:3747399/150‑1 (MQ=255)
                                                cgccgGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTac                                                                                                    <  1:3811661/150‑1 (MQ=255)
                                                cgccgGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTac                                                                                                    <  1:821430/150‑1 (MQ=255)
                                                     gggCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTAcgacga                                                                                               >  2:2366195/1‑150 (MQ=255)
                                                         gATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTACAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATcc                                                                                           <  2:3823670/150‑1 (MQ=255)
                                                              aGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaa                                                                                      <  1:1453846/150‑1 (MQ=255)
                                                                tCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaaga                                                                                    <  2:3651308/150‑1 (MQ=255)
                                                                 cGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCTGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaagaa                                                                                   <  1:4109967/150‑1 (MQ=255)
                                                                 cGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaagaa                                                                                   <  1:3737645/150‑1 (MQ=255)
                                                                 cGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaagaa                                                                                   <  1:2157354/150‑1 (MQ=255)
                                                                  ggATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaagaag                                                                                  >  1:3295315/1‑150 (MQ=255)
                                                                  ggATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaagaag                                                                                  >  1:68001/1‑150 (MQ=255)
                                                                  ggATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaagaag                                                                                  <  2:3793733/150‑1 (MQ=255)
                                                                  ggATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCaagaag                                                                                  <  2:4266154/150‑1 (MQ=255)
                                                                   gATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGc                                                                                 >  2:484105/1‑150 (MQ=255)
                                                                     tatCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCaa                                                                               <  1:2083943/150‑1 (MQ=255)
                                                                      atCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATACTGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAAt                                                                              <  2:1793677/150‑1 (MQ=255)
                                                                         aaCCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAg                                                                           <  1:552315/150‑1 (MQ=255)
                                                                         aaCCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAg                                                                           >  2:634265/1‑150 (MQ=255)
                                                                          aCCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGt                                                                          <  2:2436842/150‑1 (MQ=255)
                                                                          aCCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGt                                                                          <  2:886873/150‑1 (MQ=255)
                                                                            cGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATTAGTAc                                                                        <  2:3385714/150‑1 (MQ=255)
                                                                            cGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTAc                                                                        <  2:3922703/150‑1 (MQ=255)
                                                                               cTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATg                                                                     <  1:4245523/150‑1 (MQ=255)
                                                                               cTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATg                                                                     <  1:1996845/150‑1 (MQ=255)
                                                                                  gCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGtt                                                                  <  1:4300773/150‑1 (MQ=255)
                                                                                     caGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGt                                                               >  1:3987693/1‑150 (MQ=255)
                                                                                      aGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGtt                                                              <  1:2899641/150‑1 (MQ=255)
                                                                                       gttgttGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGtgt                                                             >  2:1515149/1‑148 (MQ=255)
                                                                                        ttgttgTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTta                                                            >  2:2740713/1‑150 (MQ=255)
                                                                                        ttgttgTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTta                                                            <  1:484105/150‑1 (MQ=255)
                                                                                          gttgtAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTtata                                                          >  2:1586834/1‑150 (MQ=255)
                                                                                          gttgtAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTtata                                                          >  2:2543205/1‑150 (MQ=255)
                                                                                           ttgtAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTatat                                                         <  1:3707036/150‑1 (MQ=255)
                                                                                           ttgtAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGGTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTatat                                                         <  1:3704360/150‑1 (MQ=255)
                                                                                              tAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTg                                                      <  2:3892084/150‑1 (MQ=255)
                                                                                              tAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTg                                                      >  1:3732419/1‑150 (MQ=255)
                                                                                              tAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTg                                                      >  1:1221414/1‑150 (MQ=255)
                                                                                              tAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTg                                                      >  1:301649/1‑150 (MQ=255)
                                                                                                aTTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCa                                                    >  2:799105/1‑150 (MQ=255)
                                                                                                  ttCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGt                                                  <  1:4228596/150‑1 (MQ=255)
                                                                                                   tCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGAGTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTc                                                 >  2:3091906/1‑150 (MQ=255)
                                                                                                    cAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTcc                                                >  1:599752/1‑150 (MQ=255)
                                                                                                     aGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCg                                               >  1:3336021/1‑150 (MQ=255)
                                                                                                       aaaTTTAATAATGCCCGGTACTCATTTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATTGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGt                                             <  2:2118947/150‑1 (MQ=255)
                                                                                                        aaTTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATAATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGAAGTCCGGtt                                            >  2:3159711/1‑150 (MQ=255)
                                                                                                        aaTTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtt                                            >  2:1317127/1‑150 (MQ=255)
                                                                                                         aTTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGttt                                           <  1:1079931/150‑1 (MQ=255)
                                                                                                              ataatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgcttt                                      >  2:3508336/1‑150 (MQ=255)
                                                                                                              ataatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGtttgcttt                                      >  1:3645249/1‑150 (MQ=255)
                                                                                                               taatGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTGTATATTTGCAGTCCGGtttgctttg                                     >  2:2756255/1‑150 (MQ=255)
                                                                                                                  tGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCAt                                  >  2:2257665/1‑150 (MQ=255)
                                                                                                                     ccGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATAcc                               <  1:2392720/150‑1 (MQ=255)
                                                                                                                     ccGGTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATAcc                               <  2:3663146/150‑1 (MQ=255)
                                                                                                                      cGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCg                              >  2:3578323/1‑150 (MQ=255)
                                                                                                                       ggTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCgg                             <  2:1251658/150‑1 (MQ=255)
                                                                                                                        gTACTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGa                            >  1:1071419/1‑150 (MQ=255)
                                                                                                                         tACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAt                           <  2:941010/150‑1 (MQ=255)
                                                                                                                           cTCATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttt                         >  1:161273/1‑150 (MQ=255)
                                                                                                                            tCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAtttt                        >  1:2773840/1‑150 (MQ=255)
                                                                                                                             cATGTTTTCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttc                       >  1:3534108/1‑150 (MQ=255)
                                                                                                                                 ttttCGGGTTTATGGTTTCTTATGAAATATATTGAATTCTCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctttt                   >  2:2053047/1‑150 (MQ=255)
                                                                                                                                  tttCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttcttttt                  <  2:3410135/150‑1 (MQ=255)
                                                                                                                                    tCGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctttttct                >  1:260192/1‑150 (MQ=255)
                                                                                                                                    tCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctttttct                <  2:4448150/150‑1 (MQ=255)
                                                                                                                                    gcGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGAttttctttttct                <  2:4448413/149‑1 (MQ=255)
                                                                                                                                     cGGGTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATtttctttttctt               <  1:132822/150‑1 (MQ=255)
                                                                                                                                      gggTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTa              >  1:2894293/1‑150 (MQ=255)
                                                                                                                                       ggTTTATGGTTTCTTATGTAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTAc             <  2:4489981/150‑1 (MQ=255)
                                                                                                                                        gTTTATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTAcc            >  1:2822319/1‑150 (MQ=255)
                                                                                                                                         tttATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCa           <  2:161273/150‑1 (MQ=255)
                                                                                                                                          ttATGGTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCAt          <  2:215431/150‑1 (MQ=255)
                                                                                                                                              ggTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTATATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTg      >  2:3223288/1‑150 (MQ=255)
                                                                                                                                              ggTTTCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTg      >  1:4380040/1‑150 (MQ=255)
                                                                                                                                                tttCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGaa    >  2:2929793/1‑150 (MQ=255)
                                                                                                                                                 ttCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAg   >  2:3885568/1‑150 (MQ=255)
                                                                                                                                                  tCTTATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGt  <  1:440769/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CGATAACCCTCTGGTGATCCGCGGACCTGGCGCTGGGCGCGACGTCACCGCCGGGGCGATTCAGTCGGATATCAACCGGCTGGCACAGTTGTTGTAATTTCAGAAATTTAATAATGCCCGGTACTCATGTTTTCGGGTTTATGGTTTCTAATGAAATATATTGAATTATCATAGGATTAGGCCGGATTAAGCGTTTACGACGAATCCGGCAAGAAGCAATAAGTACATGGTTAGTTTATATTTGCAGTCCGGTTTGCTTTGCATACCGGATTTTCTTTTTCTTACCATCCTGAAGT  >  CP009273/4122099‑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: