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,3520TA99.8% 568.0 / ‑6.6 172intergenic (‑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 (0/1);  major base A (92/78);  minor base G (0/1);  total (92/80)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 172 total aligned reads displayed.

ATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCA  >  NC_000913/2887256‑2887448
                                                                                                |                                                                                                
  cAGCGGGCTATTTTTCGATTCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTg                                                                                             >  2:1877233/1‑100 (MQ=255)
  cAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTg                                                                                             <  1:3518377/100‑1 (MQ=255)
   aGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTgg                                                                                            <  1:3124511/100‑1 (MQ=255)
   aGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTgg                                                                                            <  2:873141/100‑1 (MQ=255)
   aGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTgg                                                                                            >  1:1006355/1‑100 (MQ=255)
   aGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTgg                                                                                            <  2:406057/100‑1 (MQ=255)
    gCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTggg                                                                                           <  1:143544/100‑1 (MQ=255)
    gCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTggg                                                                                           <  1:3387188/100‑1 (MQ=255)
     cGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGt                                                                                          <  1:797478/100‑1 (MQ=255)
     cGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGt                                                                                          <  1:412739/100‑1 (MQ=255)
      gggCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGtt                                                                                         >  1:1298597/1‑100 (MQ=255)
      gggCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCAAATAGTTTGGGtt                                                                                         <  1:875136/100‑1 (MQ=255)
       ggCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTc                                                                                        <  2:1731082/100‑1 (MQ=255)
       ggCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTc                                                                                        <  2:633446/100‑1 (MQ=255)
        gCTATTTTTCGATGCTGCTCATAAATATCGAATTACATTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCg                                                                                       <  2:989554/100‑1 (MQ=255)
          tATTTTTCGATGCTGCTCATATATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCTAAGAGGAGACTAACATTCGGATAGTTTGGGTTCGa                                                                                      >  2:3870231/1‑99 (MQ=255)
           aTTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACCCATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAAc                                                                                    <  2:1981502/100‑1 (MQ=255)
           aTTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAAc                                                                                    >  1:2191069/1‑100 (MQ=255)
           aTTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAAc                                                                                    <  1:2274731/100‑1 (MQ=255)
           aTTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAAc                                                                                    >  2:1188830/1‑100 (MQ=255)
            tttttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGaa                                                                                     <  2:252119/98‑1 (MQ=255)
              tttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCt                                                                                 >  2:3195682/1‑100 (MQ=255)
              tttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCt                                                                                 >  2:436035/1‑100 (MQ=255)
              tttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCt                                                                                 >  1:2428556/1‑100 (MQ=255)
               ttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTg                                                                                <  1:1832200/100‑1 (MQ=255)
                tCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTgg                                                                               <  2:440127/100‑1 (MQ=255)
                 cGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGc                                                                              >  2:1386493/1‑100 (MQ=255)
                   aTGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCt                                                                            >  2:1049191/1‑100 (MQ=255)
                    tgctgcTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTc                                                                           >  1:152278/1‑100 (MQ=255)
                     gctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCa                                                                          <  1:2630709/100‑1 (MQ=255)
                       tgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAgg                                                                        <  2:2747281/100‑1 (MQ=255)
                       tgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAgg                                                                        <  1:1593536/100‑1 (MQ=255)
                        gctCATAAATATCTAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACGTTCGAATAGTTTGGGTTCGAACGCTGGCCTCAg                                                                         <  2:331987/98‑1 (MQ=255)
                         ctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGtt                                                                      >  2:3061573/1‑100 (MQ=255)
                             tAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATa                                                                  >  1:3109069/1‑100 (MQ=255)
                             tAAATATCGAACTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATGGTTTGGGTTCGAACGCTGGCCTCAGGTTGAt                                                                   <  2:3977792/99‑1 (MQ=255)
                              aaaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAg                                                                 <  1:1141121/100‑1 (MQ=255)
                               aaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGa                                                                >  1:299474/1‑100 (MQ=255)
                                      aaTTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGATTAACATTCGCATAGTTTGCGTTCGAACGCTGGGCCCAGGTTGATAGAACTATc                                                          >  2:139219/1‑99 (MQ=255)
                                      aaTTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCg                                                         <  2:1448460/100‑1 (MQ=255)
                                       aTTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTAGCCTCAGGTTGATAGAAATATCGc                                                        <  2:4111609/100‑1 (MQ=255)
                                         tACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCt                                                      <  1:2745159/100‑1 (MQ=255)
                                            aGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACGAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTggg                                                   <  2:2589340/100‑1 (MQ=255)
                                             gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgggg                                                  >  2:1950563/1‑100 (MQ=255)
                                             gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgggg                                                  >  2:207350/1‑100 (MQ=255)
                                             gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgggg                                                  >  1:438936/1‑100 (MQ=255)
                                                atatACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCtt                                               >  2:857143/1‑100 (MQ=255)
                                                 tataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCttt                                              >  2:1777954/1‑100 (MQ=255)
                                                  ataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCtttt                                             <  1:1791856/100‑1 (MQ=255)
                                                   taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                            <  1:1290952/100‑1 (MQ=255)
                                                   taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                            >  2:22824/1‑100 (MQ=255)
                                                   taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                            <  1:1703827/100‑1 (MQ=255)
                                                   taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                            >  1:3768025/1‑100 (MQ=255)
                                                     cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                          >  1:847416/1‑100 (MQ=255)
                                                     cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                          <  2:1108624/100‑1 (MQ=255)
                                                     cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                          <  1:3979602/100‑1 (MQ=255)
                                                       caTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCa                                        >  1:258812/1‑100 (MQ=255)
                                                       caTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCa                                        >  2:1361682/1‑100 (MQ=255)
                                                       caTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCa                                        <  1:1950576/100‑1 (MQ=255)
                                                         tGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATc                                      >  1:3554137/1‑100 (MQ=255)
                                                           tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg                                    <  1:2185445/100‑1 (MQ=255)
                                                           tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg                                    >  1:2971644/1‑100 (MQ=255)
                                                                ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                               >  2:3734411/1‑100 (MQ=255)
                                                                ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                               >  2:3378804/1‑100 (MQ=255)
                                                                 cTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCt                              >  1:2705536/1‑100 (MQ=255)
                                                                 cTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCt                              >  1:2435980/1‑100 (MQ=255)
                                                                   cccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcg                            <  2:1317568/100‑1 (MQ=255)
                                                                    ccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc                           >  1:1963832/1‑100 (MQ=255)
                                                                     cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgcg                          >  1:1871041/1‑100 (MQ=255)
                                                                     cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgcg                          <  1:2110510/100‑1 (MQ=255)
                                                                     cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcacg                          >  1:1582178/1‑97 (MQ=255)
                                                                        cAAAGTAGAGACTACTATTCGAATAGCGTGGGTTAGAACGCTGGCCTCAGGTTGATAGAAAAATCGCCTGGGGCTTTTGTCCATCTGGGACCAcgcg                          <  2:3387161/97‑1 (MQ=255)
                                                                        cAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAAt                       <  2:1006349/100‑1 (MQ=255)
                                                                        cAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCATGGGGCTTTTGTCCATCGGGaa                                 >  2:2529671/1‑90 (MQ=255)
                                                                        cAAAGAAGAGACTAAAATTCGGATAGTTTGGGTTTCAACGCTTGCCGCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGaa                        <  2:2277260/99‑1 (MQ=255)
                                                                         aaaGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATg                      <  1:2367248/100‑1 (MQ=255)
                                                                          aagaagAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGc                     <  1:2230467/100‑1 (MQ=255)
                                                                           agaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCt                    >  2:2918127/1‑100 (MQ=255)
                                                                           agaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCt                    <  1:1815632/100‑1 (MQ=255)
                                                                            gaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCtt                   >  1:2537311/1‑100 (MQ=255)
                                                                             aagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTa                  <  2:3554114/100‑1 (MQ=255)
                                                                               gagaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAAc                >  2:697762/1‑100 (MQ=255)
                                                                               gagaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAAc                <  1:3255301/100‑1 (MQ=255)
                                                                                 gaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGc              >  2:590464/1‑100 (MQ=255)
                                                                                 gaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGc              <  1:136718/100‑1 (MQ=255)
                                                                                   cTAACATTCGAATAGTTTGTGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTCGAACCTCGCGATTGCTTAAc                >  2:1236692/1‑96 (MQ=255)
                                                                                   cTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCa            <  1:1331255/100‑1 (MQ=255)
                                                                                   cTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCa            >  1:1395785/1‑100 (MQ=255)
                                                                                     aaCATTCGAATATTATGGCATCGAACGCTGCCCACAGGTTTATAGAAAGAGCGCCCTGTGCTTTTGTCAATCTGGAAGCTCGCGAATGCTTAACGCCAg           <  2:1548506/99‑1 (MQ=255)
                                                                                       cATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACa        >  1:3261925/1‑100 (MQ=255)
                                                                                       cATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACa        >  1:611832/1‑100 (MQ=255)
                                                                                       cATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACa        <  2:3103058/100‑1 (MQ=255)
                                                                                       cATTCGAATAGTTTGGGTTCGAACGCCGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACa        >  1:2520960/1‑100 (MQ=255)
                                                                                         ttCGAATAGTTTGGGTTCGAACGCTGTACTCAGGTTGATAGAAATTTCGCCTGGGTCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGc      <  2:4070897/100‑1 (MQ=255)
                                                                                          tCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGTTTAACGCCAGACAGcc     <  1:1121165/100‑1 (MQ=255)
                                                                                          tCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTCACGCCAGACAgac     <  1:2805534/100‑3 (MQ=255)
                                                                                          tCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGc      <  1:2062863/99‑1 (MQ=255)
                                                                                           cGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCt    <  2:3719033/100‑1 (MQ=255)
                                                                                            gAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTc   >  2:3088879/1‑100 (MQ=255)
                                                                                             aaTAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCa  >  2:1338556/1‑100 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: