Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,881,050 T→C 83.0% intergenic (‑2/+1) casD ← / ← casC CRISP RNA (crRNA) containing Cascade antiviral complex protein/CRISP RNA (crRNA) containing Cascade antiviral complex protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,881,0500TC83.0% 2166.8 / inf 880intergenic (‑2/+1)casD/casCCRISP RNA (crRNA) containing Cascade antiviral complex protein/CRISP RNA (crRNA) containing Cascade antiviral complex protein
Reads supporting (aligned to +/- strand):  ref base T (69/81);  new base C (359/370);  total (428/452)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

Only 100 of 883 total aligned reads displayed.

TCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCC  >  NC_000913/2880903‑2881195
                                                                                                                                                   |                                                                                                                                                 
  aCGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCa                                                                                                                                                <  1:353054/149‑1 (MQ=255)
       ggATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCt                                                                                                                                           >  4:237792/1‑149 (MQ=255)
        gATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTc                                                                                                                                          >  3:4819413/1‑149 (MQ=255)
         aTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCg                                                                                                                                         >  1:2025045/1‑149 (MQ=255)
         aTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATACTCACGCCTCg                                                                                                                                         <  2:6125731/149‑1 (MQ=255)
             cAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCat                                                                                                                                     >  1:342778/1‑149 (MQ=255)
             cAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCat                                                                                                                                     <  3:6914097/149‑1 (MQ=255)
             cAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCat                                                                                                                                     >  4:6007789/1‑149 (MQ=255)
             cAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCat                                                                                                                                     >  1:5652886/1‑149 (MQ=255)
                  cAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGAACAAATAAGATCTCATGCTCACGCCTCGCCAttatta                                                                                                                                <  4:6736717/149‑1 (MQ=255)
                          gAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCa                                                                                                                        <  2:4740901/149‑1 (MQ=255)
                          gAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCa                                                                                                                        <  4:4634711/149‑1 (MQ=255)
                           aGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAg                                                                                                                       <  3:5227799/149‑1 (MQ=255)
                              aGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGAt                                                                                                                    >  3:2498119/1‑149 (MQ=255)
                               gCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGAtt                                                                                                                   <  2:1094614/149‑1 (MQ=255)
                                cccTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGAttt                                                                                                                  >  1:3907490/1‑149 (MQ=255)
                                cccTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGAttt                                                                                                                  >  3:3040203/1‑149 (MQ=255)
                                cccTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGAttt                                                                                                                  >  1:282572/1‑149 (MQ=255)
                                 ccTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGAtttt                                                                                                                 >  4:1634283/1‑149 (MQ=255)
                                  cTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTa                                                                                                                >  2:6926961/1‑149 (MQ=255)
                                   taataaCCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTaa                                                                                                               <  1:4182547/149‑1 (MQ=255)
                                   taataaCCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTaa                                                                                                               <  1:2302368/149‑1 (MQ=255)
                                         cccGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTc                                                                                                         <  1:533617/149‑1 (MQ=255)
                                         cccGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTc                                                                                                         >  4:1791819/1‑149 (MQ=255)
                                           cGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTa                                                                                                       <  1:5864304/149‑1 (MQ=255)
                                            gCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTa                                                                                                                <  2:6538301/139‑1 (MQ=255)
                                             cTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTaaa                                                                                                     <  4:1372803/149‑1 (MQ=255)
                                             cTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTaaa                                                                                                     <  4:3670232/149‑1 (MQ=255)
                                             cTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTaaa                                                                                                     <  4:443575/149‑1 (MQ=255)
                                              ttCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAg                                                                                                    <  3:1826792/149‑1 (MQ=255)
                                                cGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTa                                                                                                  <  2:3676860/149‑1 (MQ=255)
                                                cGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTa                                                                                                  <  2:2090944/149‑1 (MQ=255)
                                                   gTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGc                                                                                               >  1:4391892/1‑149 (MQ=255)
                                                    tCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCa                                                                                              >  1:6198249/1‑149 (MQ=255)
                                                      ggAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCAtt                                                                                            >  4:7001682/1‑149 (MQ=255)
                                                        aaaTCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTg                                                                                          >  2:997835/1‑149 (MQ=255)
                                                             ttCCGGTAGGTCGCGTTCCTTGAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTaa                                                                                     <  2:5080585/149‑1 (MQ=255)
                                                              tCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAAc                                                                                    <  2:2880425/149‑1 (MQ=255)
                                                                cGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCATGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCACGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAAcgt                                                                                  <  2:2574931/149‑3 (MQ=255)
                                                                 ggTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGtt                                                                                                          <  1:6565829/124‑1 (MQ=255)
                                                                  gTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCGAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGa                                                                                <  4:2608578/149‑1 (MQ=255)
                                                                  gTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGa                                                                                >  1:1035757/1‑149 (MQ=255)
                                                                  gTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGa                                                                                <  4:3343137/149‑1 (MQ=255)
                                                                   tAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAg                                                                               >  2:3869488/1‑149 (MQ=255)
                                                                   tAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAg                                                                               <  3:776002/149‑1 (MQ=255)
                                                                     ggTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCa                                                                             <  2:5872364/149‑1 (MQ=255)
                                                                     ggTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCa                                                                             <  1:4550769/149‑1 (MQ=255)
                                                                      gTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAACAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAg                                                                            >  3:1203336/1‑149 (MQ=255)
                                                                       tCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAAc                                                                                                              <  3:2212919/114‑1 (MQ=255)
                                                                       tCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGt                                                                           >  1:6800134/1‑149 (MQ=255)
                                                                       tCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGt                                                                           >  3:3834931/1‑149 (MQ=255)
                                                                         gcgTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCCTTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAg                                                                                                 <  2:5180691/125‑1 (MQ=255)
                                                                            ttCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTg                                                                      <  2:6330614/149‑1 (MQ=255)
                                                                             tCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTgg                                                                     <  2:2289091/149‑1 (MQ=255)
                                                                             tCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTgg                                                                     <  1:2495328/149‑1 (MQ=255)
                                                                              ccTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTggg                                                                    >  1:1650052/1‑149 (MQ=255)
                                                                                ttCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGc                                                                              >  2:4126939/1‑137 (MQ=255)
                                                                                ttCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTc                                                                                                                                          >  1:773231/1‑77 (MQ=255)
                                                                                ttCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTaa                                                                                                               >  3:1651667/1‑104 (MQ=255)
                                                                                ttCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTc                                                                  >  4:2105532/1‑149 (MQ=255)
                                                                                ttCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTc                                                                  >  4:5365343/1‑149 (MQ=255)
                                                                                ttCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTc                                                                  >  3:422117/1‑149 (MQ=255)
                                                                                  cAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTa                                                                >  2:778490/1‑149 (MQ=255)
                                                                                   aaaGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTAc                                                               >  1:1289127/1‑149 (MQ=255)
                                                                                   aaaGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTAc                                                               >  3:4021626/1‑149 (MQ=255)
                                                                                    aaGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACa                                                              >  3:3216782/1‑149 (MQ=255)
                                                                                        tCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAg                                                          >  2:4532311/1‑149 (MQ=255)
                                                                                         cGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGa                                                         >  2:5424702/1‑149 (MQ=255)
                                                                                               cccccAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCt                                                   <  4:6564208/149‑1 (MQ=255)
                                                                                                ccccAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTg                                                  >  2:6804181/1‑149 (MQ=255)
                                                                                                   cAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAAt                                               <  3:5975900/149‑1 (MQ=255)
                                                                                                   cAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAAt                                               <  3:5736608/149‑1 (MQ=255)
                                                                                                    aGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAAtt                                              <  1:5167668/149‑1 (MQ=255)
                                                                                                      gCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCa                                                                                                                                      >  2:1655000/1‑59 (MQ=255)
                                                                                                      gCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTgc                                            <  1:6392971/149‑1 (MQ=255)
                                                                                                        ttGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTgcgc                                          >  4:3089965/1‑149 (MQ=255)
                                                                                                         tGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCgca                                         <  1:3602580/149‑1 (MQ=255)
                                                                                                            aTTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCgcagca                                      <  4:4937908/149‑1 (MQ=255)
                                                                                                                 cccAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTcc                                 <  3:442736/149‑1 (MQ=255)
                                                                                                                   cAGCAAGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGt                               >  3:277052/1‑149 (MQ=255)
                                                                                                                      cAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCa                            <  4:4020/149‑1 (MQ=255)
                                                                                                                       aaGCCGCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTATGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAg                           >  1:4910521/1‑149 (MQ=255)
                                                                                                                        aGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGa                          <  1:69083/149‑1 (MQ=255)
                                                                                                                            gCAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGa                          <  3:3533231/145‑1 (MQ=255)
                                                                                                                             cAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCata                     <  2:1344146/149‑1 (MQ=255)
                                                                                                                             cAAGATCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCata                     <  3:1466999/149‑1 (MQ=255)
                                                                                                                                 aTCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCa                 >  4:6163482/1‑149 (MQ=255)
                                                                                                                                 aTCAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCa                 >  4:1158481/1‑149 (MQ=255)
                                                                                                                                   cAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCAtt               <  2:3467511/149‑1 (MQ=255)
                                                                                                                                   cAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCAtt               <  3:671305/149‑1 (MQ=255)
                                                                                                                                   cAAATAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCAtt               <  4:5765193/149‑1 (MQ=255)
                                                                                                                                     aaTAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTgg             >  3:6495231/1‑149 (MQ=255)
                                                                                                                                      aTAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAg                                                                               >  3:6115713/1‑82 (MQ=255)
                                                                                                                                      aTAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATCTTAACTGTTCTaaa                                                                                                     <  3:536993/60‑1 (MQ=255)
                                                                                                                                       tAAGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCa           >  4:5106040/1‑149 (MQ=255)
                                                                                                                                        aaGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGa                          <  3:4619193/133‑1 (MQ=255)
                                                                                                                                         aGATCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAAc         <  4:2849939/149‑1 (MQ=255)
                                                                                                                                           aTCTCATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAAcgc       >  1:3120176/1‑149 (MQ=255)
                                                                                                                                              tcATGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTACAGACCATATCCATTGGCAACGCGAt    <  2:3148212/149‑1 (MQ=255)
                                                                                                                                                aTGCTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATcc  >  4:5594214/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGCAGCTCCGTTCAGACCATATCCATTGGCAACGCGATCC  >  NC_000913/2880903‑2881195

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: