Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 303,121 Δ1 bp coding (131/615 nt) yagU → DUF1440 family inner membrane acid resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913303,1190C.98.8% 402.5 / ‑0.8 84coding (129/615 nt)yagUDUF1440 family inner membrane acid resistance protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (43/40);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

TACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATCCCAATGCGGCTGTTTA  >  NC_000913/302975‑303265
                                                                                                                                                 |                                                                                                                                                  
tACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTgg                                                                                                                                             >  2:740672/1‑151 (MQ=255)
 aCAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTgg                                                                                                                                             <  1:1349125/150‑1 (MQ=255)
 aCAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTgg                                                                                                                                             <  2:1007160/150‑1 (MQ=255)
  cAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGat                                                                                                                                           <  2:69977/151‑1 (MQ=255)
  cAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGat                                                                                                                                           >  2:1030673/1‑151 (MQ=255)
    aGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGatat                                                                                                                                         <  2:973820/151‑1 (MQ=255)
      cAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGt                                                                                                                                       >  1:723548/1‑151 (MQ=255)
       aGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGtt                                                                                                                                      <  2:992321/151‑1 (MQ=255)
       aGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGtt                                                                                                                                      <  1:1120351/151‑1 (MQ=255)
       aGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGtt                                                                                                                                      <  2:1146788/151‑1 (MQ=255)
          gTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTaa                                                                                                                                   <  2:1192796/151‑1 (MQ=255)
            aGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATg                                                                                                                                 >  1:627299/1‑151 (MQ=255)
             gctatGAATATATGTGAAAAAACTCCACCGAATCTCATACGTTATGGTCTTGCTGTATTCATTGTGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑CGAGGGGGGCTGAAGTTTCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATg                                                                                                                                 <  1:330798/148‑1 (MQ=255)
              atatGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCa                                                                                                                               <  2:408124/151‑1 (MQ=255)
                  gAATATATTTGAACAAACTCCACCGAACCGCATTCGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCTTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCgt                                                                                                                           <  2:260770/151‑1 (MQ=255)
                  gAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCgt                                                                                                                           <  1:204688/151‑1 (MQ=255)
                  gAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGACTGAAGTTCCATTGCCGCCACGTCG‑CCGTTGGATATGTTTAATGCAGCgt                                                                                                                           <  1:382392/151‑1 (MQ=255)
                    atatatTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCgtg                                                                                                                          >  2:599704/1‑150 (MQ=255)
                      atatTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTgg                                                                                                                       >  2:296601/1‑151 (MQ=255)
                      atatTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTgg                                                                                                                       >  1:61628/1‑151 (MQ=255)
                         tttGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGccc                                                                                                                    >  2:658093/1‑151 (MQ=255)
                          ttGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGccc                                                                                                                    <  1:1018323/150‑1 (MQ=255)
                           tGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCgg                                                                                                                  >  2:1333394/1‑151 (MQ=255)
                           tGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCg                                                                                                                   <  2:1066886/150‑1 (MQ=255)
                           tGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCg                                                                                                                   <  2:866373/150‑1 (MQ=255)
                            gAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCgg                                                                                                                  <  1:434523/150‑1 (MQ=255)
                             aaCAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGACGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCg                                                                                                                   >  1:309418/1‑148 (MQ=255)
                             aaCAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCgg                                                                                                                  >  1:49448/1‑149 (MQ=255)
                             aaCAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGaa                                                                                                                <  2:303903/151‑1 (MQ=255)
                             aaCAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGaa                                                                                                                <  2:202537/151‑1 (MQ=255)
                              aCAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAAt                                                                                                               <  1:579068/151‑1 (MQ=255)
                                   cTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGCTTAATGCAGCGTGTGGCCCGGAATCAtt                                                                                                           >  2:73615/1‑150 (MQ=255)
                                       aCCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCa                                                                                                      <  2:1302277/151‑1 (MQ=255)
                                        ccGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCa                                                                                                      <  1:47228/150‑1 (MQ=255)
                                        ccGAACCGCAGACGGTCTGGTCTTGCTGCATTCATTGGGCTGATTGCGGGCGGTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATc                                                                                                       >  1:805628/1‑149 (MQ=255)
                                          gAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAggg                                                                                                   >  1:538239/1‑151 (MQ=255)
                                            aCCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCt                                                                                                 <  1:130849/151‑1 (MQ=255)
                                             ccGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTg                                                                                                >  2:224550/1‑151 (MQ=255)
                                              cGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGc                                                                                               <  2:861699/151‑1 (MQ=255)
                                               gCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑CGTGGGGGGCTGAAGTTCCATTGCCGCCCCGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAACCATTAATCAGGTCTGc                                                                                               >  1:829956/1‑150 (MQ=255)
                                                   aCGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGcc                                                                                          >  1:334376/1‑151 (MQ=255)
                                                        aTGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAat                                                                                      <  2:863025/150‑1 (MQ=255)
                                                         tGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAattg                                                                                    >  2:966596/1‑151 (MQ=255)
                                                          ggTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAattga                                                                                   <  1:102977/151‑1 (MQ=255)
                                                               tgctgcATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCt                                                                              >  2:1043284/1‑151 (MQ=255)
                                                                gctgcATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTc                                                                             <  2:683845/151‑1 (MQ=255)
                                                                gctgcATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTc                                                                             <  2:572501/151‑1 (MQ=255)
                                                                     attcattGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTcgc                                                                           <  1:738625/148‑1 (MQ=255)
                                                                     attcattGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTcgc                                                                           >  2:739250/1‑148 (MQ=255)
                                                                     attcattGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTa                                                                        <  1:168494/151‑1 (MQ=255)
                                                                         attGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAAttt                                                                    >  1:1389186/1‑151 (MQ=255)
                                                                         attGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAAttt                                                                    >  1:506535/1‑151 (MQ=255)
                                                                                tGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAAt                                                             >  2:166633/1‑151 (MQ=255)
                                                                                  tttGCTGGCGTTGGTTCAGCACTAGTGC‑AGTGGGGCGATGACGTTGCATTGCCGCCACGTAGCCC‑GTGGTTATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCGGCTGGCGAAATTGATTGCTCGCGTAATTTTCTCAAt                                                             >  1:485336/2‑149 (MQ=255)
                                                                                  aTTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATcc                                                           >  2:968443/1‑151 (MQ=255)
                                                                                    tGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCAc                                                         >  2:510608/1‑151 (MQ=255)
                                                                                     gCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCAcc                                                        >  2:423935/1‑151 (MQ=255)
                                                                                     gCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCAcc                                                        <  2:844759/151‑1 (MQ=255)
                                                                                      cTGCCGTTCTTTCCGCATTCGTGAGAGT‑GGTGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGTCTATGTTTAATTCAGCGTGTGGCCCTTAATCATTCATCAGGGCTGCAGGACAAATTGATTGCTCGCGTAATTTGCTAAATCCAc                                                         <  1:6094/149‑1 (MQ=255)
                                                                                         gCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGtat                                                    >  2:100833/1‑151 (MQ=255)
                                                                                           gttgttTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTatat                                                  >  1:58699/1‑151 (MQ=255)
                                                                                             tgttTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATAttt                                                >  2:236813/1‑151 (MQ=255)
                                                                                              gttTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATAttt                                                <  2:513040/150‑1 (MQ=255)
                                                                                                  ccGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTttctt                                           >  2:884857/1‑151 (MQ=255)
                                                                                                      aTTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCg                                         >  2:13415/1‑149 (MQ=255)
                                                                                                      aTTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCg                                         >  1:1080880/1‑149 (MQ=255)
                                                                                                       ttCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCg                                         >  2:1280311/1‑148 (MQ=255)
                                                                                                        tCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGAc                                     >  2:584849/1‑151 (MQ=255)
                                                                                                            gA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGt                                 <  2:178260/151‑1 (MQ=255)
                                                                                                                gTGGGGGGCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTgg                              <  1:752567/151‑1 (MQ=255)
                                                                                                                       gCTGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGa                        >  1:370581/1‑150 (MQ=255)
                                                                                                                         tGAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGAc                       >  1:760775/1‑149 (MQ=255)
                                                                                                                          gAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGAc                       >  2:40773/1‑148 (MQ=255)
                                                                                                                          gAAGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACa                      <  1:42283/149‑1 (MQ=255)
                                                                                                                           aaGTTCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGAt                   <  2:873048/151‑1 (MQ=255)
                                                                                                                           aaGATCCATTTCCGCCACGCAG‑CCGGTGGATATGTTTAATGCAACGTGTGGCCCGGAATCATTAATCACGGCTGTAGGCCAAACTGATTGCTCGCGTAATTTTCTCATTACACCGTATATTTTTCTTCGGGACTGGTCGGCGCTg                         >  1:821090/1‑145 (MQ=255)
                                                                                                                               tCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATccc                >  1:1093582/1‑150 (MQ=255)
                                                                                                                               tCCATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATccc                >  1:436788/1‑150 (MQ=255)
                                                                                                                                 cATTGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATCCCAAt             >  2:885506/1‑151 (MQ=255)
                                                                                                                                    tGCCGCCACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATCCCAATGCg          >  2:1075333/1‑151 (MQ=255)
                                                                                                                                     gccgccACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATCCCAATGCgg         <  1:740046/151‑1 (MQ=255)
                                                                                                                                          cACGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATCCCAATGCGGCTGtt    >  2:583708/1‑151 (MQ=255)
                                                                                                                                           aCGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATCCCAATGCGGCTGttt   <  1:988850/151‑1 (MQ=255)
                                                                                                                                            cGTAG‑CCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATCCCAATGCGGCTGTTta  <  1:1234621/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TACAAGCAGGGTAGATATGAATATATTTGAACAAACTCCACCGAACCGCAGACGTTATGGTCTTGCTGCATTCATTGGGCTGATTGCTGGCGTTGTTTCCGCATTCGTGA‑AGTGGGGGGCTGAAGTTCCATTGCCGCCACGTAGCCCGGTGGATATGTTTAATGCAGCGTGTGGCCCGGAATCATTAATCAGGGCTGCAGGCCAAATTGATTGCTCGCGTAATTTTCTCAATCCACCGTATATTTTTCTTCGAGACTGGTTGGGGCTGACAGATCCCAATGCGGCTGTTTA  >  NC_000913/302975‑303265

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