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,7960TG99.5% 547.3 / ‑4.4 207L1202R (CTG→CGG) rhsBRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (98/108);  total (98/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

Only 100 of 211 total aligned reads displayed.

TATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATA  >  NC_000913/3622699‑3622892
                                                                                                 |                                                                                                
 atgatgAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTa                                                                                               <  1:2122735/100‑1 (MQ=25)
 atgatgAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTa                                                                                               >  1:2701564/1‑100 (MQ=25)
 atgatgAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGAAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTa                                                                                               >  2:147432/1‑100 (MQ=255)
   gatgaATGTGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTAtt                                                                                             >  2:1764572/1‑100 (MQ=12)
   gatgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTAtt                                                                                             <  2:2482849/100‑1 (MQ=18)
    atgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTa                                                                                            >  1:2606985/1‑100 (MQ=18)
     tgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTa                                                                                            <  2:135933/99‑1 (MQ=18)
     tgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTAc                                                                                           <  1:1756312/100‑1 (MQ=18)
      gaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACa                                                                                          >  2:1425694/1‑100 (MQ=18)
          ggggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCg                                                                                      <  1:780429/100‑1 (MQ=18)
          ggggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCg                                                                                      <  2:548621/100‑1 (MQ=18)
           gggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGc                                                                                     >  2:2416800/1‑100 (MQ=17)
           gggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGc                                                                                     <  1:1815599/100‑1 (MQ=17)
            ggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGcc                                                                                    >  1:317342/1‑100 (MQ=18)
              cAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCAc                                                                                  <  1:2489367/100‑1 (MQ=18)
               aaCCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCAcc                                                                                 <  2:400357/100‑1 (MQ=18)
                aCCTGCTGAATGAAGAGAACCCGCATCAGCTGCCGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGGTGAGGAGTCCGGCCGGTAGTACAgccgcc                                                                                    >  2:3686292/1‑96 (MQ=11)
                 cttgctgAATGAAGATAGCCGGCTTGAGCTTGATCGGCTTTTCCGGCTGCCGGGGCAGCAGTATGATGATGAGTCCGGCCGGTATTACAACCTcc                                                                                    <  2:3720099/93‑1 (MQ=2)
                    gctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtat                                                                            >  2:303114/1‑100 (MQ=18)
                       gaatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtattat                                                                         <  1:844610/100‑1 (MQ=18)
                        aatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATg                                                                        >  1:3466204/1‑100 (MQ=17)
                        aatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATg                                                                        <  1:3841813/100‑1 (MQ=18)
                        aatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATg                                                                        >  1:2291078/1‑100 (MQ=18)
                         atgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGa                                                                       <  2:545981/100‑1 (MQ=18)
                           gaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAcc                                                                     >  1:437144/1‑100 (MQ=18)
                             agagAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCg                                                                   >  1:1464222/1‑100 (MQ=18)
                             agagAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCg                                                                   >  1:3099462/1‑100 (MQ=18)
                              gagaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGc                                                                  >  2:1165412/1‑100 (MQ=18)
                               agaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCt                                                                 >  2:1016747/1‑100 (MQ=17)
                                gaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTg                                                                >  1:1495578/1‑100 (MQ=18)
                                  aCCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCCGCAGTATGGTGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCa                                                              >  1:1374490/1‑100 (MQ=255)
                                  aCCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCa                                                              <  1:3125810/100‑1 (MQ=18)
                                   cccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCa                                                              <  1:3583712/99‑1 (MQ=17)
                                   cccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAg                                                             >  2:3182228/1‑100 (MQ=18)
                                   cccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAg                                                             <  2:2982445/100‑1 (MQ=18)
                                     cGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAggg                                                           >  1:105526/1‑100 (MQ=18)
                                      gCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATCATGACCCGCTGCAgggg                                                          <  1:1669950/100‑1 (MQ=255)
                                       cATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGc                                                         <  2:3371529/100‑1 (MQ=17)
                                       cATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGc                                                         <  1:985724/100‑1 (MQ=18)
                                        aTCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCg                                                        <  1:701508/100‑1 (MQ=18)
                                         tCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGa                                                       <  1:1984200/100‑1 (MQ=18)
                                          cAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGCGTCCGGCGGGTATTCCAATGGCCACCGCTATTCTGACCCGCTGCAGGGGCGa                                                       >  2:769148/1‑99 (MQ=11)
                                          cAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGat                                                      <  2:1401980/100‑1 (MQ=17)
                                          cAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGat                                                      >  1:3948724/1‑100 (MQ=18)
                                          cAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGat                                                      <  2:2624172/100‑1 (MQ=18)
                                           aGCTGGAGCAGGTTATCCGGCTGGCGGGGCAGAAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCg                                                        <  2:1302138/97‑1 (MQ=11)
                                           aGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata                                                     >  1:3639347/1‑100 (MQ=18)
                                           aGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata                                                     <  1:1518642/100‑1 (MQ=18)
                                           aGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata                                                     >  1:2016700/1‑100 (MQ=18)
                                            gCTGCAGCAGCTTATCCGGCGGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata                                                     <  2:3681530/99‑1 (MQ=255)
                                             cTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatata                                                   <  2:2376593/100‑1 (MQ=18)
                                             cTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatata                                                   <  1:2208208/100‑1 (MQ=18)
                                             cTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatata                                                   <  1:2238655/100‑1 (MQ=18)
                                              tGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatatat                                                  >  1:3779161/1‑100 (MQ=18)
                                              tGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatatat                                                  <  1:245084/100‑1 (MQ=18)
                                                cagcagCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCa                                                >  2:481675/1‑100 (MQ=18)
                                                  gcagcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACt                                              >  1:3186011/1‑100 (MQ=18)
                                                    agcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCa                                            >  2:2966097/1‑100 (MQ=18)
                                                     gcTTATCCGGCTGCCGGGGCATCATTATGATGAGGAGTCCGGCCGGTTTTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAg                                           <  1:3449446/100‑1 (MQ=11)
                                                      cTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAgg                                          >  1:3375743/1‑100 (MQ=17)
                                                        tATCCGGCTGCCGGGGCAGCAGTATGATGCGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAgg                                          >  2:2183922/1‑98 (MQ=12)
                                                        tATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGAt                                        <  1:803313/100‑1 (MQ=18)
                                                        tATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCAATTATGACCCGCTGCAGGGGCGATATATCACTCAGGan                                        <  2:1221082/100‑2 (MQ=255)
                                                          tCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATcc                                      <  2:1491204/100‑1 (MQ=17)
                                                             ggCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGAt                                   <  2:2724886/100‑1 (MQ=14)
                                                                 gCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGa                               <  1:2906399/100‑1 (MQ=11)
                                                                  ccGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGAc                              <  1:1907363/100‑1 (MQ=14)
                                                                    ggggCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTg                            >  1:4044941/1‑100 (MQ=18)
                                                                      ggCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGaa                          >  2:4044901/1‑100 (MQ=18)
                                                                        cagcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgg                        >  1:1195064/1‑100 (MQ=25)
                                                                        cagcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgg                        >  2:2878662/1‑100 (MQ=25)
                                                                         agcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAggg                       <  1:1320439/100‑1 (MQ=255)
                                                                          gcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggg                      <  2:1296974/100‑1 (MQ=255)
                                                                          gcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggg                      >  1:4095207/1‑100 (MQ=255)
                                                                           cagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAggggg                     >  1:948038/1‑100 (MQ=255)
                                                                            agTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAggggg                     <  2:1104309/99‑1 (MQ=25)
                                                                             gTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACGCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgggggg                    >  2:1728707/1‑99 (MQ=255)
                                                                              tATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGAt                  <  2:596339/100‑1 (MQ=255)
                                                                               atgatgAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATg                 <  2:3254920/100‑1 (MQ=255)
                                                                               atgatgAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATg                 >  1:3238167/1‑100 (MQ=255)
                                                                                tgatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATgg                >  1:621552/1‑100 (MQ=255)
                                                                                 gatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGa               >  2:3608113/1‑100 (MQ=255)
                                                                                 gatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGa               >  1:100459/1‑100 (MQ=255)
                                                                                  atgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGaa              >  1:1416018/1‑100 (MQ=255)
                                                                                  atgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGaa              >  2:1245638/1‑100 (MQ=255)
                                                                                   tgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAc             >  2:3480832/1‑100 (MQ=255)
                                                                                     aggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCt           >  1:1614813/1‑100 (MQ=255)
                                                                                      ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg          <  2:914967/100‑1 (MQ=255)
                                                                                      ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg          >  1:2045573/1‑100 (MQ=255)
                                                                                      ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg          <  2:3917414/100‑1 (MQ=255)
                                                                                      ggagCCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGAGGCAATATATCCCTCAGGATCAGCAAGGACTGAAGGGGGGACGGAACCTg          >  2:1610844/1‑100 (MQ=14)
                                                                                         gTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTAt       <  2:3749111/100‑1 (MQ=255)
                                                                                          tCCGGCCGGTATTACAGCCGCCACCGCTATTATGACCCGCTGCAGGGGCGATAGATCTCTCAGGTTCCGATTAGCCTGAAGGGGGGATGGAACCTGTATg      <  2:3961982/100‑1 (MQ=14)
                                                                                          tCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATg      <  1:2960303/100‑1 (MQ=255)
                                                                                          tCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATg      <  2:3237853/100‑1 (MQ=255)
                                                                                          tCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATg      >  2:180411/1‑100 (MQ=255)
                                                                                           ccggccggTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATgg     >  1:3743886/1‑100 (MQ=255)
                                                                                           ccggccggTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATgg     >  2:3522988/1‑100 (MQ=255)
                                                                                            cGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGa    <  2:4128752/100‑1 (MQ=255)
                                                                                              gCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGata  >  1:1775992/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TATGATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATA  >  NC_000913/3622699‑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: