Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,887,352 T→A 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.3% 443.4 / ‑5.0 138intergenic (‑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/0);  major base A (74/61);  minor base G (0/1);  total (76/62)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 138 total aligned reads displayed.

ATCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCA  >  NC_000913/2887256‑2887448
                                                                                                |                                                                                                
atCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGtt                                                                                               >  2:3183344/1‑100 (MQ=255)
 tCAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGttt                                                                                              <  2:874330/100‑1 (MQ=255)
 tCAGCGGGCTATTTTTCGATGCTGCTCATAAATAGCGAATTATAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGttt                                                                                              <  1:1548109/100‑1 (MQ=255)
  cAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTAAAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTg                                                                                             <  1:841304/100‑1 (MQ=255)
   aGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTgg                                                                                            >  2:1874744/1‑100 (MQ=255)
    gCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTggg                                                                                           <  1:3004379/100‑1 (MQ=255)
        gCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTATGGGt                                                                                          <  1:83644/97‑1 (MQ=255)
         cTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGa                                                                                      >  1:2507571/1‑100 (MQ=255)
           aTTTTTCGATGCTACTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGCAGCGACTAACATTCGAATAGTTTGGGTTCGAAc                                                                                    >  1:2759082/1‑100 (MQ=255)
               ttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTg                                                                                <  1:2635405/100‑1 (MQ=255)
                tCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTgg                                                                               <  1:328116/100‑1 (MQ=255)
                tCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTgg                                                                               <  1:327072/100‑1 (MQ=255)
                  gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                             <  2:2300423/100‑1 (MQ=255)
                  gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                             <  2:3083146/100‑1 (MQ=255)
                  gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                             <  1:2455220/100‑1 (MQ=255)
                   aTGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCCCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                             >  1:1151185/1‑99 (MQ=255)
                    tgctgcTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTc                                                                           >  2:1362746/1‑100 (MQ=255)
                      ctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAg                                                                         >  1:1176275/1‑100 (MQ=255)
                       tgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAgg                                                                        <  1:2389827/100‑1 (MQ=255)
                        tctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGt                                                                       <  2:141538/99‑1 (MQ=255)
                         ctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGta                                                                      >  1:141536/1‑99 (MQ=255)
                          tCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGtt                                                                      <  2:185380/99‑1 (MQ=255)
                          tCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTg                                                                     >  2:2242033/1‑100 (MQ=255)
                           cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGa                                                                    >  2:1742307/1‑100 (MQ=255)
                            aTAAATATCGAATTACAGTGATATACACATTTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGAt                                                                   <  2:1719304/100‑1 (MQ=255)
                            aTAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGAt                                                                   <  2:2925533/100‑1 (MQ=255)
                            aTAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGAt                                                                   <  1:2532120/100‑1 (MQ=255)
                             tAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATa                                                                  >  2:1164377/1‑100 (MQ=255)
                             tAAAGATCGCATTACAGTGATAGACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACAGGCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGAt                                                                   <  1:2531941/99‑1 (MQ=255)
                              aaaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAg                                                                 <  1:2734447/100‑1 (MQ=255)
                              aaaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAg                                                                 >  1:983531/1‑100 (MQ=255)
                                 tatCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGaaa                                                              >  2:2989444/1‑100 (MQ=255)
                                  atCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAat                                                             >  2:1231652/1‑100 (MQ=255)
                                   tCGAATTACAGTGATATATACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAata                                                            <  2:211938/100‑1 (MQ=255)
                                   tCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAata                                                            >  2:2999594/1‑100 (MQ=255)
                                     gAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATc                                                          >  1:2029250/1‑100 (MQ=255)
                                      aaTTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATc                                                          <  2:3300433/99‑1 (MQ=255)
                                        ttACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGcc                                                       >  2:202695/1‑100 (MQ=255)
                                        gtACAGTGATATACACATGTCTTGCCTGCCGGCAAAGAAGAGACTAACATTCGAATAGTGTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGcc                                                       <  1:2011424/99‑1 (MQ=255)
                                         tACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCt                                                      >  2:3344453/1‑100 (MQ=255)
                                         tACAGAGATATAAACATGTCTTGCCTAAATGACAGGAAGAGACTAACATTCGGATAGTTTGGGTTCGAACGCTGGCAAGATGTTTATAGAGCTATCgtc                                                       <  1:2379107/99‑3 (MQ=255)
                                          aCAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTg                                                     >  1:2240210/1‑100 (MQ=255)
                                          aCAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTg                                                     <  2:576711/100‑1 (MQ=255)
                                          aCAGTGATATACACATGTCTCGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTg                                                     >  1:382456/1‑100 (MQ=255)
                                           cAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgg                                                    >  1:3399601/1‑100 (MQ=255)
                                            aGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTggg                                                   <  1:1627664/100‑1 (MQ=255)
                                             gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgggg                                                  >  1:171092/1‑100 (MQ=255)
                                              tGATATACACATGTCTTGTCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGc                                                 <  1:131866/100‑1 (MQ=255)
                                              tGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGc                                                 <  1:1012158/100‑1 (MQ=255)
                                                 tataCACATGTCTTGCCTCCCTGCAAAGAAGGGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTTGGGCttt                                              >  1:1089049/1‑100 (MQ=255)
                                                 tataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCttt                                              >  2:2598553/1‑100 (MQ=255)
                                                 tataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCttt                                              <  2:1295678/100‑1 (MQ=255)
                                                  ataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACAGTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCttt                                              <  1:65901/99‑1 (MQ=255)
                                                   taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                            <  2:2684148/100‑1 (MQ=255)
                                                   taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                            >  2:1920365/1‑100 (MQ=255)
                                                    acacATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTTGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGt                                           <  1:1869607/100‑1 (MQ=255)
                                                    acacATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGt                                           <  1:888939/100‑1 (MQ=255)
                                                     cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                          >  2:1536145/1‑100 (MQ=255)
                                                     cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                          >  2:357828/1‑100 (MQ=255)
                                                       caTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATAGCCTGGGGCTTTTGTc                                          >  1:1513307/1‑98 (MQ=255)
                                                         tGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATc                                      >  2:66607/1‑100 (MQ=255)
                                                          gTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCt                                     >  1:710573/1‑100 (MQ=255)
                                                            cTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGACCTCCGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg                                    >  1:612453/1‑99 (MQ=255)
                                                             ttGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGa                                  >  1:1914377/1‑100 (MQ=255)
                                                             ttGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGa                                  <  2:256699/100‑1 (MQ=255)
                                                             ttGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGa                                  <  1:3044753/100‑1 (MQ=255)
                                                             ttGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGa                                  <  2:272213/100‑1 (MQ=255)
                                                              tGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGGTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGaa                                 >  1:1202472/1‑100 (MQ=255)
                                                              tGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGCGTTCGAACGCTGTCCTCAGGTTGATAGAAATATCGCCTCGGGCTTTTGTCCATCTGGaa                                 <  1:861535/100‑1 (MQ=255)
                                                              tGCCTCCATGCGAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGaa                                 >  1:1823686/1‑100 (MQ=255)
                                                               gCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAc                                <  1:2885523/100‑1 (MQ=255)
                                                                ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                               >  2:815261/1‑100 (MQ=255)
                                                                ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                               >  1:2593419/1‑100 (MQ=255)
                                                                   cccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAc                                >  1:3386309/1‑96 (MQ=255)
                                                                    ccTGCAAAGAAGAGACTAACCTTCGAACAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc                           >  1:2118196/1‑100 (MQ=255)
                                                                    ccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc                           <  1:2680600/100‑1 (MQ=255)
                                                                     cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTGGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTTTCCATCTGGAACCTcgcg                          >  1:1855959/1‑100 (MQ=255)
                                                                      tGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTAGAACATCACGa                         >  1:246580/1‑100 (MQ=255)
                                                                       gCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGaa                        <  2:2088269/100‑1 (MQ=255)
                                                                        cAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAAt                       >  1:1023047/1‑100 (MQ=255)
                                                                        cAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAAt                       >  2:2453000/1‑100 (MQ=255)
                                                                         aaaGAAGAGACTAACATTTGAATCGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCACCGGGGGCTTTTGTCCATCTGGAACCTCACGAAt                       >  1:2065124/1‑99 (MQ=255)
                                                                         aaaGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATg                      <  2:3399653/100‑1 (MQ=255)
                                                                         aaaGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATg                      >  1:811900/1‑100 (MQ=255)
                                                                             aagaGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTa                  <  2:2506380/100‑1 (MQ=255)
                                                                               gagaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAAc                <  1:1363117/100‑1 (MQ=255)
                                                                                agaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACg               <  1:1943606/100‑1 (MQ=255)
                                                                                 gaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGc              >  1:480526/1‑100 (MQ=255)
                                                                                   cTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCa            <  2:899308/100‑1 (MQ=255)
                                                                                   cTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCa            <  2:1230478/100‑1 (MQ=255)
                                                                                   cTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCa            >  2:1763163/1‑100 (MQ=255)
                                                                                    tAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAg           >  1:862926/1‑100 (MQ=255)
                                                                                      acacTCAAATAGTTCGGGTTTGAGCGCTGGCCCCAGTTTGATAGAAATATCGCCTTGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGAc         <  2:2519533/100‑1 (MQ=255)
                                                                                      aCATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGAc         <  2:1452537/100‑1 (MQ=255)
                                                                                        aTTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAg       <  1:1245294/100‑1 (MQ=255)
                                                                                         ttCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGc      >  2:796255/1‑100 (MQ=255)
                                                                                          tcaaaTAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGcc     <  2:3397953/97‑1 (MQ=255)
                                                                                            gAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTc   <  1:319198/100‑1 (MQ=255)
                                                                                             aaTAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCa  >  1:366621/1‑100 (MQ=255)
                                                                                             aaTAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCACGAATGCTTAACGCCAGACAGCCTc   >  1:13884/1‑99 (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: