Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,054,888 A→T 8.2% V186E (GTA→GAA)  ygfB ← UPF0149 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,054,8880AT 8.2% 1694.9 / 32.7 549V186E (GTA→GAA) ygfBUPF0149 family protein
Reads supporting (aligned to +/- strand):  ref base A (186/310);  new base T (46/3);  total (236/313)
Fisher's exact test for biased strand distribution p-value = 2.41e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 557 total aligned reads displayed.

TACTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGA  >  NC_000913/3054739‑3055037
                                                                                                                                                     |                                                                                                                                                     
tACTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTa                                                                                                                                                       <  2:5671901/150‑1 (MQ=255)
   ttCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACtt                                                                                                                                                    >  2:229417/1‑150 (MQ=255)
      tGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTg                                                                                                                                                 <  1:7086888/150‑1 (MQ=255)
       ggTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTgg                                                                                                                                                <  1:2139923/150‑1 (MQ=255)
       ggTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTgg                                                                                                                                                <  2:3421403/150‑1 (MQ=255)
       ggTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTgg                                                                                                                                                <  2:3641042/150‑1 (MQ=255)
        gTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGTAATTCTGgc                                                                                                                                               >  1:7358432/1‑150 (MQ=255)
         tGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGgcg                                                                                                                                              >  1:4995365/1‑150 (MQ=255)
           cagcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGcgcg                                                                                                                                            <  2:2859083/150‑1 (MQ=255)
           cagcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGcgcg                                                                                                                                            <  1:3019174/150‑1 (MQ=255)
            agcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTGTTGTACTTCTGGGGcgg                                                                                                                                           >  2:8082302/1‑150 (MQ=255)
            agcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTGACGTTTTTTGTTGTTAGTGGAGAGTCGGGTGTTGTACGTCTGGCGcgg                                                                                                                                           >  2:4666019/1‑150 (MQ=255)
             gcaAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTGATTGTGGCGcggg                                                                                                                                          >  1:9979140/1‑149 (MQ=255)
              caaccctCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATCTCTCACTCATAACACTCTCCTTCCGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGcggtc                                                                                                                                         <  2:9531259/144‑1 (MQ=255)
              caAAAATCCGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGTCGcgggc                                                                                                                                         >  2:1044866/1‑148 (MQ=255)
              caAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTGGACTTATGGCGcgggc                                                                                                                                         >  2:1284401/1‑148 (MQ=255)
              caAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTGTTAGTGTCGAGTCGGTTTTTGTTGATATGGCGcggtc                                                                                                                                         >  2:9718023/1‑150 (MQ=255)
                   aTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTg                                                                                                                                    <  1:7924805/150‑1 (MQ=255)
                   aTCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTg                                                                                                                                    <  1:10020896/150‑1 (MQ=255)
                    tCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGc                                                                                                                                   <  1:4775663/150‑1 (MQ=255)
                     cAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGATGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCg                                                                                                                                  <  1:8688149/150‑1 (MQ=255)
                       gcgcGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTGACGTGGGTTGTTGGGCGGGGAGAGTCGGGTGGTGGGCGGCTGGCGCGGGCGGGTGCggc                                                                                                                                >  2:7235704/1‑149 (MQ=255)
                         gcggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTAATTATTGCGCGGTCGGTTGCGGATg                                                                                                                              >  1:1289096/1‑150 (MQ=255)
                         gcggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATg                                                                                                                              <  2:6675537/150‑1 (MQ=255)
                          cggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGa                                                                                                                             <  1:5173880/150‑1 (MQ=255)
                           ggcgCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTCTTATTGGGCGGTCGGGTGCGGGTGAg                                                                                                                            >  1:9698543/1‑150 (MQ=255)
                              tcTGCCGGGGTGGATTTGCTCCACCAGGTCCTGACGGCGAAGCTGAATCTCTTGGCGTGATATCTCAATCATATCTGTGTCATTTCGTTTTTTGTTTTTAGTGTCGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGTTGAGTaa                                                                                                                         <  2:757590/149‑1 (MQ=255)
                               cTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTaaa                                                                                                                        <  1:6797817/150‑1 (MQ=255)
                                tGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAg                                                                                                                       <  2:9437627/150‑1 (MQ=255)
                                tGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTAATTGTGGCGCGGTCGGTTGCGGGTGTGTAAAg                                                                                                                       >  1:8572039/1‑150 (MQ=255)
                                  ccGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGgt                                                                                                                     <  2:65769/150‑1 (MQ=255)
                                  ccGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTATGGCGCGGGCGGTTTCGGGTGAGTAAAGgt                                                                                                                     >  1:2424959/1‑150 (MQ=255)
                                          cATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGc                                                                                                             <  2:2644700/150‑1 (MQ=255)
                                          cATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCGTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTCGTCCGGCGCGGGCGGGTGCGGCTGGGTAAAGGGGGCGTGGc                                                                                                             >  2:8722650/1‑150 (MQ=255)
                                                 tCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTGTGGACGTCTGGCGCGGTCGGGTGCGGCGGAGTCAAGGGGTCGGGGCAGAACAg                                                                                                      >  2:6424643/1‑150 (MQ=255)
                                                    accaGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAgcgc                                                                                                   <  2:5943220/150‑1 (MQ=255)
                                                    accaGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAgcgc                                                                                                   <  2:2543159/150‑1 (MQ=255)
                                                     ccaGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGcgcg                                                                                                  <  1:6711362/150‑1 (MQ=255)
                                                     ccaGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTGTTTGTAGTTCTGGCGCGGGGGGGTGCGGCTGGGTAAAGGTGTCGTGGCATAACAGcgcg                                                                                                  >  1:5339655/1‑150 (MQ=255)
                                                       aGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGc                                                                                                <  1:10180340/150‑1 (MQ=255)
                                                         ggCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTACTTCTGGCGCGGGCGGGTGCGGGTGTGTAAAGGTGTCGTGGGCTAACAGCGCGGcca                                                                                              >  1:1128255/1‑148 (MQ=255)
                                                         ggCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGGGTAGAGTCGGTTTTTGTAATTATGGCGCGGTCCGTTGCGGCGGAGTGAAGGTGTCGTGGCCTAAAAGCGCGGCaa                                                                                              >  2:9143652/1‑150 (MQ=255)
                                                          gCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTTTAGTTATGGCGCGGTGGGTGGCGGGTGTGTTAAGGTGTCGTGGGATAACAGCGCGGCacc                                                                                             >  1:2226682/1‑148 (MQ=255)
                                                           ccTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAaca                                                                                            <  2:9100544/150‑1 (MQ=255)
                                                           ccTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAaca                                                                                            <  2:8612835/150‑1 (MQ=255)
                                                             tGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTAATTATGGCGCGGTCGGTTGCGGGTGAGTAAAGGTTTCGTGGCATAACAGCGCGGCAACCCg                                                                                          >  2:402860/1‑150 (MQ=255)
                                                               aCGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGaa                                                                                        <  1:9151901/150‑1 (MQ=255)
                                                               aCGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGaa                                                                                        <  1:2305116/150‑1 (MQ=255)
                                                                cGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAAc                                                                                       <  2:7326439/150‑1 (MQ=255)
                                                                 ggCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGGTTTGTTCATCAGGCGCGGGGGGGTGCGGGGGCGGGAAGGTGGCGGGGCAGAACAGCCCGGCCCCACGCAcc                                                                                      >  1:6326411/1‑149 (MQ=255)
                                                                     aCGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTAtt                                                                                  <  2:2588614/150‑1 (MQ=255)
                                                                        cTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTATTATTGGGCGGTCGGTTGCGGGGTTGGTAAGGTTTGGGGGCAAAAAAGCGCGGCAACACCAACGTaataga                                                                               >  1:6161844/1‑145 (MQ=255)
                                                                        cTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCGTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCga                                                                               <  2:3100572/150‑1 (MQ=255)
                                                                         tGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCgat                                                                              <  1:4279103/150‑1 (MQ=255)
                                                                          gAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTAATTCAGGCGCGGTCGGTTGCGGGTGTGTAAAGGTGTGGTGGCATAACAGCGCGGCAACACCAACGGATTAgatg                                                                             >  1:5985789/1‑150 (MQ=255)
                                                                          gAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCgatg                                                                             <  1:2651541/150‑1 (MQ=255)
                                                                           aaaCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTTCTTGTTGGGCGGTGGGTTTCGGATGAGGAAAGGTGGCGTGGCCTAACAGCGCGGCCACACCAACCTATTAgatga                                                                            >  1:9013168/1‑150 (MQ=255)
                                                                           aaaCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGGTTTTTTACTTCTGGCGCGGGGGGGTTGGGCTGAGTAAAAGTGGGGTGGGATAAAAGCGCGGCCACACCAACCGATTCggagt                                                                            >  1:3890434/1‑146 (MQ=255)
                                                                            aaCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCgatgat                                                                           <  2:6213372/150‑1 (MQ=255)
                                                                            aaCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCgatgat                                                                           <  1:8768758/150‑1 (MQ=255)
                                                                            aaCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCgatgat                                                                           <  2:7783722/150‑1 (MQ=255)
                                                                            aaCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCgatgat                                                                           <  1:8297518/150‑1 (MQ=255)
                                                                            aaCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCgatgat                                                                           <  1:6538193/150‑1 (MQ=255)
                                                                             aCTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGAtc                                                                          <  1:6597885/150‑1 (MQ=255)
                                                                               tctTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCAGGCGCGGTCGGTTTCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGAtcac                                                                        >  1:1600119/1‑148 (MQ=255)
                                                                                ctTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATctct                                                                       <  1:5686118/150‑1 (MQ=255)
                                                                                ctTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATctct                                                                       <  2:9327808/150‑1 (MQ=255)
                                                                                  tGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTc                                                                     <  2:541447/150‑1 (MQ=255)
                                                                                     cGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACGTCTGGCGCGGTCGGTTGCGGGTGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAg                                                                  >  2:1624114/1‑150 (MQ=255)
                                                                                      gggATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTAAATATTGCGCGGGCGGGTGCGGGTGAGTAAAGGTGTGGTGGCATAAAAGCGCGGCCACACGAACCTTTTAGATTATCTCTTAAAGc                                                                 >  1:5482122/1‑150 (MQ=255)
                                                                                      gggATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGTATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGc                                                                 <  2:4543235/150‑1 (MQ=255)
                                                                                      gggATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGc                                                                 <  1:4666019/150‑1 (MQ=255)
                                                                                           atCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACAt                                                            <  2:7033243/150‑1 (MQ=255)
                                                                                              tcactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTACTTATGGCGCCGGCGGTTGCGGATGAGGAAAGGTGTCGTGGCATAACAGCGCGGCAACACCAACGTATTAGATGATATATTCAAGAGACAttta                                                         >  2:2201899/1‑149 (MQ=255)
                                                                                               cactcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACGTCAGGCGCGGTCGGGTGCGGATGAGTAAAGGTGTCGGGGCAGAACAGCCCGGCAACACGAACCTAATCGCTGATCACTTCCACCGCCAATTCa                                                        >  2:6785466/1‑150 (MQ=255)
                                                                                                actcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTTTTATTATTGGGCGGGCGGGTGCGGGGGGGTTAAGGTGTGGGGGCATAAAAGCGCGGCCACACCAAACTATTCGGTTATCTATTAAAGCGACAGTTaaa                                                       >  2:52266/1‑150 (MQ=255)
                                                                                                actcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCaa                                                       <  1:3279411/150‑1 (MQ=255)
                                                                                                actcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCaa                                                       <  1:7477481/150‑1 (MQ=255)
                                                                                                actcaTAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACATCTGGCGCGGGCGGTTGCGGCTGGGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCaa                                                       >  2:1513162/1‑150 (MQ=255)
                                                                                                           tctcCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGa                                            <  1:5572482/150‑1 (MQ=255)
                                                                                                                ttACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCTGTTGCGGATGAGTAAAGGTGTCGTGTCATAACAGCGCGGCAACAAGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGAtctta                                       <  2:6851237/150‑2 (MQ=255)
                                                                                                                ttACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGAtcttc                                       <  2:10129444/150‑1 (MQ=255)
                                                                                                                gtACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGAtcttc                                       <  1:9387745/149‑1 (MQ=255)
                                                                                                                  aCGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTAATTATGGCGCGGTCGGTTGCGGGTGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGAtcatagt                                     >  2:5855869/1‑145 (MQ=255)
                                                                                                                    gTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGAtcttcgtct                                   <  2:5021733/150‑1 (MQ=255)
                                                                                                                         ttGTTTTTCGTGTAGAGTCGGTTTTTGTACTTCGGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCACCACGAACGTATTCGATGAGCTCTTCAAGCGACATTTCACGCTCTTCCTGATCttcgtcttcgtc                              <  2:4832378/150‑1 (MQ=255)
                                                                                                                             ttttAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTaa                          <  2:2004506/150‑1 (MQ=255)
                                                                                                                               ttAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAAcc                        <  2:4186070/150‑1 (MQ=255)
                                                                                                                                    gtAGAGTCGGTTTTTTTACTTCATGCGCGGTCGGGTTCGGATGAGTAAAGGTGTCGTGGCATAAAACCGCGGCAACACCAACCTATTAGAAGAAATATTCAAGCGACAATTCAAGCACTTACTGATAATAGTCTTAGTCGTAACCCAGtt                   >  2:8405372/1‑150 (MQ=255)
                                                                                                                                    gtAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGtt                   <  1:589738/150‑1 (MQ=255)
                                                                                                                                      agagTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTgc                 <  2:9199150/150‑1 (MQ=255)
                                                                                                                                         gTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCa              <  2:7966381/150‑1 (MQ=255)
                                                                                                                                          tCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCaa             >  2:3250761/1‑150 (MQ=255)
                                                                                                                                           cGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAAt            <  2:5740635/150‑1 (MQ=255)
                                                                                                                                           cGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAAt            <  1:4088761/150‑1 (MQ=255)
                                                                                                                                            ggTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATg           <  1:3409366/150‑1 (MQ=255)
                                                                                                                                             gTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGt          >  1:3758107/1‑150 (MQ=255)
                                                                                                                                             gTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGt          >  1:1811216/1‑150 (MQ=255)
                                                                                                                                             gTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGt          >  2:9968186/1‑150 (MQ=255)
                                                                                                                                                 ttGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACg      <  2:4453091/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACTTCTGGTGCAGCAAAAATCAGCGCGGCGCTGCCGGGTTGCATTTGCTCCACCAGGGCCTGACGGCGACGCTGAAACTCTTGCCGGGATATCTCACTCATAACACTCTCCTTACGTTTTTTGTTTTTAGTGTAGAGTCGGTTTTTGTACTTCTGGCGCGGTCGGTTGCGGATGAGTAAAGGTGTCGTGGCATAACAGCGCGGCAACACGAACGTATTCGATGATCTCTTCAAGCGACATTTCAAGCTCTTCCTGATCTTCGTCTTCGTCGTAACCCAGTTGCGCAATGTTACGCAGA  >  NC_000913/3054739‑3055037

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: