Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,233,642 G→A 23.6% intergenic (‑409/‑116) lysC ← / → pgi lysine‑sensitive aspartokinase 3/glucosephosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,233,6420GA23.6% 789.7 / 249.4 385intergenic (‑409/‑116)lysC/pgilysine‑sensitive aspartokinase 3/glucosephosphate isomerase
Reads supporting (aligned to +/- strand):  ref base G (163/131);  new base A (42/49);  total (205/180)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 389 total aligned reads displayed.

ACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGC  >  NC_000913/4233494‑4233789
                                                                                                                                                    |                                                                                                                                                   
          cgcgCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGtt                                                                                                                                          <  2:5422246/150‑1 (MQ=255)
          cgcgCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGtt                                                                                                                                          <  1:5833003/150‑1 (MQ=255)
           gcgcCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGttt                                                                                                                                         >  1:1112809/1‑150 (MQ=255)
           gcgcCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGttt                                                                                                                                         >  1:3374360/1‑150 (MQ=255)
             gcCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTac                                                                                                                                       >  1:3568176/1‑150 (MQ=255)
              ccTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGAGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTaca                                                                                                                                      >  1:5030744/1‑150 (MQ=255)
               ctctTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTacac                                                                                                                                     <  1:5966961/150‑1 (MQ=255)
                tcttctGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACt                                                                                                                                    >  2:2929757/1‑150 (MQ=255)
                 cttctGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTc                                                                                                                                   <  1:389293/150‑1 (MQ=255)
                 cttctGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTc                                                                                                                                   <  1:3110026/150‑1 (MQ=255)
                   tctGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCaa                                                                                                                                 <  1:1526968/150‑1 (MQ=255)
                   tctGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCaa                                                                                                                                 <  2:6588706/150‑1 (MQ=255)
                   tctGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCaa                                                                                                                                 <  1:5136420/150‑1 (MQ=255)
                   tctGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCaa                                                                                                                                 <  1:635583/150‑1 (MQ=255)
                   tctGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCaa                                                                                                                                 <  2:5664271/150‑1 (MQ=255)
                    ctGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAAc                                                                                                                                <  1:2774420/150‑1 (MQ=255)
                       gCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACAtt                                                                                                                             >  2:4726752/1‑150 (MQ=255)
                        ccTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTa                                                                                                                            <  2:3025911/150‑1 (MQ=255)
                        ccTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTa                                                                                                                            <  2:1855506/150‑1 (MQ=255)
                        ccTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTa                                                                                                                            <  2:4024339/150‑1 (MQ=255)
                        ccTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTa                                                                                                                            <  1:4100303/150‑1 (MQ=255)
                          tGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACg                                                                                                                          <  1:6500267/150‑1 (MQ=255)
                          tGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACg                                                                                                                          >  1:3682989/1‑150 (MQ=255)
                          tGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACg                                                                                                                          >  2:176206/1‑150 (MQ=255)
                               tAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAAc                                                                                                                     >  1:2445088/1‑150 (MQ=255)
                               tAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAAc                                                                                                                     >  2:6256585/1‑150 (MQ=255)
                               tAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAAc                                                                                                                     <  1:4349319/150‑1 (MQ=255)
                                   gcAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCa                                                                                                                 <  2:214703/150‑1 (MQ=255)
                                    cAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCAc                                                                                                                <  2:3553684/150‑1 (MQ=255)
                                      ggTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTa                                                                                                              >  2:6447359/1‑150 (MQ=255)
                                        tagtaCATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTaaa                                                                                                            <  1:5913017/150‑1 (MQ=255)
                                              aTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATc                                                                                                      >  1:260289/1‑150 (MQ=255)
                                              aTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATc                                                                                                      >  1:2899264/1‑150 (MQ=255)
                                              aTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATc                                                                                                      >  2:678718/1‑150 (MQ=255)
                                               tttATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATca                                                                                                     <  2:3418999/150‑1 (MQ=255)
                                                ttATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACAATcac                                                                                                    >  2:1139392/1‑150 (MQ=255)
                                                 tataAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATcaca                                                                                                   >  2:5580573/1‑150 (MQ=255)
                                                        aaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCt                                                                                            <  2:4266958/150‑1 (MQ=255)
                                                        aaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCt                                                                                            >  2:1351290/1‑150 (MQ=255)
                                                        aaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCt                                                                                            >  2:4931055/1‑150 (MQ=255)
                                                          aGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCtgt                                                                                          <  1:4289004/150‑1 (MQ=255)
                                                              tGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACt                                                                                      >  2:3703311/1‑150 (MQ=255)
                                                                  cGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGgcg                                                                                  <  1:512330/150‑1 (MQ=255)
                                                                  cGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGgcg                                                                                  <  2:2821904/150‑1 (MQ=255)
                                                                  cGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGgcg                                                                                  <  1:5050764/150‑1 (MQ=255)
                                                                    gggCGGTGGTCCACGAGGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCt                                                                                <  2:5853406/150‑1 (MQ=255)
                                                                     ggCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTa                                                                               >  2:4977697/1‑150 (MQ=255)
                                                                     ggCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTa                                                                               >  2:511710/1‑150 (MQ=255)
                                                                      gCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTAc                                                                              >  1:1313990/1‑150 (MQ=255)
                                                                      gCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTAc                                                                              <  1:2951726/150‑1 (MQ=255)
                                                                        ggTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACaa                                                                            >  2:5505200/1‑150 (MQ=255)
                                                                         gTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAAt                                                                           <  1:3532665/150‑1 (MQ=255)
                                                                            gTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCtt                                                                        <  1:6558608/150‑1 (MQ=255)
                                                                              cAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTcc                                                                      >  1:6606668/1‑150 (MQ=255)
                                                                               aaCGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCa                                                                     >  1:5554394/1‑150 (MQ=255)
                                                                                     ggggTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTca                                                               <  2:1632970/150‑1 (MQ=255)
                                                                                       ggTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTcaca                                                             >  1:1009878/1‑150 (MQ=255)
                                                                                         tCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAAt                                                           <  2:4035580/150‑1 (MQ=255)
                                                                                            tGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtct                                                        <  2:1184357/150‑1 (MQ=255)
                                                                                             gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtctc                                                       >  1:1865520/1‑150 (MQ=255)
                                                                                             gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtctc                                                       <  2:5030744/150‑1 (MQ=255)
                                                                                             gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtctc                                                       >  1:1407038/1‑150 (MQ=255)
                                                                                              cGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCa                                                      >  2:4608300/1‑150 (MQ=255)
                                                                                               ggATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCaa                                                     >  2:79598/1‑150 (MQ=255)
                                                                                                gATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCaaa                                                    >  1:1810689/1‑150 (MQ=255)
                                                                                                gATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCaaa                                                    <  1:3238028/150‑1 (MQ=255)
                                                                                                    tttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCa                                                <  1:1984995/150‑1 (MQ=255)
                                                                                                    tttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCa                                                <  1:865591/150‑1 (MQ=255)
                                                                                                     ttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCag                                               >  2:5733156/1‑150 (MQ=255)
                                                                                                     ttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCag                                               <  1:167435/150‑1 (MQ=255)
                                                                                                        aTCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaag                                            >  2:236440/1‑150 (MQ=255)
                                                                                                            aCTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTa                                        >  2:6065456/1‑150 (MQ=255)
                                                                                                              tCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTAtt                                      <  2:5866549/150‑1 (MQ=255)
                                                                                                              tCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTAtt                                      >  1:4742515/1‑150 (MQ=255)
                                                                                                               ccTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTg                                     <  1:4824151/150‑1 (MQ=255)
                                                                                                               ccTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTg                                     <  1:4931055/150‑1 (MQ=255)
                                                                                                                  ggCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTa                                  <  2:1478373/150‑1 (MQ=255)
                                                                                                                        cAGTAGGTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGaaa                            <  1:1592945/150‑1 (MQ=255)
                                                                                                                         agtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCAAATGaaaa                           >  1:2282620/1‑150 (MQ=255)
                                                                                                                          gtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGaaaaa                          <  1:6110372/150‑1 (MQ=255)
                                                                                                                          gtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGaaaaa                          <  2:3514541/150‑1 (MQ=255)
                                                                                                                          gtagtTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGaaaaa                          <  1:4078462/150‑1 (MQ=255)
                                                                                                                               tCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCa                     <  2:3285520/150‑1 (MQ=255)
                                                                                                                                 aGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAAt                   >  1:3007126/1‑150 (MQ=255)
                                                                                                                                 aGCTAATAAATGCTTCACTACGCTAAGGGTTTACTCTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAAt                   <  2:2590415/150‑1 (MQ=255)
                                                                                                                                 aGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACTCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAAt                   <  1:783646/150‑1 (MQ=255)
                                                                                                                                 aGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAAt                   <  1:4494557/150‑1 (MQ=255)
                                                                                                                                  gCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATc                  >  1:698069/1‑150 (MQ=255)
                                                                                                                                   cTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATcc                 >  2:1065799/1‑150 (MQ=255)
                                                                                                                                    taataaATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCa                >  2:5189117/1‑150 (MQ=255)
                                                                                                                                      ataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAAc              >  1:582312/1‑150 (MQ=255)
                                                                                                                                       taaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACg             <  2:2067343/150‑1 (MQ=255)
                                                                                                                                         aaTGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCa           >  1:5754203/1‑150 (MQ=255)
                                                                                                                                           tGCTTCACTGCGCTAAGGGTTTACACTCAACATTCCGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGa         <  2:1280653/150‑1 (MQ=255)
                                                                                                                                            gCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGAc        <  2:2728503/150‑1 (MQ=255)
                                                                                                                                              ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCg      >  2:5444099/1‑150 (MQ=255)
                                                                                                                                              ttCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCg      >  1:3971788/1‑150 (MQ=255)
                                                                                                                                                 aCTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTg   <  1:5306054/150‑1 (MQ=255)
                                                                                                                                                 aCTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTg   <  1:2149281/150‑1 (MQ=255)
                                                                                                                                                  cTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGc  <  2:3770455/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGC  >  NC_000913/4233494‑4233789

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: