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.2% 560.7 / 0.1 169intergenic (‑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/2);  new base A (86/80);  total (87/82)
Fisher's exact test for biased strand distribution p-value = 6.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.36e-01

Only 100 of 169 total aligned reads displayed.

ATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAA  >  NC_000913/2887256‑2887449
                                                                                                |                                                                                                 
  cAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTg                                                                                              <  1:3668243/100‑1 (MQ=255)
  cAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTg                                                                                              <  1:905976/100‑1 (MQ=255)
    gCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTggg                                                                                            <  1:1406704/100‑1 (MQ=255)
       ggCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTc                                                                                         <  1:1778831/100‑1 (MQ=255)
        gCTATTTTTCGATGCTGCTCATGAATAACGAAATACAGTGTTATACTCATGTCTTGCATCCCAGCAAAGAAGAGCCTAACATGCGAATAGTTGGGGGTCg                                                                                        <  2:3589720/100‑1 (MQ=255)
        gCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCg                                                                                        <  2:2420152/100‑1 (MQ=255)
         cTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGa                                                                                       >  1:33726/1‑100 (MQ=255)
         cTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGa                                                                                       <  2:582057/100‑1 (MQ=255)
           aTTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAAc                                                                                     >  1:1930087/1‑100 (MQ=255)
            tttttCTATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACg                                                                                    <  2:424149/100‑1 (MQ=255)
            tttttCGATACTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGaa                                                                                      <  2:518302/98‑1 (MQ=255)
                tCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTgg                                                                                <  2:3295050/100‑1 (MQ=255)
                tCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTgg                                                                                <  1:472577/100‑1 (MQ=255)
                 cGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGc                                                                               >  1:2621183/1‑100 (MQ=255)
                  gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                              >  1:1761429/1‑100 (MQ=255)
                  gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                              <  2:1094510/100‑1 (MQ=255)
                    tgctgcTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCtn                                                                            <  2:3632347/100‑2 (MQ=255)
                      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAg                                                                          <  2:3520922/100‑1 (MQ=255)
                       tgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAgg                                                                         >  1:3183647/1‑100 (MQ=255)
                        gctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGt                                                                        <  2:878201/100‑1 (MQ=255)
                         ctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGtt                                                                       <  2:733160/100‑1 (MQ=255)
                          tCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTg                                                                      >  2:3708570/1‑100 (MQ=255)
                           cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGa                                                                     >  2:608116/1‑100 (MQ=255)
                             tAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATa                                                                   >  1:1113904/1‑100 (MQ=255)
                              aaaTATCGAATTACAGTGATATACACATGTCTTTCCTGCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAg                                                                  <  2:2875215/100‑1 (MQ=255)
                              aaaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAg                                                                  <  1:4093709/100‑1 (MQ=255)
                               aaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGa                                                                 >  2:1335930/1‑100 (MQ=255)
                                atatCGAATTACAGTGATATACACATGTCGTGCCTCCCGGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGaa                                                                <  2:307661/100‑1 (MQ=255)
                                  atCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAat                                                              >  1:2158699/1‑100 (MQ=255)
                                   tCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAata                                                             >  1:3176496/1‑100 (MQ=255)
                                    ngAATTACAGTGATATACACATGTCTAGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAatat                                                            >  1:2541296/2‑100 (MQ=255)
                                    cGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAatat                                                            >  2:3329333/1‑100 (MQ=255)
                                     gAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAatat                                                            >  2:1982670/1‑99 (MQ=255)
                                      aaTTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCg                                                          <  1:3929411/100‑1 (MQ=255)
                                        ttACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGcc                                                        <  1:1603791/100‑1 (MQ=255)
                                           cAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgg                                                     >  1:2169038/1‑100 (MQ=255)
                                             gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgggg                                                   >  2:674538/1‑100 (MQ=255)
                                             gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTggg                                                    >  1:3545883/1‑99 (MQ=255)
                                               gATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCt                                                 <  1:1657275/100‑1 (MQ=255)
                                                 tataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCttt                                               <  2:3183627/100‑1 (MQ=255)
                                                   taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTCACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                             >  2:3608682/1‑100 (MQ=255)
                                                   taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGCATAGTTTGAGTTCGAACGCTGGCCTCAGGTTGATAGACATATCGcc                                                        >  2:2899619/1‑89 (MQ=255)
                                                    acacATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGt                                            >  2:1559357/1‑100 (MQ=255)
                                                     cacaTTTCATGCCTACCTGATAAGAAGAGACTAATATTCGAATAGTCTGGGTTCGAACGCTGGCCTCAGGTTGATATAAATATCTTCTGGGACTTgtgtt                                           <  2:135604/100‑2 (MQ=255)
                                                     cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                           >  2:2230916/1‑100 (MQ=255)
                                                     cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                           >  2:2655726/1‑100 (MQ=255)
                                                     cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                           >  1:3917253/1‑100 (MQ=255)
                                                       caTGTCTTGCCTCTCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCa                                         <  2:1994814/100‑1 (MQ=255)
                                                       caTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCa                                         >  2:3359040/1‑100 (MQ=255)
                                                         tGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATc                                       >  2:3943584/1‑100 (MQ=255)
                                                          gTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCt                                      >  1:2728790/1‑100 (MQ=255)
                                                          gTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCt                                      <  2:4116600/100‑1 (MQ=255)
                                                          gTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGCCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATc                                       >  2:2994240/1‑99 (MQ=255)
                                                           tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg                                     >  1:2472415/1‑100 (MQ=255)
                                                           tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg                                     >  1:1086251/1‑100 (MQ=255)
                                                             ttGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGa                                   >  1:2419469/1‑100 (MQ=255)
                                                              tGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGaa                                  >  2:392464/1‑100 (MQ=255)
                                                               gCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAc                                 >  2:1157260/1‑100 (MQ=255)
                                                               gCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAc                                 <  1:1943301/100‑1 (MQ=255)
                                                                ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                                >  1:2860773/1‑100 (MQ=255)
                                                                ccTCCCTACAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                                >  1:4026752/1‑100 (MQ=255)
                                                                 cTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCt                               >  2:4120564/1‑100 (MQ=255)
                                                                 cTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCt                               >  1:2009618/1‑100 (MQ=255)
                                                                  tCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCTTCTGGAACCTc                              >  1:2016118/1‑100 (MQ=255)
                                                                   cccTGCAAAGAAGAGACTAACATTCGAATATTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcg                             <  2:1152141/100‑1 (MQ=255)
                                                                   cccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcg                             <  1:1180983/100‑1 (MQ=255)
                                                                    ccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc                            <  2:519860/100‑1 (MQ=255)
                                                                    ccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc                            >  1:1089459/1‑100 (MQ=255)
                                                                    ccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc                            >  1:1724179/1‑100 (MQ=255)
                                                                    ccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc                            <  2:1962998/100‑1 (MQ=255)
                                                                    ccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc                            >  1:3091901/1‑100 (MQ=255)
                                                                    ccTGCAAAGAAGAGACTAACATTCGAATAGTTTCGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc                            >  2:3856526/1‑100 (MQ=255)
                                                                     cTGCACAGAAGAGACTAACATTCGAATAGTTTGTGTTCGAACGATGGCCTCAGTTTGATAGAGATATCGTCTGGCGCTTTTGTACATCTGGAACCCcgcg                           >  2:2519223/1‑100 (MQ=255)
                                                                     cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgcg                           <  1:3087043/100‑1 (MQ=255)
                                                                      tGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGa                          >  2:2658342/1‑100 (MQ=255)
                                                                       gCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGaa                         >  1:1164644/1‑100 (MQ=255)
                                                                         aaaGCAGAGACTAACATTTGAATAGTTTGGGTTCGAAAGCTGCACTCAGGTTTATAGCAATATCGCCTGGGGCTTTTGTCCATCAGGaa                                  >  2:3511943/1‑89 (MQ=255)
                                                                         aaaGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATg                       <  2:3450021/100‑1 (MQ=255)
                                                                         aaaGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATg                       <  1:3531046/100‑1 (MQ=255)
                                                                          aagaagAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGc                      >  2:2282571/1‑100 (MQ=255)
                                                                          aagaagAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGc                      <  1:2435603/100‑1 (MQ=255)
                                                                          aagaagAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGc                      <  1:3996884/100‑1 (MQ=255)
                                                                           agaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCt                     <  2:2848526/100‑1 (MQ=255)
                                                                           agaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCt                     >  1:991258/1‑100 (MQ=255)
                                                                            gaagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCtt                    >  1:2889830/1‑100 (MQ=255)
                                                                              agagACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTaa                  >  1:1363388/1‑100 (MQ=255)
                                                                               gggagtAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAAGATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAAc                 <  2:3492773/95‑1 (MQ=255)
                                                                                agaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACg                <  2:2326428/100‑1 (MQ=255)
                                                                                 gaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGc               >  1:2842773/1‑100 (MQ=255)
                                                                                  aCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc              <  1:3132505/100‑1 (MQ=255)
                                                                                   cTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCa             >  1:1470726/1‑100 (MQ=255)
                                                                                      aCATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGAc          <  2:2595595/100‑1 (MQ=255)
                                                                                      aCATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGAc          >  2:1741623/1‑100 (MQ=255)
                                                                                       cATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACa         >  1:2596410/1‑100 (MQ=255)
                                                                                        aTTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGTTTGATAGAAATATCGCCTAGGGCTTTTGTCCATCTGGAACCTCGCCAATGCTTAACGCAAGAc          >  2:4048556/1‑98 (MQ=255)
                                                                                        aTTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAg        >  1:2661666/1‑100 (MQ=255)
                                                                                         ttCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGc       >  1:229642/1‑100 (MQ=255)
                                                                                         ttCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGc       >  1:1569907/1‑100 (MQ=255)
                                                                                          tCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGcc      <  1:392467/100‑1 (MQ=255)
                                                                                              aTTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCaa  <  1:619009/100‑1 (MQ=255)
                                                                                                |                                                                                                 
ATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAA  >  NC_000913/2887256‑2887449

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: