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.4% 416.0 / ‑4.7 162L1202R (CTG→CGG) rhsBRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (74/86);  total (75/87)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

Only 100 of 165 total aligned reads displayed.

GATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATA  >  NC_000913/3622702‑3622892
                                                                                              |                                                                                                
gatgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTAtt                                                                                             >  1:3087671/1‑100 (MQ=33)
gatgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTAtt                                                                                             >  1:1572921/1‑100 (MQ=17)
gatgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTAtt                                                                                             >  2:2738066/1‑100 (MQ=18)
 atgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTa                                                                                            <  2:3149631/100‑1 (MQ=18)
 atgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTa                                                                                            <  1:1933597/100‑1 (MQ=18)
  tgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTAc                                                                                           <  2:2959418/100‑1 (MQ=18)
  tgaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTAc                                                                                           <  1:1129441/100‑1 (MQ=18)
   gaATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACa                                                                                          >  2:730226/1‑100 (MQ=18)
    aaTGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACaa                                                                                         >  2:3308854/1‑100 (MQ=18)
       ggggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCg                                                                                      <  1:1131351/100‑1 (MQ=18)
        gggCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGc                                                                                     >  1:3380411/1‑100 (MQ=18)
          gCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCa                                                                                   >  2:505911/1‑100 (MQ=18)
          gCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCa                                                                                   >  2:2391223/1‑100 (MQ=18)
            aaCCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCAcc                                                                                 <  1:2680851/100‑1 (MQ=18)
            aaCCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCAcc                                                                                 <  1:3331001/100‑1 (MQ=18)
             aCCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCAcc                                                                                 <  1:696176/99‑1 (MQ=18)
               ctgctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTAATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCAc                                                                                  >  1:2092021/1‑96 (MQ=12)
                tgctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCta                                                                             >  2:1039678/1‑100 (MQ=18)
                  ctgAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtatt                                                                           >  2:102712/1‑100 (MQ=18)
                   tgaatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtatta                                                                          <  1:1728933/100‑1 (MQ=18)
                    gaatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtcttat                                                                         >  1:930234/1‑100 (MQ=12)
                    gaatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCtattat                                                                         <  1:101164/100‑1 (MQ=18)
                     aatgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATg                                                                        <  2:2799303/100‑1 (MQ=18)
                      atgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGa                                                                       <  2:797843/100‑1 (MQ=18)
                       tgaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAc                                                                      >  2:1281911/1‑100 (MQ=18)
                        gaaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAcc                                                                     <  1:3101427/100‑1 (MQ=17)
                         aaGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGAccc                                                                    >  1:2515055/1‑100 (MQ=17)
                           gagaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGc                                                                  <  2:3252610/100‑1 (MQ=18)
                           gagaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGc                                                                  <  1:95879/100‑1 (MQ=18)
                            agaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCt                                                                 >  2:2661819/1‑100 (MQ=18)
                             gaACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTg                                                                >  2:1382357/1‑100 (MQ=18)
                              aaCCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGc                                                               <  1:2716505/100‑1 (MQ=18)
                               aCCCGCATCAGCTGCAGCAGTTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCa                                                              <  2:2185977/100‑1 (MQ=255)
                               aCCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCa                                                              <  1:2391206/100‑1 (MQ=18)
                                cccGCATCAGCTGCAGCCGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAg                                                             >  1:1279590/1‑100 (MQ=12)
                                 ccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAgg                                                            >  2:376176/1‑100 (MQ=18)
                                 ccGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAgg                                                            >  2:3293598/1‑100 (MQ=18)
                                    cATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGc                                                         <  1:1400869/100‑1 (MQ=17)
                                     aTCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCg                                                        <  1:3141994/100‑1 (MQ=17)
                                       cAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGat                                                      <  1:493759/100‑1 (MQ=18)
                                        aGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata                                                     <  2:1164276/100‑1 (MQ=17)
                                        aGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata                                                     >  1:45951/1‑100 (MQ=18)
                                        aGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata                                                     <  2:2205580/100‑1 (MQ=18)
                                        aGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata                                                     <  1:3351497/100‑1 (MQ=17)
                                         gCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGCTGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatat                                                    <  1:3224286/100‑1 (MQ=12)
                                         gCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatat                                                    >  1:2339724/1‑100 (MQ=17)
                                         gCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGata                                                     >  1:3172372/1‑99 (MQ=17)
                                          cTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatata                                                   >  1:2387824/1‑100 (MQ=18)
                                          cTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGGATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGatata                                                   <  1:206399/100‑1 (MQ=12)
                                            gcagcaGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATc                                                 >  1:2156569/1‑100 (MQ=18)
                                            gcagcaGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATc                                                 <  2:2549188/100‑1 (MQ=17)
                                             cagcagCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCa                                                >  2:2181622/1‑100 (MQ=18)
                                              agcagcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCAc                                               <  2:820487/100‑1 (MQ=18)
                                               gcagcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACt                                              >  1:457532/1‑100 (MQ=18)
                                                cagcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTc                                             <  2:1955168/100‑1 (MQ=18)
                                                 agcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCa                                            <  2:1679749/100‑1 (MQ=18)
                                                  gcTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAg                                           <  2:285871/100‑1 (MQ=18)
                                                   cTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAgg                                          >  1:747215/1‑100 (MQ=18)
                                                      aTCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATc                                       <  2:1284346/100‑1 (MQ=18)
                                                       tCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACAGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATcc                                      >  1:66274/1‑100 (MQ=12)
                                                         cGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGAt                                        <  1:2163475/96‑1 (MQ=17)
                                                          ggCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGTCCCGCTGCAGGGGCGATATATCACTCAGGATCCGAt                                   <  2:2387841/100‑1 (MQ=9)
                                                          ggCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGAt                                   >  2:3222388/1‑100 (MQ=14)
                                                           gCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGAtt                                  <  1:1780155/100‑1 (MQ=11)
                                                              gCCGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGa                               <  1:1415140/100‑1 (MQ=11)
                                                               ccGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGAc                              >  2:2444275/1‑100 (MQ=14)
                                                               ccGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGAc                              <  2:747218/100‑1 (MQ=14)
                                                                cGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGAc                              >  2:1959929/1‑99 (MQ=14)
                                                                cGGGGCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGACACCGCTATTATGACCCGCTGCAGGGGCGCTATATCACTCAGGATCCGATTGGAc                              >  2:1262831/1‑99 (MQ=11)
                                                                 ggggCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTg                            <  2:1997709/100‑1 (MQ=17)
                                                                 ggggCAGCAGTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTg                            <  2:2142314/100‑1 (MQ=17)
                                                                     cagcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGGAgg                        >  1:2533231/1‑100 (MQ=255)
                                                                     cagcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgg                        <  1:3010335/100‑1 (MQ=25)
                                                                     cagcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgg                        >  2:1552224/1‑100 (MQ=25)
                                                                      agcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGCAggg                       >  1:1985979/1‑100 (MQ=255)
                                                                      agcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAgg                        >  1:2704816/1‑99 (MQ=25)
                                                                       gcagTATGATGAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATACATCACTCAGCATCCGATTGGACTGAAggag                      >  1:3241854/1‑98 (MQ=11)
                                                                            atgatgAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGAt                  <  1:1521479/99‑1 (MQ=25)
                                                                            atgatgAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATg                 <  2:2270977/100‑1 (MQ=255)
                                                                            atgatgAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATg                 <  1:3240135/100‑1 (MQ=255)
                                                                            atgatgAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATg                 >  2:3352646/1‑100 (MQ=255)
                                                                            atgatgAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATg                 <  1:15413/100‑1 (MQ=25)
                                                                            atgatgAGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATg                 <  2:2042999/100‑1 (MQ=255)
                                                                             tgatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATgg                >  1:434575/1‑100 (MQ=255)
                                                                              gatgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGa               >  2:2575465/1‑100 (MQ=25)
                                                                               atgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGaa              >  2:972936/1‑100 (MQ=25)
                                                                               atgaGGAGTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGa               >  2:1663836/1‑99 (MQ=255)
                                                                                 gaggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAcc            >  2:1524862/1‑100 (MQ=255)
                                                                                  aggaggCCGGCTCGTTTGACAGCCGTCACTGGTTTTATGGCGCGGTGCAGGGGCGCTGTATCACTCAGGATCTGTTTGGGCTGTAGGGGGGATGGAACCt           <  1:2076318/100‑1 (MQ=255)
                                                                                   ggtgTCCGGCCGGTATTACATCCTCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAAcc            <  1:2821208/95‑1 (MQ=14)
                                                                                   ggagTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTg          >  1:879040/1‑100 (MQ=255)
                                                                                     agTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTa        <  2:2349540/100‑1 (MQ=255)
                                                                                     agTCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTa        <  2:2054737/100‑1 (MQ=255)
                                                                                       tCCGGCCGGTATTACAATTGCCACTGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATg      <  2:846525/100‑1 (MQ=14)
                                                                                       tCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATg      <  2:29994/100‑1 (MQ=255)
                                                                                       tCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATg      <  2:385107/100‑1 (MQ=255)
                                                                                       tCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATg      <  1:372225/100‑1 (MQ=255)
                                                                                       tCCGGCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATg      <  1:1039675/100‑1 (MQ=255)
                                                                                        ccGGACGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATg      <  2:347853/99‑1 (MQ=255)
                                                                                           gCCGGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGata  <  1:2925456/100‑1 (MQ=255)
                                                                                              |                                                                                                
GATGAATGGGGCAACCTGCTGAATGAAGAGAACCCGCATCAGCTGCAGCAGCTTATCCGGCTGCCGGGGCAGCAGTATGATGAGGAGTCCGGCCTGTATTACAACCGCCACCGCTATTATGACCCGCTGCAGGGGCGATATATCACTCAGGATCCGATTGGACTGAAGGGGGGATGGAACCTGTATGGATA  >  NC_000913/3622702‑3622892

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