Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,887,352 T→A 100% intergenic (‑133/+226) ygcB ← / ← cysH Cascade complex anti‑viral R‑loop helicase‑annealase Cas3/phosphoadenosine phosphosulfate reductase; PAPS reductase, thioredoxin dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,887,3520TA98.7% 397.0 / ‑4.5 121intergenic (‑133/+226)ygcB/cysHCascade complex anti‑viral R‑loop helicase‑annealase Cas3/phosphoadenosine phosphosulfate reductase; PAPS reductase, thioredoxin dependent
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (57/62);  total (58/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

ATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCA  >  NC_000913/2887256‑2887448
                                                                                                |                                                                                                
 tCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGttt                                                                                              >  1:575506/1‑100 (MQ=255)
  cAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTg                                                                                             <  1:1745850/100‑1 (MQ=255)
  cAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTg                                                                                             <  1:1495222/100‑1 (MQ=255)
  cAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTg                                                                                             <  1:2333029/100‑1 (MQ=255)
   aGCGGGCTATTTTTTTATTTTTTTTATACATATTTTATTACAGTGATATACATATTTCTTTTATCCCTGTCAATTAGAGCCTAACATTTGAATAGttttg                                                                                            >  2:1767964/1‑98 (MQ=255)
   aGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACCTTCGAATAGTTTgg                                                                                            >  2:1132281/1‑100 (MQ=255)
   aGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTgg                                                                                            >  1:988076/1‑100 (MQ=255)
   aGCGGGCTATTTTTCGATGCTGCTCATAAATATCCAATTCCAGTGATATACCCATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTg                                                                                             >  1:1211523/1‑99 (MQ=255)
      gggCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGtt                                                                                         >  1:227020/1‑100 (MQ=255)
       ggCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTc                                                                                        <  1:1327552/100‑1 (MQ=255)
        gCTATTTTTCGATGCTGCTCATAAATATCGTATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTc                                                                                        <  2:1396958/99‑1 (MQ=255)
         cTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCg                                                                                       <  1:664300/99‑1 (MQ=255)
         cTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGa                                                                                      >  1:1938416/1‑100 (MQ=255)
         cTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGa                                                                                      <  2:609344/100‑1 (MQ=255)
         cTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGa                                                                                      >  1:453499/1‑100 (MQ=255)
         cTAATTTTCGATGTTGCTCATAAATAACGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGATACTAATATTCGAATAGTTTGGGTTCGa                                                                                      <  2:1181175/100‑1 (MQ=255)
              tttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCt                                                                                 <  2:411463/100‑1 (MQ=255)
                tCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTgg                                                                               <  1:2213837/100‑1 (MQ=255)
                 cGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTgg                                                                               >  1:1953478/1‑99 (MQ=255)
                 cGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGc                                                                              >  1:2223417/1‑100 (MQ=255)
                  gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                             <  1:1504372/100‑1 (MQ=255)
                  gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                             <  2:601350/100‑1 (MQ=255)
                  gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                             >  1:194991/1‑100 (MQ=255)
                  gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGc                                                                              <  2:954877/99‑1 (MQ=255)
                    tgctgcTCATAAATATCTAATTACAGTGTTATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTc                                                                           <  2:1065811/100‑1 (MQ=255)
                     gcttctCATAAATATCGAATTACAGTGATAGACACATGTCGTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCa                                                                          >  1:1039660/1‑100 (MQ=255)
                     gctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCa                                                                          >  2:878009/1‑100 (MQ=255)
                      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAg                                                                         <  1:2011308/100‑1 (MQ=255)
                      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAg                                                                         >  2:429568/1‑100 (MQ=255)
                        gcttctAAATATTGAATTAAATTGATATATACATGTTTTGCCTATTTGTAAAGAAGAGACTAACATTAGAATAGTTTGGGCTCGAAGGTTGGCCTCa                                                                          >  2:819823/6‑97 (MQ=255)
                        gctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTAGCCTCa                                                                          >  2:2259293/1‑97 (MQ=255)
                           cATAAATATCGAATTACAGTGATATACACGTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGATGGCCTCAGGTTGa                                                                    >  1:37208/1‑100 (MQ=255)
                           cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGa                                                                    <  1:1542623/100‑1 (MQ=255)
                           cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGa                                                                    >  1:1579657/1‑100 (MQ=255)
                           cATAAATATCAAATTACAGATATATACACATGTCTTACCTCCCTATAAAGAAGAGACTAACATTCGAATAGTTAGGATTCGAACGCAGGCCTCAGGTTGa                                                                    >  2:2090604/1‑100 (MQ=255)
                            aTAAATTTCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGAt                                                                   <  1:1915172/100‑1 (MQ=255)
                             tAAATATAGAATTTCATTGTTATACACATGTTTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGTTTCGAATGCTGGCCTCAGGTTGAt                                                                   <  2:287012/99‑1 (MQ=255)
                                 tatCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGaan                                                              <  1:1996219/100‑2 (MQ=255)
                                  atCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAat                                                             >  1:1203457/1‑100 (MQ=255)
                                   tCGAATTACAGTGATATACACATGTCTTGCCTCCCTTCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGaaa                                                              >  2:1539024/1‑98 (MQ=255)
                                   tCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAata                                                            <  2:960528/100‑1 (MQ=255)
                                     gAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATc                                                          >  1:2039422/1‑100 (MQ=255)
                                      aaTTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCg                                                         <  1:959468/100‑1 (MQ=255)
                                      aaTTACAGTGATATACACATGTCTTGCCTCCCTGAAAAGAAGAGACTAACATTCGAATAGGTTGGGTTCGAACGCTAGACTAAGGTTAATAGAAatat                                                           >  2:2148288/1‑98 (MQ=255)
                                        ttACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGcc                                                       <  2:300123/100‑1 (MQ=255)
                                        ttACAGTGATATAAACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGcc                                                       >  1:1624874/1‑100 (MQ=255)
                                          aCAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTg                                                     <  1:911426/100‑1 (MQ=255)
                                            aGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTggg                                                   <  1:1252650/100‑1 (MQ=255)
                                             gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgggg                                                  >  1:1279922/1‑100 (MQ=255)
                                              tGATATACACATGTCTTGCCTCCCCGCAAAGAAGAGACTAACATTCGAATAGTGTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgggg                                                  <  2:38192/99‑1 (MQ=255)
                                               gATATCCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgg                                                    >  2:1289139/1‑96 (MQ=255)
                                               gATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCt                                                <  1:1061592/100‑1 (MQ=255)
                                                 tataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCttt                                              >  2:974524/1‑100 (MQ=255)
                                                 tataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCttt                                              >  1:1715858/1‑100 (MQ=255)
                                                  ataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCtttt                                             <  2:1949811/100‑1 (MQ=255)
                                                   taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                            >  1:2338518/1‑100 (MQ=255)
                                                   taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                            >  1:165830/1‑100 (MQ=255)
                                                   taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGACTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                            <  2:944608/100‑1 (MQ=255)
                                                   taCACATGTCTTGCCTCCCTGAAAAAAAGAAACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGATACTTTTg                                            >  2:1017788/1‑100 (MQ=255)
                                                    acacATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGt                                           <  1:1840539/100‑1 (MQ=255)
                                                    acacATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGt                                           <  2:1958237/100‑1 (MQ=255)
                                                     cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                          <  1:818603/100‑1 (MQ=255)
                                                     cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                          <  1:206148/100‑1 (MQ=255)
                                                        aTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCAt                                       >  1:1494610/1‑100 (MQ=255)
                                                          gTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGATATATCGCATGGGGCTTTTGTCCAt                                       >  2:236412/1‑98 (MQ=255)
                                                           tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg                                    >  1:1918127/1‑100 (MQ=255)
                                                           tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg                                    <  2:1287882/100‑1 (MQ=255)
                                                           tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg                                    <  2:2243125/100‑1 (MQ=255)
                                                                ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                               <  1:150592/100‑1 (MQ=255)
                                                                ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                               >  1:959874/1‑100 (MQ=255)
                                                                ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                               >  2:128632/1‑100 (MQ=255)
                                                                    ccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATATAAATATCAACTGGGGCTTTTGTCCATCTGGAACCTcg                            >  1:2021439/1‑99 (MQ=255)
                                                                    ccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc                           <  2:142920/100‑1 (MQ=255)
                                                                     cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgcg                          <  1:411198/100‑1 (MQ=255)
                                                                     cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgcg                          <  1:98784/100‑1 (MQ=255)
                                                                      tGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGa                         >  1:1839666/1‑100 (MQ=255)
                                                                       gCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGaa                        >  1:1051105/1‑100 (MQ=255)
                                                                       gCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGaa                        <  1:457311/100‑1 (MQ=255)
                                                                        cAAAGAAGCGACTAACATTCGAATAGTTTGGGTTCTGTCGCTGGCCTTAGGTTGATAGAAATATCGCCTGGGGTTTTTGTCCATCTGGAACCTcg                            <  2:1018930/95‑1 (MQ=255)
                                                                        cAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAAt                       >  2:1938814/1‑100 (MQ=255)
                                                                         aaaGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATg                      >  2:891293/1‑100 (MQ=255)
                                                                         aaaGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATg                      <  2:2213414/100‑1 (MQ=255)
                                                                          aagaagAGGCTAACATTCGAATAGTTTAGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGc                     <  1:888598/100‑1 (MQ=255)
                                                                          aagaagAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGc                     <  1:1373466/100‑1 (MQ=255)
                                                                          aagaagAGACTAACATTCGAATAGATTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGaa                        >  2:2054829/1‑97 (MQ=255)
                                                                           agaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCt                    <  1:1388381/100‑1 (MQ=255)
                                                                               gagaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAAc                <  2:1039650/100‑1 (MQ=255)
                                                                                agaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACg               <  2:2039400/100‑1 (MQ=255)
                                                                                 gaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGc              >  1:1942651/1‑100 (MQ=255)
                                                                                 gaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGc              >  2:1522835/1‑100 (MQ=255)
                                                                                  aCTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc             <  2:1953456/100‑1 (MQ=255)
                                                                                   cTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCa            <  2:194988/100‑1 (MQ=255)
                                                                                     aaCATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCa            <  1:543158/98‑1 (MQ=255)
                                                                                      aCATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGAc         <  1:596471/100‑1 (MQ=255)
                                                                                      aCATTCGAATAGTTTGGGTTCGAACGCTAACCTCAGGTTGATAGAAATATCGACTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGAc         >  2:1047124/1‑100 (MQ=255)
                                                                                       cATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACa        >  1:780595/1‑100 (MQ=255)
                                                                                        aTTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAg       >  2:1629620/1‑100 (MQ=255)
                                                                                         ttCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCGTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGc      >  2:2248913/1‑100 (MQ=255)
                                                                                           cGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCGTCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCt    <  1:1692659/100‑1 (MQ=255)
                                                                                             aaTAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCa  <  2:1461019/100‑1 (MQ=255)
                                                                                                |                                                                                                
ATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCA  >  NC_000913/2887256‑2887448

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: