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,3520TA97.0% 325.8 / ‑0.9 98intergenic (‑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 (2/1);  new base A (52/43);  total (54/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTC  >  NC_000913/2887255‑2887447
                                                                                                 |                                                                                               
tatCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGt                                                                                               >  1:1057845/1‑100 (MQ=255)
 atCAGCGGGCTATGTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGtt                                                                                              >  2:1586461/1‑100 (MQ=255)
   cAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTg                                                                                            >  1:715204/1‑100 (MQ=255)
    aGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTgg                                                                                           >  1:617152/1‑100 (MQ=255)
     gCGGGCTATTTTGCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTggg                                                                                          <  2:2279864/100‑1 (MQ=255)
       gggCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGtt                                                                                        <  1:1579306/100‑1 (MQ=255)
       gggCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGtt                                                                                        <  1:1244513/100‑1 (MQ=255)
         gCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCg                                                                                      <  1:486258/100‑1 (MQ=255)
          cTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGa                                                                                     <  2:1854065/100‑1 (MQ=255)
          cTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGa                                                                                     >  1:1605693/1‑100 (MQ=255)
           tATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGaa                                                                                    <  1:1816532/100‑1 (MQ=255)
              ttttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGc                                                                                 <  2:560143/100‑1 (MQ=255)
                ttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTg                                                                               <  1:872630/100‑1 (MQ=255)
                ttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTg                                                                               >  1:183289/1‑100 (MQ=255)
                  cGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTgg                                                                              <  1:221767/99‑1 (MQ=255)
                   gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                            <  2:1603386/100‑1 (MQ=255)
                   gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                            <  2:1868816/100‑1 (MQ=255)
                     tgctgcTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTc                                                                          >  1:2427730/1‑100 (MQ=255)
                      gctgctCATAAATATCGAATTACAGTGATCTACACATGTTTTGCCTCCCTGCAAAGAAGAGACTAAGATTCGAATAGTTTGGGTTCGAACGCTGGCCTCa                                                                         <  2:1543905/100‑1 (MQ=255)
                      gctgctCAGAAATATCGAAGTACAGTGATATACACATGTCTTGCCGCCTTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCa                                                                         <  2:830082/100‑1 (MQ=255)
                       ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAg                                                                        >  1:2289579/1‑100 (MQ=255)
                       ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAg                                                                        <  1:879938/100‑1 (MQ=255)
                            cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTAGCCTCAGGTTg                                                                    >  2:2127756/1‑99 (MQ=255)
                                aaTATCGAATTACAGTGATATACACATGTCTTGGCTGCGTGGAAAGAAGAGACTAACATGCGAATAGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGa                                                               <  2:750995/100‑1 (MQ=255)
                                  tatCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGaaa                                                             >  1:2026559/1‑100 (MQ=255)
                                    tCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAata                                                           <  1:1584086/100‑1 (MQ=255)
                                    tCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAata                                                           <  1:1679899/100‑1 (MQ=255)
                                     cGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAatat                                                          >  1:2263877/1‑100 (MQ=255)
                                      gAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATc                                                         >  2:378844/1‑100 (MQ=255)
                                      gAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATc                                                         <  2:2289565/100‑1 (MQ=255)
                                            cAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgg                                                   <  2:508247/100‑1 (MQ=255)
                                            cAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgg                                                   >  2:1470680/1‑100 (MQ=255)
                                              gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgggg                                                 >  2:742967/1‑100 (MQ=255)
                                              gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgggg                                                 >  2:1335154/1‑100 (MQ=255)
                                                gATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCt                                               <  2:16868/100‑1 (MQ=255)
                                                  tataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCttt                                             <  2:525354/100‑1 (MQ=255)
                                                   ataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCtttt                                            <  2:647713/100‑1 (MQ=255)
                                                    tgcacaTGTCTTGCCTCCCTGCAAAGAAGAGACTACCATTCCAATAGTTTGGGTTCCAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                           <  1:1793957/98‑1 (MQ=38)
                                                    taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                           >  2:2127185/1‑100 (MQ=255)
                                                    taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                           >  2:1091288/1‑100 (MQ=255)
                                                    taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                           >  1:2212233/1‑100 (MQ=255)
                                                    taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                           >  2:645762/1‑100 (MQ=255)
                                                    taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                           >  1:100108/1‑100 (MQ=255)
                                                     acacATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGt                                          <  2:1097948/100‑1 (MQ=255)
                                                     acacATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGt                                          >  1:1031453/1‑100 (MQ=255)
                                                      nacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                         >  2:1996722/2‑100 (MQ=255)
                                                      cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                         >  2:692939/1‑100 (MQ=255)
                                                      cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                         >  1:910574/1‑100 (MQ=255)
                                                      cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                         <  1:239329/100‑1 (MQ=255)
                                                        caTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCa                                       >  1:900074/1‑100 (MQ=255)
                                                           gTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCt                                    <  1:268893/100‑1 (MQ=255)
                                                            tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg                                   <  1:2268663/100‑1 (MQ=255)
                                                            tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg                                   <  1:2444643/100‑1 (MQ=255)
                                                            tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg                                   >  2:2499378/1‑100 (MQ=255)
                                                              ttGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGa                                 <  2:128181/100‑1 (MQ=255)
                                                               tGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGaa                                >  1:549612/1‑100 (MQ=255)
                                                                 ccTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGCCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCGATCTGGCAcc                              >  2:1220327/1‑100 (MQ=38)
                                                                 ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTTGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                              <  2:1916354/100‑1 (MQ=255)
                                                                 ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                              >  1:1823381/1‑100 (MQ=255)
                                                                    cccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcg                           >  2:1130348/1‑100 (MQ=255)
                                                                    cccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcg                           >  1:1671849/1‑100 (MQ=255)
                                                                    cccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcg                           >  1:1187033/1‑100 (MQ=255)
                                                                     ccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc                          <  2:1472398/100‑1 (MQ=255)
                                                                      cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgcg                         >  2:1055582/1‑100 (MQ=255)
                                                                      cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgcg                         <  1:630514/100‑1 (MQ=255)
                                                                       tGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGa                        >  2:2361027/1‑100 (MQ=255)
                                                                        gCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGaa                       >  2:1208321/1‑100 (MQ=255)
                                                                        gCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGaa                       >  2:333796/1‑100 (MQ=255)
                                                                         cAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGGTTTTGTCCATTTGGAACCTcgcg                         >  2:2063251/1‑97 (MQ=255)
                                                                         cAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAAt                      >  1:1307913/1‑100 (MQ=255)
                                                                         cAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAAAGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTACATCTAGAACCTCGCGAAt                      >  2:986205/1‑100 (MQ=255)
                                                                          aaaGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATg                     <  2:2018647/100‑1 (MQ=255)
                                                                           acgaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGc                    <  2:1557097/98‑1 (MQ=255)
                                                                           aagaagAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGc                    >  2:178229/1‑100 (MQ=255)
                                                                           aagaagAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGc                    <  2:1441890/100‑1 (MQ=255)
                                                                            agaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCt                   <  1:118568/100‑1 (MQ=255)
                                                                            agaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCt                   >  1:2025455/1‑100 (MQ=255)
                                                                             gaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCtt                  >  1:686355/1‑100 (MQ=255)
                                                                                gagaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAAc               <  1:217927/100‑1 (MQ=255)
                                                                                 agaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACg              <  1:375076/100‑1 (MQ=255)
                                                                                    cTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc            >  2:544405/1‑99 (MQ=255)
                                                                                    cTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCa           >  1:1280508/1‑100 (MQ=255)
                                                                                    cTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCa           <  1:409443/100‑1 (MQ=255)
                                                                                      aaCATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAg          >  2:1249908/1‑99 (MQ=255)
                                                                                      aaCATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGa         >  2:229605/1‑100 (MQ=255)
                                                                                       tcATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGAc        <  1:198302/99‑1 (MQ=255)
                                                                                       aCATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGAc        >  1:2127564/1‑100 (MQ=255)
                                                                                       aCATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGAc        <  1:1070678/100‑1 (MQ=255)
                                                                                        cATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACa       >  2:198300/1‑100 (MQ=255)
                                                                                        cATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGAGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACa       >  2:1637006/1‑100 (MQ=255)
                                                                                          ttCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGCCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGc     >  1:739589/1‑100 (MQ=255)
                                                                                           tCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGcc    <  1:1483264/100‑1 (MQ=255)
                                                                                            cGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCt   >  1:1589042/1‑100 (MQ=255)
                                                                                            cGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCt   >  1:1203526/1‑100 (MQ=255)
                                                                                             gAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTc  >  1:1451165/1‑100 (MQ=255)
                                                                                             gAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTc  <  1:122247/100‑1 (MQ=255)
                                                                                             gAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTc  >  1:624614/1‑100 (MQ=255)
                                                                                             gAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTc  <  1:116344/100‑1 (MQ=255)
                                                                                                 |                                                                                               
TATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTC  >  NC_000913/2887255‑2887447

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: