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,7960TG100.0% 263.3 / NA 98L1202R (CTG→CGG) rhsBRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (48/50);  total (48/50)

GAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATA  >  NC_000913/3622705‑3622892
                                                                                           |                                                                                                
naaTGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACa                                                                                          >  1:1828959/2‑100 (MQ=17)
gaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGTCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACa                                                                                          >  1:995059/1‑100 (MQ=255)
gaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACa                                                                                          <  2:860256/100‑1 (MQ=18)
gaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACa                                                                                          <  2:394890/100‑1 (MQ=18)
    ggggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCg                                                                                      <  1:447685/100‑1 (MQ=18)
        cAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCAc                                                                                  <  1:2055849/100‑1 (MQ=18)
         aaCCTGTTGAATGAAGAGATCCCGAATCAGCTGCAGCAGCTTATCCGGCTGCCGGGTCAGCAGTATGATGAGGAGTCTGGCCGGTATTACAACCGCCAcc                                                                                 <  2:81643/100‑1 (MQ=11)
           ccTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGc                                                                               <  1:1548906/100‑1 (MQ=18)
             tgctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCta                                                                             >  1:946715/1‑100 (MQ=17)
             tgctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAACAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCta                                                                             >  2:1235279/1‑100 (MQ=12)
               ctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtatt                                                                           <  1:1329273/100‑1 (MQ=18)
                tgaatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtatta                                                                          <  2:1986055/100‑1 (MQ=18)
                 gaatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtattat                                                                         >  2:508972/1‑100 (MQ=17)
                 gaatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCACtattat                                                                         >  2:1399209/1‑100 (MQ=255)
                   ctgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGa                                                                       <  1:247928/99‑1 (MQ=17)
                    tgaaGAGAGCCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAc                                                                      <  2:1164451/100‑1 (MQ=255)
                    tgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAc                                                                      <  1:2442844/100‑1 (MQ=18)
                     gaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGaag                                                                     >  2:247925/1‑98 (MQ=18)
                      aaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAccc                                                                    >  2:32740/1‑100 (MQ=18)
                        gagaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGTGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGc                                                                  <  2:753661/100‑1 (MQ=12)
                           aaCCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGc                                                               >  1:2226715/1‑100 (MQ=18)
                           aaCCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGc                                                               <  2:949316/100‑1 (MQ=17)
                           aaCCCGCATCAACTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGc                                                               <  1:1488100/100‑1 (MQ=12)
                            aCCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCa                                                              >  1:543298/1‑100 (MQ=17)
                             cccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAg                                                             >  1:965120/1‑100 (MQ=18)
                              ccGCATCAGCCGCCGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAgg                                                            <  2:452832/100‑1 (MQ=255)
                                gCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAgggg                                                          <  1:1527632/100‑1 (MQ=17)
                                gCATCAGCTGCAGCAGCTTATCCGGCTGCCGGCGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAgggg                                                          <  2:497125/100‑1 (MQ=12)
                                   tCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGa                                                       <  2:1509376/100‑1 (MQ=18)
                                    cAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGan                                                      <  2:1972597/100‑2 (MQ=18)
                                    cAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGac                                                      <  2:146264/100‑2 (MQ=18)
                                     aGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata                                                     <  2:181685/100‑1 (MQ=17)
                                     aGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata                                                     <  1:41304/100‑1 (MQ=17)
                                      gCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatat                                                    >  2:2377767/1‑100 (MQ=18)
                                      gCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatat                                                    >  2:951180/1‑100 (MQ=17)
                                      gCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGat                                                      >  2:2394395/1‑98 (MQ=17)
                                        tgctgcAGCTTATCCGGCTGCCGGGGCTGCAGTATGATGAGGGGTCCGGCCGGTATTGCTACCGCGAGCGCTATTATGAGGCGCTGCAGGGGCGat                                                      >  2:117199/1‑96 (MQ=11)
                                        tGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatatat                                                  >  1:1085710/1‑100 (MQ=18)
                                        tGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatatat                                                  >  2:831701/1‑100 (MQ=18)
                                         gcagcaGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATc                                                 <  1:186831/100‑1 (MQ=18)
                                           agcagcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCAc                                               <  2:839863/100‑1 (MQ=18)
                                           agcagcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCAc                                               <  1:1187804/100‑1 (MQ=18)
                                            gcagcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACt                                              >  1:1764531/1‑100 (MQ=18)
                                               gcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAg                                           <  1:567609/100‑1 (MQ=18)
                                                  tATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGAt                                        >  1:348954/1‑100 (MQ=18)
                                                  tATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGAt                                        >  1:2363865/1‑100 (MQ=18)
                                                   atactgCTGCCCCGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTTTTATGACCACCCCCAGGGCAGATATATCAccc                                             >  2:203381/6‑92 (MQ=2)
                                                   atacGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTAAGGa                                         >  1:848185/4‑98 (MQ=255)
                                                    tCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCTACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATcc                                      <  2:1424979/100‑1 (MQ=255)
                                                    tCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATcc                                      <  2:13290/100‑1 (MQ=18)
                                                    tCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATcc                                      <  1:1880376/100‑1 (MQ=18)
                                                    tCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATcc                                      <  2:514944/100‑1 (MQ=17)
                                                    tCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATcc                                      <  1:1140816/100‑1 (MQ=18)
                                                    tCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATc                                       >  1:1229400/1‑99 (MQ=17)
                                                      cGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGa                                    >  2:2147274/1‑100 (MQ=17)
                                                       ggCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGa                                    >  2:277857/1‑99 (MQ=14)
                                                        gCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGAtt                                  <  1:1346282/100‑1 (MQ=11)
                                                           gCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGa                               >  1:2319219/1‑100 (MQ=11)
                                                            ccGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGAc                              >  1:784518/1‑100 (MQ=14)
                                                            ccGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGAc                              <  2:2478169/100‑1 (MQ=14)
                                                            ccGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGAc                              >  1:1170093/1‑100 (MQ=14)
                                                             cGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACt                             >  2:184164/1‑100 (MQ=17)
                                                                ggCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGa                           <  2:638927/99‑1 (MQ=17)
                                                                ggCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGAGTGGACTGaa                          >  1:1526686/1‑100 (MQ=12)
                                                                 gcagcaGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAg                         <  2:1722229/100‑1 (MQ=25)
                                                                 gcagcaGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAg                         >  2:2342968/1‑100 (MQ=25)
                                                                 gcagcaGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAg                         >  1:279045/1‑100 (MQ=25)
                                                                    gcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggg                      >  1:1173953/1‑100 (MQ=255)
                                                                     cagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAggggg                     <  1:199352/100‑1 (MQ=38)
                                                                     cagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAggggg                     >  1:1979830/1‑100 (MQ=25)
                                                                      agTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggggg                    <  1:286295/100‑1 (MQ=255)
                                                                      agTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggggg                    <  1:1612080/100‑1 (MQ=255)
                                                                       gTATGGGGAGGAGTCCGGCCGGGATGACACCCGGCACCGCGATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggg                      <  2:2143386/97‑1 (MQ=11)
                                                                       gTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTACTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggggg                    <  1:1656848/99‑1 (MQ=255)
                                                                          tgatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGCTCCGATTGGACTGAAGGGGGGATgg                >  2:1823360/1‑100 (MQ=255)
                                                                          tgatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATgn                <  1:1562829/100‑2 (MQ=255)
                                                                           gatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGa               >  1:533407/1‑100 (MQ=255)
                                                                            atgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGaa              >  2:608912/1‑100 (MQ=255)
                                                                             tgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAc             >  2:2319476/1‑100 (MQ=255)
                                                                             tgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAc             <  1:2342983/100‑1 (MQ=255)
                                                                             tgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAc             >  1:2222397/1‑100 (MQ=255)
                                                                              gaggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAcc            <  2:1595422/100‑1 (MQ=255)
                                                                                ggatgccggccggTAGTACAAGCGCCACGGCTATTATGAGCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAcc            <  2:2032506/93‑1 (MQ=14)
                                                                                ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCt           >  2:912038/1‑99 (MQ=255)
                                                                                ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg          >  1:316276/1‑100 (MQ=255)
                                                                                ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg          >  1:1197601/1‑100 (MQ=255)
                                                                                ggagTCCGGCCGGGATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAcc            <  2:2126101/98‑1 (MQ=255)
                                                                                  agTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTa        >  1:371891/1‑100 (MQ=255)
                                                                                   gTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTAt       >  2:1160598/1‑100 (MQ=255)
                                                                                   gTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTAt       >  2:2468394/1‑100 (MQ=255)
                                                                                    tCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATg      >  1:765085/1‑100 (MQ=255)
                                                                                      cGGCCGGTGTTACGACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGa    <  2:1031379/100‑1 (MQ=255)
                                                                                      cGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGa    <  2:371887/100‑1 (MQ=255)
                                                                                      cGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGa    <  1:1679172/100‑1 (MQ=255)
                                                                                       ggCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGat   >  1:953652/1‑100 (MQ=255)
                                                                                       ggCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGat   <  1:2018204/100‑1 (MQ=255)
                                                                                       ggCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGat   >  1:2418484/1‑100 (MQ=255)
                                                                                       ggCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGat   <  2:2430021/100‑1 (MQ=255)
                                                                                        gcttGTATTATATCCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATAAGTCAGGCTCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGat   <  1:1380758/96‑1 (MQ=17)
                                                                                        gCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGata  <  1:2178038/100‑1 (MQ=255)
                                                                                           |                                                                                                
GAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATA  >  NC_000913/3622705‑3622892

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