Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,233,642 G→A 22.3% 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,6420GA22.3% 540.2 / 146.4 251intergenic (‑409/‑116)lysC/pgilysine‑sensitive aspartokinase 3/glucosephosphate isomerase
Reads supporting (aligned to +/- strand):  ref base G (100/95);  new base A (19/37);  total (119/132)
Fisher's exact test for biased strand distribution p-value = 2.35e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 257 total aligned reads displayed.

ACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCC  >  NC_000913/4233494‑4233790
                                                                                                                                                    |                                                                                                                                                    
aCTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTgc                                                                                                                                                     <  1:1166165/150‑1 (MQ=255)
 cTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTgcg                                                                                                                                                    >  1:2707867/1‑150 (MQ=255)
 cTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTacg                                                                                                                                                    >  1:2007955/1‑150 (MQ=255)
  ttGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTgcgc                                                                                                                                                   <  1:4008213/150‑1 (MQ=255)
     ggCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTaa                                                                                                                                                >  1:3658461/1‑150 (MQ=255)
       cAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAgg                                                                                                                                              >  2:3981868/1‑150 (MQ=255)
       cAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAgg                                                                                                                                              >  2:690334/1‑150 (MQ=255)
       cAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAgg                                                                                                                                              <  1:2852316/150‑1 (MQ=255)
       cAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAgg                                                                                                                                              <  1:193212/150‑1 (MQ=255)
         aCGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGt                                                                                                                                            <  1:1291841/150‑1 (MQ=255)
           gcgcCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGttt                                                                                                                                          >  2:2711945/1‑150 (MQ=255)
            cgcCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTa                                                                                                                                         <  2:3589902/150‑1 (MQ=255)
              ccTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTaca                                                                                                                                       <  2:2293424/150‑1 (MQ=255)
                tcttctGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACt                                                                                                                                     >  2:1553067/1‑150 (MQ=255)
                tcttctGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACt                                                                                                                                     >  2:3534267/1‑150 (MQ=255)
                tcttctGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACt                                                                                                                                     >  2:2187794/1‑150 (MQ=255)
                tcttctGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGATTCACTGCGCTAAGGGTTTACACt                                                                                                                                     >  2:2527978/1‑150 (MQ=255)
                   tctGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCaa                                                                                                                                  <  2:1166119/150‑1 (MQ=255)
                   tctGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCaa                                                                                                                                  <  2:2442914/150‑1 (MQ=255)
                   tctGGCCTGCGCGAGGGCAGGTAGTAGATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCTACGTTGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCaa                                                                                                                                  <  2:3809336/150‑1 (MQ=255)
                   tctGGCCTGCGCGAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCaa                                                                                                                                  <  2:2785080/150‑1 (MQ=255)
                      taccTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAg                                                                                                                                               >  2:1048990/3‑134 (MQ=255)
                      taccTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAg                                                                                                                                               <  1:1048990/132‑1 (MQ=255)
                       gCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACAtt                                                                                                                              >  1:595640/1‑150 (MQ=255)
                         cTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTAc                                                                                                                            >  1:2046806/1‑150 (MQ=255)
                          tGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACg                                                                                                                           >  1:4048014/1‑150 (MQ=255)
                          tGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACg                                                                                                                           >  1:3808943/1‑150 (MQ=255)
                          tGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGATCAGCTAAAAAATTCTTCACAGCGCTAAGGTTTTACACTCAACAATaac                                                                                                                           >  1:312764/1‑148 (MQ=255)
                          agcgcTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACg                                                                                                                           <  2:1344375/149‑1 (MQ=255)
                            cgcTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCt                                                                                                                         <  2:1174253/150‑1 (MQ=255)
                              cTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTaa                                                                                                                       <  1:3412236/150‑1 (MQ=255)
                               tAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAAc                                                                                                                      <  1:1003732/150‑1 (MQ=255)
                                   gcAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTAAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCa                                                                                                                  <  1:26749/150‑1 (MQ=255)
                                         agtaCATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTaaaa                                                                                                            <  1:2524043/150‑1 (MQ=255)
                                           taCATTTATAAATAAAGGGTGAGCGGGGTGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAAcc                                                                                                          <  1:3085999/150‑1 (MQ=255)
                                           taCATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAAcc                                                                                                          >  1:378751/1‑150 (MQ=255)
                                             cATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCAt                                                                                                        <  1:101889/150‑1 (MQ=255)
                                              aTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATc                                                                                                       >  1:1818933/1‑150 (MQ=255)
                                              aTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATc                                                                                                       >  1:1575478/1‑150 (MQ=255)
                                               tttATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATca                                                                                                      >  1:442106/1‑150 (MQ=255)
                                               tttATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATca                                                                                                      <  2:3786767/150‑1 (MQ=255)
                                                 tataAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATcaca                                                                                                    <  2:1963556/150‑1 (MQ=255)
                                                 tataAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATcaca                                                                                                    <  2:4117912/150‑1 (MQ=255)
                                                  ataaataaAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACAt                                                                                                   >  1:1760231/1‑150 (MQ=255)
                                                      ataaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACAttttt                                                                                               <  2:188234/150‑1 (MQ=255)
                                                       taaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTc                                                                                              >  1:449192/1‑150 (MQ=255)
                                                       taaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTc                                                                                              >  2:592057/1‑150 (MQ=255)
                                                         aaGGGTGAGCGGGGCGGTTGTCAACGTTGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCtg                                                                                            <  1:450509/150‑1 (MQ=255)
                                                             gTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGAc                                                                                        >  1:3469733/1‑150 (MQ=255)
                                                             gTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGAc                                                                                        >  1:3358137/1‑150 (MQ=255)
                                                                aGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTgg                                                                                     <  1:3293929/150‑1 (MQ=255)
                                                                  cGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGgcg                                                                                   <  1:1162575/150‑1 (MQ=255)
                                                                        ggTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACaa                                                                             <  1:1646592/150‑1 (MQ=255)
                                                                          ttGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATc                                                                           >  1:848718/1‑150 (MQ=255)
                                                                          ttGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATc                                                                           <  1:495440/150‑1 (MQ=255)
                                                                            gTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACAATCACATTATTCTGTGACTGGCGCTACAATCtt                                                                         >  2:3677350/1‑150 (MQ=255)
                                                                             tCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTc                                                                        <  1:2346337/150‑1 (MQ=255)
                                                                               aaCGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCa                                                                      <  1:1594290/150‑1 (MQ=255)
                                                                                  gATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAg                                                                   >  2:3100091/1‑150 (MQ=255)
                                                                                    tGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTc                                                                 >  1:1231887/1‑150 (MQ=255)
                                                                                        gTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACaa                                                             >  2:2964858/1‑150 (MQ=255)
                                                                                             gCGGTTTTTTCATACACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCATCATTACGATAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCAAAAGTCACAATtctc                                                        <  2:2219370/150‑1 (MQ=255)
                                                                                             gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtctc                                                        <  2:3186021/150‑1 (MQ=255)
                                                                                             gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtctc                                                        >  1:190199/1‑150 (MQ=255)
                                                                                             gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtctc                                                        >  1:2859001/1‑150 (MQ=255)
                                                                                             gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtctc                                                        <  1:3571089/150‑1 (MQ=255)
                                                                                             gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtctc                                                        <  2:482540/150‑1 (MQ=255)
                                                                                             gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtctc                                                        <  1:927335/150‑1 (MQ=255)
                                                                                             gCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATtctc                                                        <  2:344263/150‑1 (MQ=255)
                                                                                              cGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCATAGTCACAATTCTCa                                                       <  1:735981/150‑1 (MQ=255)
                                                                                               ggATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACTCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCaa                                                      >  1:1410026/1‑150 (MQ=255)
                                                                                               ggATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCaa                                                      >  2:2107355/1‑150 (MQ=255)
                                                                                               ggATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCaa                                                      >  1:2903447/1‑150 (MQ=255)
                                                                                               ggATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCaa                                                      >  2:2107406/1‑150 (MQ=255)
                                                                                               ggATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCaa                                                      >  2:260859/1‑150 (MQ=255)
                                                                                                   ttttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATc                                                  <  2:1291926/150‑1 (MQ=255)
                                                                                                    tttCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCa                                                 <  1:2527978/150‑1 (MQ=255)
                                                                                                         tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                            >  2:188861/1‑150 (MQ=255)
                                                                                                         tCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCagaaga                                            <  1:3662740/150‑1 (MQ=255)
                                                                                                            aCTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTa                                         <  2:2064590/150‑1 (MQ=255)
                                                                                                            aCTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTa                                         >  2:3840993/1‑150 (MQ=255)
                                                                                                              tCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTAtt                                       >  1:1517330/1‑150 (MQ=255)
                                                                                                                     ggTCAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATg                                >  2:4048436/1‑150 (MQ=255)
                                                                                                                        cAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGaaa                             <  2:1998163/150‑1 (MQ=255)
                                                                                                                        cAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGaaa                             >  1:3277312/1‑150 (MQ=255)
                                                                                                                        cAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGaaa                             <  2:2822472/150‑1 (MQ=255)
                                                                                                                        cAGTAGTTCAGCTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGaaa                             <  2:1002314/150‑1 (MQ=255)
                                                                                                                              ttCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATc                       <  2:83009/150‑1 (MQ=255)
                                                                                                                               tCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCa                      <  1:212246/150‑1 (MQ=255)
                                                                                                                                 aGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAAt                    >  2:484768/1‑150 (MQ=255)
                                                                                                                                  gCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATc                   <  2:2801702/150‑1 (MQ=255)
                                                                                                                                   cTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATcc                  >  1:1423181/1‑150 (MQ=255)
                                                                                                                                   cTAATAAATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATcc                  >  1:360139/1‑150 (MQ=255)
                                                                                                                                    taataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCa                 >  2:2658671/1‑150 (MQ=255)
                                                                                                                                     aataaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCaa                >  1:1316382/1‑150 (MQ=255)
                                                                                                                                      ataaATGCTTCACTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATACAAc               >  2:3623607/1‑150 (MQ=255)
                                                                                                                                       taaATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACg              >  1:2185202/1‑150 (MQ=255)
                                                                                                                                                 aCTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTg    <  2:3399042/150‑1 (MQ=255)
                                                                                                                                                 aCTACGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTg    >  2:3479604/1‑150 (MQ=255)
                                                                                                                                                   tGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGcc  >  2:3283785/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGGCGCTACAATCTTCCAAAGTCACAATTCTCAAAATCAGAAGAGTATTGCTAATGAAAAACATCAATCCAACGCAGACCGCTGCC  >  NC_000913/4233494‑4233790

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: