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.2% 379.5 / ‑4.3 117intergenic (‑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);  new base A (52/64);  total (53/64)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

Only 100 of 117 total aligned reads displayed.

CAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTC  >  NC_000913/2887258‑2887447
                                                                                              |                                                                                               
cAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTATAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTg                                                                                            <  1:412677/100‑1 (MQ=255)
cAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTg                                                                                            >  2:2591304/1‑100 (MQ=255)
 aGCGGGCTATTTTTTGATGCTGCTTATAAACATCGAATTACATTGATATACACATGTCTTGCCTCCATGCAAAGAAGAGACTTACCATCGAATAGTTTg                                                                                            <  2:427402/99‑1 (MQ=255)
  gCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTggg                                                                                          >  1:947673/1‑100 (MQ=255)
   cGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTCCAGTGATATACACCTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGCATAGTTTGAGt                                                                                         <  2:892813/100‑1 (MQ=255)
    gggCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGtt                                                                                        >  1:870977/1‑100 (MQ=255)
    gggCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGtt                                                                                        <  2:1726377/100‑1 (MQ=255)
      gCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCg                                                                                      <  1:1936735/100‑1 (MQ=255)
       cTATTTTTCGATGCTGCTCATAAATATCGAATTCCAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAATTTGGGTTCGa                                                                                     >  1:482896/1‑100 (MQ=255)
             ttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATTTCTTGTCTCCCTGCAAAGAAGAGACTAACATTCTAATAGTTTGGGTTCGAACGCTg                                                                               <  2:627138/100‑1 (MQ=255)
             ttCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTg                                                                               <  2:1746049/100‑1 (MQ=255)
              tCGATGCTGCTCATAAATATCGAATTATAGTGATATACAGATGTCTTGCCTTCCTGCAAAGAAGAGACTAACTTTCGAATAGTTTGGGTTCGAACGCTgg                                                                              <  2:1804015/100‑1 (MQ=255)
              tCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTgg                                                                              >  2:2359231/1‑100 (MQ=255)
              tCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTgg                                                                              <  1:1006199/100‑1 (MQ=255)
               cGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGc                                                                             <  2:1395321/100‑1 (MQ=255)
               cGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGc                                                                             >  1:2849507/1‑100 (MQ=255)
               cGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTACCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTAGAACGCTGGc                                                                             >  2:1927416/1‑100 (MQ=255)
               cGATGCTGCTCATAAATATCAAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGc                                                                             <  2:1363473/100‑1 (MQ=255)
                gATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTgg                                                                              >  1:919538/1‑98 (MQ=255)
                 aTGCTGCTCATAAATATCGAATTACAGTGATATACACATTTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGcc                                                                            <  1:251163/99‑1 (MQ=255)
                      gctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGt                                                                      >  1:1259612/1‑100 (MQ=255)
                      gctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGt                                                                      <  1:1751248/100‑1 (MQ=255)
                      gctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGt                                                                      <  2:2938273/100‑1 (MQ=255)
                       ctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGtt                                                                     >  1:1903057/1‑100 (MQ=255)
                        tCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTg                                                                    <  2:1506906/100‑1 (MQ=255)
                        tCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTg                                                                    <  2:1506861/100‑1 (MQ=255)
                         cATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTg                                                                    <  2:2896150/99‑1 (MQ=255)
                          aTAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGAt                                                                  <  1:1907862/100‑1 (MQ=255)
                           tAAATATCGAATTACAGTGTTATACACATGTCTTGCCTTCCTGCAAAGAAGAGACTAACATTCGAATAGTTGGGGTTCGAACGCTGGCCTCAGGTTGAt                                                                  <  1:760078/99‑1 (MQ=255)
                           tAAATATCGAATTACAGTGATATACACATGTCTTTACTCCCTGCAAAGAAGAGACTAAAATTCGAATAGTTGGGGATCGAACGCTGGCCTCAGGTTGATa                                                                 >  2:2150524/1‑100 (MQ=255)
                           tAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCCGCAAAGAAGAGACAAAAATTCAAATATATTGGGTTCGAACGCTGGCCTCAGGTTGATa                                                                 >  1:1821781/1‑100 (MQ=255)
                            aaaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAg                                                                <  2:530408/100‑1 (MQ=255)
                               tttctaaTTACTGTGATATACACATGTCGTGCTGCTCTGCAAATAAGTGGCTTACATTCGGTTAATTTGGGTTCGATCGCTGTCCTCAGGTTGATAGaa                                                              <  2:449114/94‑1 (MQ=255)
                                atCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAat                                                            >  2:794123/1‑100 (MQ=255)
                                atCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCATCAGGTTGATAGAAat                                                            >  1:854246/1‑100 (MQ=255)
                                 tCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAata                                                           >  2:237844/1‑100 (MQ=255)
                                  cGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAatat                                                          >  2:934959/1‑100 (MQ=255)
                                   gAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATc                                                         <  1:1068324/100‑1 (MQ=255)
                                   gAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATc                                                         <  1:1522577/100‑1 (MQ=255)
                                     aTTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGc                                                       <  2:1919771/100‑1 (MQ=255)
                                     aTTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGc                                                       <  2:1919735/100‑1 (MQ=255)
                                      ttACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGcc                                                      <  1:2570571/100‑1 (MQ=255)
                                         cAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgg                                                   >  1:831528/1‑100 (MQ=255)
                                           gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGGCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgggg                                                 >  1:434515/1‑100 (MQ=255)
                                           gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgggg                                                 <  1:2763072/100‑1 (MQ=255)
                                           gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTgggg                                                 >  2:1564999/1‑100 (MQ=255)
                                           gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTAGGGTTCGAACGCTGGCATCAGGTTGATAGAAATATCGCCTggg                                                  >  2:2442305/1‑99 (MQ=255)
                                             gATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCt                                               >  1:2417021/1‑100 (MQ=255)
                                             gATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCt                                               <  2:831290/100‑1 (MQ=255)
                                               tataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCttt                                             <  2:2832886/100‑1 (MQ=255)
                                                 taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCtttt                                            >  2:512391/1‑99 (MQ=255)
                                                 taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTg                                           >  2:2219187/1‑100 (MQ=255)
                                                  acacATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGt                                          <  2:1428029/100‑1 (MQ=255)
                                                  acacATGTCTCGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGACCGATGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGt                                          >  2:2138691/1‑100 (MQ=255)
                                                   cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                         >  2:3041061/1‑100 (MQ=255)
                                                   cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                         >  1:332985/1‑100 (MQ=255)
                                                   cacaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                         >  2:3026488/1‑100 (MQ=255)
                                                    acaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTcc                                        >  1:1111486/1‑100 (MQ=255)
                                                     caTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCa                                       >  1:720072/1‑100 (MQ=255)
                                                       tGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATc                                     <  2:482893/100‑1 (MQ=255)
                                                       tGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATc                                     <  2:602125/100‑1 (MQ=255)
                                                        gTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCt                                    >  2:2250376/1‑100 (MQ=255)
                                                        gTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCt                                    <  1:2297386/100‑1 (MQ=255)
                                                         tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg                                   <  1:2822802/100‑1 (MQ=255)
                                                         tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg                                   <  1:2826535/100‑1 (MQ=255)
                                                         tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTg                                   <  1:2826835/100‑1 (MQ=255)
                                                          nttGCCTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACACTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTgg                                  >  1:1167584/2‑100 (MQ=255)
                                                              ccTCCCTGCAAAGAAGAGACTAACATTCGAATTGCTTAGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                              >  1:1935065/1‑100 (MQ=255)
                                                              ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTTGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                              <  1:967313/100‑1 (MQ=255)
                                                              ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                              >  1:2058732/1‑100 (MQ=255)
                                                              ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                              <  2:677220/100‑1 (MQ=255)
                                                              ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                              >  1:1571687/1‑100 (MQ=255)
                                                              ccTCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAcc                              <  2:2633322/100‑1 (MQ=255)
                                                               cTCCCTGTAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCt                             <  1:674661/100‑1 (MQ=255)
                                                                tCCCTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTc                            <  1:378550/100‑1 (MQ=255)
                                                                 cccTGCAAAGAAGAGACTAACATTCGAATAGTTTGTGTTCGAACGCTGGCATCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcg                           >  2:857899/1‑100 (MQ=255)
                                                                 cccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcg                           >  2:117672/1‑100 (MQ=255)
                                                                  ccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgc                          <  2:947667/100‑1 (MQ=255)
                                                                  ccTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGACCTCAGGTTGATAGAAATATCACCTGGGGCTTTTGTCCATCTGGAACCTcgc                          >  1:32784/1‑100 (MQ=255)
                                                                   cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgcg                         <  1:848795/100‑1 (MQ=255)
                                                                   cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgcg                         >  2:1156230/1‑100 (MQ=255)
                                                                   cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTcgcg                         <  1:192190/100‑1 (MQ=255)
                                                                   cTGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATAGCCTGGGGCTTTTGTCCATCTGGAACCTcgcg                         >  1:999892/1‑100 (MQ=255)
                                                                    tGCAAAGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGa                        >  1:2367744/1‑100 (MQ=255)
                                                                     gCAAAGCAGAGACTAACATTCGAATAGTTTGGGTTCGAAAGCTGGCCACAGGTTGATAGAAATATCGACTGGGGATTTTGACGATCTGGAAACTCGCGaa                       >  2:2975340/1‑100 (MQ=255)
                                                                     gCAAAGAAGAGACTAACATTGGAATAGATTGGGATCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATATGGAACCTcccaa                        >  2:689596/1‑95 (MQ=255)
                                                                       aaaGAAGAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATg                     <  2:807149/100‑1 (MQ=255)
                                                                        aagaagAGACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGc                    <  1:121374/100‑1 (MQ=255)
                                                                          gaagaGACTAACATTCGAATAGTTGGGGTTCGAACGCTGGCCTCATGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCtt                  <  2:1259600/100‑1 (MQ=255)
                                                                            agagACTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTaa                <  1:2854974/100‑1 (MQ=255)
                                                                             gagaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAAc               <  1:954496/100‑1 (MQ=255)
                                                                               gaCTACCATTAGAATAGTTTGGGATAGAAGGCTGGCCTCAGGTTGATAGAAATATGGCCTGGGGCTTTTGTCGATCTGGAACTTAGCGAATGCTTAACGc             >  2:18560/1‑100 (MQ=255)
                                                                               gaCTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGc             >  1:2645121/1‑100 (MQ=255)
                                                                                 cTAACATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCa           <  2:389137/100‑1 (MQ=255)
                                                                                    aCATTCGAATAGTTTGGGTTCGATCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGAc        <  1:169151/100‑1 (MQ=255)
                                                                                    aCATTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGAc        <  1:1093289/100‑1 (MQ=255)
                                                                                      aTTCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAg      <  2:1956006/100‑1 (MQ=255)
                                                                                       ttCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGc     <  2:1305812/100‑1 (MQ=255)
                                                                                        tCGAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGcc    <  2:1290184/100‑1 (MQ=255)
                                                                                          gAATAGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTc  <  2:2714884/100‑1 (MQ=255)
                                                                                              |                                                                                               
CAGCGGGCTATTTTTCGATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTC  >  NC_000913/2887258‑2887447

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: