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.

ATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGAT  >  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)
                         tgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAc                                                                     >  2:539290/1‑100 (MQ=18)
                         tgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCGATTATGa                                                                       >  2:218846/1‑99 (MQ=12)
                           aaGAGAACCCGCTTCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGACGAGTCCGGCCGGTATTACATCCGCCACCGCTATTATGAc                                                                     <  2:2810852/98‑1 (MQ=11)
                           aaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAccc                                                                   >  2:939219/1‑100 (MQ=17)
                           aaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAccc                                                                   >  1:2241971/1‑100 (MQ=18)
                           aaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAccc                                                                   >  1:1287707/1‑100 (MQ=18)
                           aaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAATCCGGCCGGTATTACAACCGACACCGCTAAAATGAccc                                                                   >  1:1002077/1‑100 (MQ=11)
                             aagaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGc                                                                 >  2:1149540/2‑100 (MQ=17)
                               gaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtattatagatcgga                                                                 >  1:64883/1‑91 (MQ=11)
                               gaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTg                                                               >  1:2850447/1‑100 (MQ=18)
                                  cccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAg                                                            >  2:2141623/1‑100 (MQ=17)
                                  cccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAg                                                            >  1:1055481/1‑100 (MQ=18)
                                  cccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAg                                                            >  1:2218480/1‑100 (MQ=18)
                                   ccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAgg                                                           >  2:1400799/1‑100 (MQ=17)
                                   ccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAgg                                                           >  1:2000332/1‑100 (MQ=18)
                                   ccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCAGCTGCAg                                                            >  1:1552025/1‑99 (MQ=16)
                                     gCATAAGCTGCAGCAGGTTATTCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAgggg                                                         <  1:426686/100‑1 (MQ=11)
                                     atctCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAgggg                                                         <  1:1702765/97‑1 (MQ=255)
                                      cATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGc                                                        <  1:403274/100‑1 (MQ=17)
                                       aTCAGCTTCAGCAGCTTATCCGGCTGCCGGGGCAGCATTATTATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCg                                                       <  1:2538343/100‑1 (MQ=255)
                                       aTCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCg                                                       <  2:596876/100‑1 (MQ=17)
                                       aTCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCg                                                       <  1:439827/100‑1 (MQ=17)
                                         cAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGAGat                                                     >  2:1702745/1‑100 (MQ=9)
                                         cAGCTGCAGCAGCGTATCCGGCTGCCGGGCCAGCATTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGat                                                     <  2:1486952/100‑1 (MQ=11)
                                          aGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata                                                    >  2:2033416/1‑100 (MQ=17)
                                          aGCTGCAGCAGCTCATCCGGCTGCCTGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGGAGGGGCGa                                                      <  2:1751721/98‑1 (MQ=11)
                                            cTGCAGCAGCTTATGCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTTTTATGACCCGCTGCAGGGGCGata                                                    <  2:864754/98‑1 (MQ=255)
                                             tGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatata                                                  <  2:2808715/99‑1 (MQ=18)
                                              gcagcaGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATc                                                >  1:2983987/1‑100 (MQ=18)
                                              gcagcaGCGTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatatatt                                                <  2:2864468/100‑2 (MQ=12)
                                               cagcagCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCa                                               <  1:2496521/100‑1 (MQ=17)
                                                 gcagcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACt                                             >  1:2737324/1‑100 (MQ=18)
                                                   agcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAGCCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTc                                            <  1:2150938/99‑1 (MQ=12)
                                                   agcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCa                                           >  1:2641986/1‑100 (MQ=17)
                                                    gcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGATATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGataagatcggaag                                          >  1:1749404/1‑90 (MQ=11)
                                                     cTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAgg                                         <  2:176461/100‑1 (MQ=18)
                                                        aTCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATc                                      >  1:2700891/1‑100 (MQ=17)
                                                         tCCGGCTGCCGGGGCAGGAGTATGATGAGAGGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGTGGCGATATATCACTCAGGATcc                                     <  2:2309481/100‑1 (MQ=11)
                                                                 ccGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGAc                             >  1:2815111/1‑100 (MQ=14)
                                                                 ccGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCAATCAGGATACGATTGGAc                             >  1:1174936/1‑100 (MQ=255)
                                                                  cGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACt                            >  1:3021221/1‑100 (MQ=17)
                                                                      gcagcaGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCACTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAg                        <  1:2101913/100‑1 (MQ=255)
                                                                       cagcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgg                       >  1:1967951/1‑100 (MQ=25)
                                                                         gcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggg                     >  1:1636182/1‑100 (MQ=255)
                                                                          cagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAggggg                    >  1:2905372/1‑100 (MQ=25)
                                                                          cagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAggggg                    <  1:1947038/100‑1 (MQ=38)
                                                                          cagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGATGCAGGAGCGATAAATCACTCAGGATCCGATTGGACTGAAggggg                    >  1:112231/1‑100 (MQ=11)
                                                                           agTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggggg                   >  2:2103644/1‑100 (MQ=38)
                                                                           agTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggggg                   <  2:2983954/100‑1 (MQ=255)
                                                                           agTATGATGAGGAGTCCGGCCGGTATGACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggggg                   <  2:2861624/100‑1 (MQ=255)
                                                                             tATGATGAGGAGTACGGCCGGTATTACAACCACCACCGCTATGATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGAt                 >  1:2046761/1‑100 (MQ=11)
                                                                              atgatgAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATg                <  2:1642834/100‑1 (MQ=255)
                                                                               tgatgaGGGGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATgg               >  2:798818/1‑100 (MQ=255)
                                                                               tgatgaGGGGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATg                >  2:2903789/1‑99 (MQ=255)
                                                                               tgatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATgg               >  1:2406943/1‑100 (MQ=255)
                                                                                gatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGa              >  1:987480/1‑100 (MQ=255)
                                                                                gatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGa              >  2:1152561/1‑100 (MQ=25)
                                                                                  tgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGaa             >  1:687969/1‑99 (MQ=255)
                                                                                   gaggagTCCTGCCGGTATTACAACCTCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAcc           <  2:876306/100‑1 (MQ=255)
                                                                                    aggagTCCGGCCGGTATTACAACCGTCACCGCTATTATTTTTCGTTTTAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAcc           <  1:876320/99‑1 (MQ=14)
                                                                                    aggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCt          >  2:1950718/1‑100 (MQ=255)
                                                                                    aggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCt          >  1:2839898/1‑100 (MQ=255)
                                                                                     ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg         >  1:1169403/1‑100 (MQ=255)
                                                                                     ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg         >  1:2486970/1‑100 (MQ=255)
                                                                                     ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg         <  2:2479990/100‑1 (MQ=255)
                                                                                     ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg         <  1:1222825/100‑1 (MQ=255)
                                                                                       agTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCCGTTGCAGGGGCAATAATTAACTAAGAACCCGATGGGCTTGAGggggggtggaca             >  2:785326/1‑89 (MQ=2)
                                                                                       agTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTa       >  1:451216/1‑100 (MQ=255)
                                                                                       agTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTa       >  1:1063794/1‑100 (MQ=255)
                                                                                        gTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTAt      >  1:1024506/1‑100 (MQ=255)
                                                                                        gTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTAt      >  2:923660/1‑100 (MQ=255)
                                                                                          ccggccggTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATgg    >  1:2285821/1‑100 (MQ=255)
                                                                                           nggCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGa   >  2:2620187/2‑100 (MQ=255)
                                                                                           cGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGGTTGGACTGAAGGGGGGATGGAACCTGTATGGa   <  1:2887248/100‑1 (MQ=255)
                                                                                                |                                                                                                
ATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGAT  >  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: