Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,560,233 G→A Q150* (CAG→TAG)  yjiH ← nucleoside recognition pore and gate family putative inner membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,560,2330GA98.5% 249.0 / NA 65Q150* (CAG→TAG) yjiHnucleoside recognition pore and gate family putative inner membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (38/26);  minor base T (0/1);  total (38/27)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCAT  >  NC_000913/4560085‑4560379
                                                                                                                                                    |                                                                                                                                                  
aGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTacg                                                                                                                                                  <  2:756755/151‑1 (MQ=255)
 ggATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTacgc                                                                                                                                                 >  2:448707/1‑151 (MQ=255)
 ggATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTacgc                                                                                                                                                 <  1:716585/151‑1 (MQ=255)
   aTGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGcc                                                                                                                                                >  2:810945/1‑150 (MQ=255)
      aCGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGct                                                                                                                                            <  1:1101561/151‑1 (MQ=255)
       cGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGctc                                                                                                                                           >  1:1143667/1‑151 (MQ=255)
        cgCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGctct                                                                                                                                          >  1:368192/2‑151 (MQ=255)
           tAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTg                                                                                                                                       >  2:954483/1‑151 (MQ=255)
            aaaGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTgg                                                                                                                                      >  2:437394/1‑151 (MQ=255)
            aaaGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTgg                                                                                                                                      >  1:1106324/1‑151 (MQ=255)
            aaaGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTgg                                                                                                                                      <  2:590365/151‑1 (MQ=255)
            aaaGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGGAGTTGGTGATGCTGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCtt                                                                                                                                         >  1:1089294/1‑148 (MQ=255)
                  aCGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATAc                                                                                                                                >  1:652222/1‑151 (MQ=255)
                  aCGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATAc                                                                                                                                <  1:817018/151‑1 (MQ=255)
                      ttACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACcggc                                                                                                                            >  1:249049/1‑151 (MQ=255)
                       tACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACcggcg                                                                                                                           >  2:1131854/1‑151 (MQ=255)
                       tACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACcggc                                                                                                                            <  1:770489/150‑1 (MQ=255)
                        aCGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACcggcgg                                                                                                                          <  2:276768/151‑1 (MQ=255)
                         cGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACcggcggc                                                                                                                         <  2:733591/151‑1 (MQ=255)
                         cGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACcggcggc                                                                                                                         <  1:683126/151‑1 (MQ=255)
                           gACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACcggcggcgg                                                                                                                       >  1:754859/1‑151 (MQ=255)
                                gCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCg                                                                                                                  <  1:496814/151‑1 (MQ=255)
                                       aaGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGtttt                                                                                                           <  1:1059156/151‑1 (MQ=255)
                                       aaGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGtttt                                                                                                           >  1:203475/1‑151 (MQ=255)
                                                acCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGtt                                                                                                    >  2:327711/1‑149 (MQ=255)
                                                  cgcgACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCg                                                                                                >  1:842444/1‑151 (MQ=255)
                                                   gcgACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGa                                                                                               <  1:844809/151‑1 (MQ=255)
                                                     gACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATc                                                                                             >  1:1110868/1‑151 (MQ=255)
                                                      aCGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCa                                                                                            <  2:842458/151‑1 (MQ=255)
                                                       cGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAg                                                                                           <  2:476409/151‑1 (MQ=255)
                                                        gCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAgc                                                                                          >  1:828996/1‑151 (MQ=255)
                                                           ggAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGcc                                                                                       >  1:828191/1‑151 (MQ=255)
                                                           ggAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGcc                                                                                       >  1:814684/1‑151 (MQ=255)
                                                             aaGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGcc                                                                                       >  1:14184/1‑149 (MQ=255)
                                                                 gaAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCg                                                                                   >  1:1099537/1‑149 (MQ=255)
                                                                 gaAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAg                                                                                 <  2:1053577/151‑1 (MQ=255)
                                                                  aaGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAg                                                                                 >  2:115861/1‑150 (MQ=255)
                                                                  aaGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGc                                                                                >  2:451030/1‑151 (MQ=255)
                                                                      agtTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGc                                                                                >  2:526538/1‑147 (MQ=255)
                                                                          ggTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCt                                                                        <  1:426442/151‑1 (MQ=255)
                                                                          ggTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCt                                                                        >  2:1064188/1‑151 (MQ=255)
                                                                           gTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCAGTAATTTCTCCTTCCTTCGCCAGCTCTTTGGTCATACCGGCGGCGGCCTCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCCGATcc                                                                          <  1:800699/148‑1 (MQ=255)
                                                                              atgatgGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGa                                                                     >  1:312744/1‑150 (MQ=255)
                                                                                gatgGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAAt                                                                   <  2:19268/150‑1 (MQ=255)
                                                                                        tACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCcaaca                                                           >  1:46338/1‑150 (MQ=255)
                                                                                         aCCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCaacaac                                                         >  2:716286/1‑151 (MQ=255)
                                                                                          ccGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACg                                                        >  1:638002/1‑151 (MQ=255)
                                                                                               gggttGATATGCGGCGAAAATCAATTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAAAGGtttt                                                   >  2:916092/4‑151 (MQ=255)
                                                                                                     atatGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGcaca                                             <  1:162963/151‑1 (MQ=255)
                                                                                                     atatGCCGCGAAAATCACTTTATCGCGGTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGGGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGcac                                              >  1:505529/1‑150 (MQ=255)
                                                                                                       atGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAg                                            <  1:409454/150‑1 (MQ=255)
                                                                                                            gcgAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCa                                      <  1:94694/151‑1 (MQ=255)
                                                                                                             cgAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCAt                                     >  1:1070082/1‑151 (MQ=255)
                                                                                                             cgAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCAt                                     >  2:50997/1‑151 (MQ=255)
                                                                                                                        ttATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACtgctgc                           <  2:312747/150‑1 (MQ=255)
                                                                                                                        ttATCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTgcg                          >  2:646738/1‑151 (MQ=255)
                                                                                                                           tCGCGTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTgcgcgg                       >  2:122309/1‑151 (MQ=255)
                                                                                                                            cgcgTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTgcgcggc                      <  2:662309/151‑1 (MQ=255)
                                                                                                                            cgcgTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTgcgcggc                      >  1:816557/1‑151 (MQ=255)
                                                                                                                             gcgTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTgcgcggcg                     >  1:864848/1‑151 (MQ=255)
                                                                                                                              cgTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTgcgcggcgc                    <  1:559423/151‑1 (MQ=255)
                                                                                                                              cgTTCGGTAATTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTgcgcggcgc                    <  1:130852/151‑1 (MQ=255)
                                                                                                                                      aaTTTCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATcca            <  2:1153638/151‑1 (MQ=255)
                                                                                                                                         ttCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATccacc          >  2:457688/1‑150 (MQ=255)
                                                                                                                                          tCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATccacca         <  1:984118/150‑1 (MQ=255)
                                                                                                                                          tCACCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATccacca         <  1:1153720/150‑1 (MQ=255)
                                                                                                                                            aCCTTCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCaa        <  2:998528/149‑1 (MQ=255)
                                                                                                                                                tCCTACGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCAt  >  2:916520/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGGATGACGGCTAAAGGTACGATTACGGACGTGCCCAGAAAGGCAAACACCGCGACGCCGGAAGAGAAGTAGTTGGTGATGATGGCGCTACCGCTGGTTTGATATGCGGCGAAAATCACTTTATCGCGTTCGGTAATTTCACCTTCCTGCGCCAGCTCTTTGGTCATACCGGCGGCGGCATCGGTGTTTTGCAGGTTAGCGATCAGCGCCAGCGAGCAGATCCCTGGAATGCCCAACAACGGTTTTAGCACAGGTGTCATTAACTGCTGCGCGGCGCGTAATCCACCAAGACCAT  >  NC_000913/4560085‑4560379

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