Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,622,796 T→G 100% 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,7960TG95.8% 338.1 / 6.2 141L1202R (CTG→CGG) rhsBRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base G (81/53);  total (84/57)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

Only 100 of 143 total aligned reads displayed.

ATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGAT  >  NC_000913/3622700‑3622891
                                                                                                |                                                                                                
atgatgAATGGGGCAACCTGCTGAATGAAGAGAGCCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGGCCGGCCGGt                                                                                                >  2:1998772/1‑99 (MQ=255)
 tgatgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTAt                                                                                              >  1:2111038/1‑100 (MQ=25)
 tgatgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTAt                                                                                              >  1:2794643/1‑100 (MQ=25)
  gatgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTAtt                                                                                             <  2:1828378/100‑1 (MQ=18)
    tgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTAc                                                                                           >  2:225148/1‑100 (MQ=18)
    tgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAAGAAACCGACCGGTATTa                                                                                            >  1:67336/1‑99 (MQ=11)
     gaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACa                                                                                          <  2:2567218/100‑1 (MQ=18)
     gaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACa                                                                                          >  2:1950858/1‑100 (MQ=18)
        tGGGGCAACCTGCTGAATGAAGAGAACTCGTATCAGTTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAAcc                                                                                       <  1:963956/100‑1 (MQ=11)
        tGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAAcc                                                                                       <  1:1037240/100‑1 (MQ=17)
          gggCAACCTGATGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCAGGGCAGCAGTATGACGAGCAGTCCGGCCGGTCTTACAccc                                                                                       >  2:2869085/1‑98 (MQ=11)
           ggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGcc                                                                                    <  2:1531476/100‑1 (MQ=18)
           ggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGcc                                                                                    <  2:530613/100‑1 (MQ=17)
            gCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCa                                                                                   >  1:2294589/1‑100 (MQ=18)
              aaCCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCAcc                                                                                 <  2:2254009/100‑1 (MQ=18)
                ccTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGc                                                                               <  2:2509045/100‑1 (MQ=18)
                ccTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGc                                                                               >  1:2509074/1‑100 (MQ=18)
                   gctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtat                                                                            >  2:2644616/1‑100 (MQ=17)
                   gctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtat                                                                            >  2:952415/1‑100 (MQ=18)
                    ctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtatt                                                                           <  1:874012/100‑1 (MQ=17)
                    ctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtatt                                                                           <  1:2933882/100‑1 (MQ=18)
                     tgaatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtatta                                                                          >  1:2025574/1‑100 (MQ=18)
                      gaatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtattat                                                                         >  2:2925277/1‑100 (MQ=17)
                      gaatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtattat                                                                         >  1:1099523/1‑100 (MQ=17)
                      gaatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtattat                                                                         <  2:2395719/100‑1 (MQ=18)
                        atgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGa                                                                       >  1:1368437/1‑100 (MQ=18)
                         tgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑c                                                                     >  2:539290/1‑100 (MQ=18)
                         tgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCGATTATGa                                                                       >  2:218846/1‑99 (MQ=12)
                           aaGAGAACCCGCTTCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGACGAGTCCGGCCGGTATTACATCCGCCACCGCTATTATGA‑c                                                                     <  2:2810852/98‑1 (MQ=11)
                           aaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑ccc                                                                   >  2:939219/1‑100 (MQ=17)
                           aaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑ccc                                                                   >  1:2241971/1‑100 (MQ=18)
                           aaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑ccc                                                                   >  1:1287707/1‑100 (MQ=18)
                           aaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAATCCGGCCGGTATTACAACCGACACCGCTAAAATGA‑ccc                                                                   >  1:1002077/1‑100 (MQ=11)
                             aagaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGc                                                                 >  2:1149540/2‑100 (MQ=17)
                               gaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtattatagatcgga                                                                 >  1:64883/1‑91 (MQ=11)
                               gaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTg                                                               >  1:2850447/1‑100 (MQ=18)
                                  cccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAg                                                            >  2:2141623/1‑100 (MQ=17)
                                  cccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAg                                                            >  1:1055481/1‑100 (MQ=18)
                                  cccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAg                                                            >  1:2218480/1‑100 (MQ=18)
                                   ccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAgg                                                           >  2:1400799/1‑100 (MQ=17)
                                   ccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAgg                                                           >  1:2000332/1‑100 (MQ=18)
                                   ccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCAGCTGCAg                                                            >  1:1552025/1‑99 (MQ=16)
                                     gCATAAGCTGCAGCAGGTTATTCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAgggg                                                         <  1:426686/100‑1 (MQ=11)
                                     atctCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAgggg                                                         <  1:1702765/97‑1 (MQ=255)
                                      cATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGc                                                        <  1:403274/100‑1 (MQ=17)
                                       aTCAGCTTCAGCAGCTTATCCGGCTGCCGGGGCAGCATTATTATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCg                                                       <  1:2538343/100‑1 (MQ=255)
                                       aTCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCg                                                       <  2:596876/100‑1 (MQ=17)
                                       aTCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCg                                                       <  1:439827/100‑1 (MQ=17)
                                         cAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGAGat                                                     >  2:1702745/1‑100 (MQ=9)
                                         cAGCTGCAGCAGCGTATCCGGCTGCCGGGCCAGCATTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGat                                                     <  2:1486952/100‑1 (MQ=11)
                                          aGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGata                                                    >  2:2033416/1‑100 (MQ=17)
                                          aGCTGCAGCAGCTCATCCGGCTGCCTGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGGAGGGGCGa                                                      <  2:1751721/98‑1 (MQ=11)
                                            cTGCAGCAGCTTATGCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTTTTATGA‑CCCGCTGCAGGGGCGata                                                    <  2:864754/98‑1 (MQ=255)
                                             tGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGatata                                                  <  2:2808715/99‑1 (MQ=18)
                                              gcagcaGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATc                                                >  1:2983987/1‑100 (MQ=18)
                                              gcagcaGCGTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGatatatt                                                <  2:2864468/100‑2 (MQ=12)
                                               cagcagCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCa                                               <  1:2496521/100‑1 (MQ=17)
                                                 gcagcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACt                                             >  1:2737324/1‑100 (MQ=18)
                                                   agcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAGCCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTc                                            <  1:2150938/99‑1 (MQ=12)
                                                   agcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCa                                           >  1:2641986/1‑100 (MQ=17)
                                                    gcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGATATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGataagatcggaag                                          >  1:1749404/1‑90 (MQ=11)
                                                     cTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAgg                                         <  2:176461/100‑1 (MQ=18)
                                                        aTCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATc                                      >  1:2700891/1‑100 (MQ=17)
                                                         tCCGGCTGCCGGGGCAGGAGTATGATGAGAGGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGTGGCGATATATCACTCAGGATcc                                     <  2:2309481/100‑1 (MQ=11)
                                                                 ccGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGAc                             >  1:2815111/1‑100 (MQ=14)
                                                                 ccGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCAATCAGGATACGATTGGAc                             >  1:1174936/1‑100 (MQ=255)
                                                                  cGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACt                            >  1:3021221/1‑100 (MQ=17)
                                                                      gcagcaGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCACTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAg                        <  1:2101913/100‑1 (MQ=255)
                                                                       cagcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgg                       >  1:1967951/1‑100 (MQ=25)
                                                                         gcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggg                     >  1:1636182/1‑100 (MQ=255)
                                                                          cagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAggggg                    >  1:2905372/1‑100 (MQ=25)
                                                                          cagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAggggg                    <  1:1947038/100‑1 (MQ=38)
                                                                          cagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGATGCAGGAGCGATAAATCACTCAGGATCCGATTGGACTGAAggggg                    >  1:112231/1‑100 (MQ=11)
                                                                           agTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggggg                   >  2:2103644/1‑100 (MQ=38)
                                                                           agTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggggg                   <  2:2983954/100‑1 (MQ=255)
                                                                           agTATGATGAGGAGTCCGGCCGGTATGACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggggg                   <  2:2861624/100‑1 (MQ=255)
                                                                             tATGATGAGGAGTACGGCCGGTATTACAACCACCACCGCTATGATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGAt                 >  1:2046761/1‑100 (MQ=11)
                                                                              atgatgAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATg                <  2:1642834/100‑1 (MQ=255)
                                                                               tgatgaGGGGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATgg               >  2:798818/1‑100 (MQ=255)
                                                                               tgatgaGGGGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATg                >  2:2903789/1‑99 (MQ=255)
                                                                               tgatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATgg               >  1:2406943/1‑100 (MQ=255)
                                                                                gatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGa              >  1:987480/1‑100 (MQ=255)
                                                                                gatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGa              >  2:1152561/1‑100 (MQ=25)
                                                                                  tgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGaa             >  1:687969/1‑99 (MQ=255)
                                                                                   gaggagTCCTGCCGGTATTACAACCTCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAcc           <  2:876306/100‑1 (MQ=255)
                                                                                    aggagTCCGGCCGGTATTACAACCGTCACCGCTATTATTT‑TTCGTTTTAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAcc           <  1:876320/99‑1 (MQ=14)
                                                                                    aggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCt          >  2:1950718/1‑100 (MQ=255)
                                                                                    aggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCt          >  1:2839898/1‑100 (MQ=255)
                                                                                     ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg         >  1:1169403/1‑100 (MQ=255)
                                                                                     ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg         >  1:2486970/1‑100 (MQ=255)
                                                                                     ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg         <  2:2479990/100‑1 (MQ=255)
                                                                                     ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg         <  1:1222825/100‑1 (MQ=255)
                                                                                       agTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCCGTTGCAGGGGCAATAATTAACTAAGAACCCGATGGGCTTGAGggggggtggaca             >  2:785326/1‑89 (MQ=2)
                                                                                       agTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTa       >  1:451216/1‑100 (MQ=255)
                                                                                       agTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTa       >  1:1063794/1‑100 (MQ=255)
                                                                                        gTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTAt      >  1:1024506/1‑100 (MQ=255)
                                                                                        gTCCGGCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTAt      >  2:923660/1‑100 (MQ=255)
                                                                                          ccggccggTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATgg    >  1:2285821/1‑100 (MQ=255)
                                                                                           nggCCGGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGa   >  2:2620187/2‑100 (MQ=255)
                                                                                           cGGCCTGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGGTTGGACTGAAGGGGGGATGGAACCTGTATGGa   <  1:2887248/100‑1 (MQ=255)
                                                                                                |                                                                                                
ATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGA‑CCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGAT  >  NC_000913/3622700‑3622891

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