Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,622,796 T→G L1202R (CTG→CGG)  rhsB → Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,622,7960TG97.7% 219.0 / ‑1.2 86L1202R (CTG→CGG) rhsBRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (50/34);  total (51/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

TGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGAT  >  NC_000913/3622701‑3622891
                                                                                               |                                                                                               
tgatgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTAt                                                                                             >  1:819148/1‑100 (MQ=25)
  atgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTAtt                                                                                            >  1:691771/1‑99 (MQ=17)
  atgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTa                                                                                           >  2:874582/1‑100 (MQ=18)
  atgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTa                                                                                           <  1:1891573/100‑1 (MQ=18)
  atgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTa                                                                                           <  1:1754760/100‑1 (MQ=18)
    gaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACa                                                                                         >  1:440461/1‑100 (MQ=18)
    gaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACa                                                                                         >  1:2161355/1‑100 (MQ=18)
    gaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACa                                                                                         >  1:286084/1‑100 (MQ=18)
       tGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAAcc                                                                                      <  1:774643/100‑1 (MQ=18)
       tGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAAcc                                                                                      >  1:2240773/1‑100 (MQ=18)
        ggggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGTCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCg                                                                                     >  2:1274565/1‑100 (MQ=255)
        ggggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCg                                                                                     >  2:915615/1‑100 (MQ=18)
         gggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCAGGGGCAGCAGTCTGATGAGGAGTACGGCCGGTATTACAACCGc                                                                                    >  1:1952447/1‑100 (MQ=11)
            cAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCAc                                                                                 <  1:615862/100‑1 (MQ=18)
             aaCCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCAcc                                                                                <  1:1322504/100‑1 (MQ=18)
             aaCCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCAcc                                                                                >  1:2215330/1‑100 (MQ=18)
              aCCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCg                                                                               >  2:556787/1‑100 (MQ=18)
                ctgctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCt                                                                             >  2:1934236/1‑100 (MQ=18)
                 tgctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCta                                                                            <  2:1680153/100‑1 (MQ=18)
                 tgctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCta                                                                            >  1:159826/1‑100 (MQ=18)
                 tgctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCta                                                                            >  2:1654864/1‑100 (MQ=18)
                  gctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtat                                                                           >  2:1814076/1‑100 (MQ=17)
                  gctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtat                                                                           >  2:339742/1‑100 (MQ=18)
                   ctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtatt                                                                          <  1:2077415/100‑1 (MQ=18)
                    tgaatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtatta                                                                         >  1:2152227/1‑100 (MQ=18)
                      aatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATg                                                                       <  2:286082/100‑1 (MQ=18)
                       atgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGa                                                                      >  2:2257246/1‑100 (MQ=18)
                        tgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAc                                                                     >  2:105364/1‑100 (MQ=18)
                        tgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAc                                                                     <  2:2161348/100‑1 (MQ=18)
                        tgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAc                                                                     >  1:783477/1‑100 (MQ=18)
                        tgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAc                                                                     >  1:261619/1‑100 (MQ=18)
                           ngagaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCg                                                                  >  2:1740439/2‑100 (MQ=18)
                           agagAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCg                                                                  >  2:919500/1‑100 (MQ=18)
                             agaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCt                                                                <  2:1284563/100‑1 (MQ=17)
                              gaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTg                                                               >  2:727058/1‑100 (MQ=18)
                               aaCCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGc                                                              <  2:634977/100‑1 (MQ=18)
                               aaCCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGc                                                              >  2:646806/1‑100 (MQ=18)
                               aaCCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGc                                                              <  2:1378936/100‑1 (MQ=18)
                               aaCCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGc                                                              <  2:1375343/100‑1 (MQ=18)
                                 cccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAg                                                            >  1:99040/1‑100 (MQ=18)
                                  ccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAgg                                                           <  2:819146/100‑1 (MQ=18)
                                   cGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAggg                                                          <  1:1482118/100‑1 (MQ=18)
                                      aTCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCg                                                       <  2:535024/100‑1 (MQ=17)
                                      aTCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCg                                                       <  1:1343672/100‑1 (MQ=18)
                                       tCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGa                                                      >  2:140959/1‑100 (MQ=18)
                                       tCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGa                                                      <  2:686404/100‑1 (MQ=18)
                                        cAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGat                                                     <  2:988969/100‑1 (MQ=18)
                                        cAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGat                                                     >  1:1158392/1‑100 (MQ=18)
                                          gCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatat                                                   <  1:2257254/100‑1 (MQ=18)
                                          gCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatat                                                   >  1:2228139/1‑100 (MQ=18)
                                            tGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatatat                                                 >  1:1171676/1‑100 (MQ=17)
                                             gcagcaGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATc                                                >  2:1239436/1‑100 (MQ=18)
                                             gcagcaGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATc                                                >  2:1239231/1‑100 (MQ=18)
                                                    cTTATCCGGCTGCC‑GGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGa                                        >  1:1930667/1‑100 (MQ=255)
                                                        tCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATcc                                     <  1:1654869/100‑1 (MQ=18)
                                                        tCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCAGCAGGGGCGATATATCACTCAGGATcc                                     >  2:973921/1‑100 (MQ=12)
                                                         ccGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCg                                    >  1:2051934/1‑100 (MQ=17)
                                                          ctgCTGCCCCGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGa                                   >  1:250554/3‑100 (MQ=11)
                                                          cGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGa                                   >  2:749/1‑100 (MQ=17)
                                                           ggCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGGTCCGAt                                  >  1:2021386/1‑100 (MQ=9)
                                                           ggCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGAt                                  <  2:99038/100‑1 (MQ=14)
                                                           ggCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGAt                                  >  1:2021286/1‑100 (MQ=14)
                                                            gCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGAtt                                 <  1:919502/100‑1 (MQ=11)
                                                               gCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGa                              >  1:1366440/1‑100 (MQ=11)
                                                                ccGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGa                              >  1:1202256/1‑99 (MQ=11)
                                                                  ggagcagcaGTATGATGAGTAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCGGCAGGGGCGATATATCACTCAGGATCCGATTGGACt                            <  2:202189/96‑1 (MQ=11)
                                                                      cagcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgg                       >  2:1290784/1‑100 (MQ=25)
                                                                        gcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggg                     >  2:28231/1‑100 (MQ=255)
                                                                        gcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggg                     <  1:339744/100‑1 (MQ=255)
                                                                           gTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGa                  >  1:2107036/1‑100 (MQ=25)
                                                                             atgatgAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATg                <  1:2073086/100‑1 (MQ=255)
                                                                             atgatgAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATg                <  2:182882/100‑1 (MQ=255)
                                                                             atgatgAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATg                <  1:45755/100‑1 (MQ=255)
                                                                             atgatgAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATg                <  1:466900/100‑1 (MQ=255)
                                                                               gatgaGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGa              >  1:1911977/1‑100 (MQ=38)
                                                                                  gaggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAc            <  1:556788/99‑1 (MQ=38)
                                                                                    ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg         <  2:2021379/100‑1 (MQ=255)
                                                                                    ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg         <  2:2021279/100‑1 (MQ=255)
                                                                                     gagTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGCGGATGGAACCtat        <  1:1587830/100‑3 (MQ=255)
                                                                                     gagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGt        <  2:716514/100‑1 (MQ=255)
                                                                                     gagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGt        >  1:1744058/1‑100 (MQ=255)
                                                                                       gTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTAt      >  2:713593/1‑100 (MQ=255)
                                                                                         ccggccggTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGTATGGAACCTGTATgg    >  1:2296136/1‑100 (MQ=255)
                                                                                         ccggccggTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATgg    >  1:519810/1‑100 (MQ=255)
                                                                                         ccggccggTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATgg    >  2:133843/1‑100 (MQ=255)
                                                                                         ccggccggTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGCATGGAACCTGTATgg    >  1:2296264/1‑100 (MQ=255)
                                                                                          cGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGa   <  2:2215322/100‑1 (MQ=255)
                                                                                           ggCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGat  >  2:101231/1‑100 (MQ=255)
                                                                                           ggCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGat  >  1:783511/1‑100 (MQ=255)
                                                                                           ggCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGat  >  2:100855/1‑100 (MQ=255)
                                                                                               |                                                                                               
TGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGAT  >  NC_000913/3622701‑3622891

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