Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,956,973:1 +A 6.8% coding (848/894 nt) ilvY ← transcriptional activator of ilvC; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,956,9731.A 6.8% 384.8 / 24.8 103coding (848/894 nt)ilvYtranscriptional activator of ilvC; autorepressor
Reads supporting (aligned to +/- strand):  ref base . (48/48);  new base A (3/4);  total (51/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

Only 100 of 103 total aligned reads displayed.

CCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGCCGAGTG  >  NC_000913/3956827‑3957120
                                                                                                                                                   |                                                                                                                                                   
ccACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCAt                                                                                                                                                  >  1:1015623/1‑150 (MQ=255)
 cACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCAtt                                                                                                                                                 >  2:343182/1‑149 (MQ=255)
     acgaAACCAATGACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGc                                                                                                                                             <  2:1452785/150‑1 (MQ=255)
     acgaAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGc                                                                                                                                             <  1:1201939/150‑1 (MQ=255)
       gaAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCg                                                                                                                                           <  2:77291/150‑1 (MQ=255)
       gaAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCg                                                                                                                                           >  1:177675/1‑150 (MQ=255)
          aCCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCtt                                                                                                                                        <  1:1502194/150‑1 (MQ=255)
           ccAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCGTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCtaa                                                                                                                                       <  1:1515202/150‑3 (MQ=255)
            cAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCtaat                                                                                                                                      <  1:1640728/150‑4 (MQ=255)
             aaTAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCtaatt                                                                                                                                     <  1:1337454/150‑5 (MQ=255)
              aTAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTc                                                                                                                                    <  2:1503339/150‑1 (MQ=255)
                aaCCCGGCGTTCAGGGTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCtaattgag                                                                                                                                  <  1:68728/150‑8 (MQ=255)
                  cccGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTgcg                                                                                                                                >  2:581662/1‑150 (MQ=255)
                  cccGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTAGCTTATCAGGCCTACCCGCTCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTgcg                                                                                                                                <  1:1287348/150‑1 (MQ=255)
                    cGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGca                                                                                                                              <  1:362999/150‑1 (MQ=255)
                       cGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGcacac                                                                                                                           >  2:138118/1‑150 (MQ=255)
                        gTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGcacaca                                                                                                                          >  1:793255/1‑150 (MQ=255)
                         ttCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTTCCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCacacac                                                                                                                         <  1:866826/150‑1 (MQ=255)
                           cAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGc                                                                                                                       <  2:246161/150‑1 (MQ=255)
                              gTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCa                                                                                                                    <  1:334420/150‑1 (MQ=255)
                              gTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCa                                                                                                                    <  2:596489/150‑1 (MQ=255)
                               ttCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATAAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAg                                                                                                                   >  1:1076076/1‑150 (MQ=255)
                                tCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGc                                                                                                                  >  2:1053130/1‑150 (MQ=255)
                                gcTTTTTGGCGGGTTAGGGAAAAATGCCTGAGAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCTCACACGCCCAGc                                                                                                                  <  1:282234/149‑1 (MQ=255)
                                 cTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCt                                                                                                                 >  2:1264598/1‑150 (MQ=255)
                                 cTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCt                                                                                                                 >  1:1357168/1‑150 (MQ=255)
                                  tttttGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTc                                                                                                                >  1:771044/1‑150 (MQ=255)
                                      tGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGATCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCAGCTTTTTCTGCGCACACACGCCCAGCTCAAAc                                                                                                            >  2:1531438/1‑150 (MQ=255)
                                       ggCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACg                                                                                                           <  2:863690/150‑1 (MQ=255)
                                       ggCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACg                                                                                                           <  2:446923/150‑1 (MQ=255)
                                            gTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAACATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCACACGGCttt                                                                                                      >  1:331243/1‑150 (MQ=255)
                                            gTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACCCGCCCAGCTCAAACGGCGtt                                                                                                      >  1:336264/1‑150 (MQ=255)
                                             ttAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGttt                                                                                                     <  1:789414/150‑1 (MQ=255)
                                              tAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGtttt                                                                                                    >  1:956399/1‑150 (MQ=255)
                                               aGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATATAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGtttta                                                                                                   >  2:585386/1‑149 (MQ=255)
                                                gggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtct                                                                                                  >  2:499660/1‑150 (MQ=255)
                                                  gAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCa                                                                                                >  2:371809/1‑150 (MQ=255)
                                                  gAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCa                                                                                                >  1:1344400/1‑150 (MQ=255)
                                                        tGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTg                                                                                          >  1:692703/1‑150 (MQ=255)
                                                         gCCTGATAGCTCTTCGCTTATCAGGCATACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTgc                                                                                         <  2:1059305/150‑1 (MQ=255)
                                                         gCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTgc                                                                                         <  2:363904/150‑1 (MQ=255)
                                                         gCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTgc                                                                                         <  2:1064881/150‑1 (MQ=255)
                                                         gCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGGTTTCTCATCGCTgc                                                                                         >  2:291486/1‑150 (MQ=255)
                                                           cTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTgcgc                                                                                       >  1:220022/1‑150 (MQ=255)
                                                            tGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGct                                                                                      >  1:962942/1‑150 (MQ=255)
                                                             gATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGctc                                                                                     <  2:1102411/150‑1 (MQ=255)
                                                              aTAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGctct                                                                                    <  2:1171402/150‑1 (MQ=255)
                                                              aTAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGctct                                                                                    <  2:1134652/150‑1 (MQ=255)
                                                                   gcttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATATTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAat                                                                               >  1:1476557/1‑150 (MQ=255)
                                                                      tcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAatcat                                                                            >  1:1280038/1‑150 (MQ=255)
                                                                        gcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAAtcatca                                                                          <  1:578042/150‑1 (MQ=255)
                                                                        gcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCATCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAAtcatca                                                                          <  1:571216/150‑1 (MQ=255)
                                                                        gcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGGTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAAtcatca                                                                          <  2:944029/150‑1 (MQ=255)
                                                                         cttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCAc                                                                         >  1:1400851/1‑150 (MQ=255)
                                                                           tATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCAcgc                                                                       >  2:1582205/1‑150 (MQ=255)
                                                                              cAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGt                                                                    <  1:499660/150‑1 (MQ=255)
                                                                               aGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGtt                                                                   >  1:1192888/1‑150 (MQ=255)
                                                                                 gCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTAc                                                                 <  1:1520943/150‑1 (MQ=255)
                                                                                  ccTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACg                                                                >  1:1186816/1‑150 (MQ=255)
                                                                                  ccTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACg                                                                >  1:1711669/1‑150 (MQ=255)
                                                                                      cccGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCAcc                                                            <  1:1327558/150‑1 (MQ=255)
                                                                                       ccGCGCGTCAACGTCATTTGTGGGTCGGCAAAATCTTCCAGAATGCCTCTATTATCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCTAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCg                                                           <  1:1146988/150‑1 (MQ=255)
                                                                                       ccGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCg                                                           >  2:1460144/1‑150 (MQ=255)
                                                                                        cgcgcgACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCgg                                                          >  2:409003/1‑150 (MQ=255)
                                                                                        cgcgcgACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCgg                                                          >  1:895361/1‑150 (MQ=255)
                                                                                        cgcgcgACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCgg                                                          >  2:603879/1‑150 (MQ=255)
                                                                                             gACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCgg                                                     <  1:371809/150‑1 (MQ=255)
                                                                                               cAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCgggg                                                   <  2:1344400/150‑1 (MQ=255)
                                                                                               cAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCgggg                                                   >  1:611067/1‑150 (MQ=255)
                                                                                                aaCGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGc                                                  >  2:1104970/1‑150 (MQ=255)
                                                                                                aaCGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGc                                                  <  2:220022/150‑1 (MQ=255)
                                                                                                  cGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTg                                                <  1:397152/150‑1 (MQ=255)
                                                                                                    tCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCACTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGtt                                              >  2:361721/1‑150 (MQ=255)
                                                                                                     cATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGttt                                             <  2:1076076/150‑1 (MQ=255)
                                                                                                      aTTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTATGCACCGGTTCGGGGCTGtttt                                            >  2:817202/1‑150 (MQ=255)
                                                                                                           tgGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGc                                       <  2:1542155/150‑1 (MQ=255)
                                                                                                               tCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCAGcaccac                                  <  2:1015623/150‑1 (MQ=255)
                                                                                                                 ggCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACt                                 >  1:1636028/1‑150 (MQ=255)
                                                                                                                  gCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACtt                                >  2:1273161/1‑150 (MQ=255)
                                                                                                                    aaaaTCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTcc                              <  2:692703/150‑1 (MQ=255)
                                                                                                                     aaaTCTTCCGGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCg                             <  2:513462/150‑1 (MQ=255)
                                                                                                                       aTCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCTTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGc                           <  1:1689834/150‑1 (MQ=255)
                                                                                                                        tCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCACTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGca                          <  1:361721/150‑1 (MQ=255)
                                                                                                                         cTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGcaa                         >  2:1063917/1‑150 (MQ=255)
                                                                                                                         cTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCACAAACGGAGTTATCTCATCGCTGCGCTCTAAAAACATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGcaa                         >  1:658118/1‑150 (MQ=255)
                                                                                                                          ttCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGcaac                        >  1:597858/1‑150 (MQ=255)
                                                                                                                           tCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGcaaca                       <  2:956399/150‑1 (MQ=255)
                                                                                                                               gAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATgc                   <  2:1192888/150‑1 (MQ=255)
                                                                                                                                 aTGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATgcca                 >  1:647971/1‑150 (MQ=255)
                                                                                                                                    ccTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCCCGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATgccccgc              >  1:1077795/1‑150 (MQ=255)
                                                                                                                                    ccTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATgcc                  <  1:1604086/146‑1 (MQ=255)
                                                                                                                                    ccTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATgcc                  >  2:1604086/1‑146 (MQ=255)
                                                                                                                                       cagttAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATgccacgccac           >  2:426025/4‑150 (MQ=255)
                                                                                                                                        aaaaaGCGGCTACATGAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAACCGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGTTTCGGGGCTGTTATCCAGCGCCAATTCCGGCAACAATGCCACGCcaca          >  2:1052872/5‑150 (MQ=255)
                                                                                                                                          ttAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGc        <  1:988028/150‑1 (MQ=255)
                                                                                                                                          aaaGCGGCTACATGAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAACCGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGc        >  1:110470/3‑150 (MQ=255)
                                                                                                                                           tAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGcc       <  1:1249195/150‑1 (MQ=255)
                                                                                                                                            aGCGGCTACATGAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGCCg      >  1:978011/1‑150 (MQ=255)
                                                                                                                                             gCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGCCGa     <  2:793255/150‑1 (MQ=255)
                                                                                                                                                gCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGCCGAGTg  <  2:481416/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGCCGAGTG  >  NC_000913/3956827‑3957120

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: