Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,884,579 T→G 5.5% R881R (AGA→CGA) ‡ ygcB ← Cascade complex anti‑viral R‑loop helicase‑annealase Cas3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,884,5790TG 5.5% 417.7 / 10.8 137R881R (AGA→CGA) ‡ygcBCascade complex anti‑viral R‑loop helicase‑annealase Cas3
Reads supporting (aligned to +/- strand):  ref base T (52/77);  new base G (2/6);  total (54/83)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

Only 100 of 140 total aligned reads displayed.

TGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGA  >  NC_000913/2884430‑2884728
                                                                                                                                                     |                                                                                                                                                       
tgtgAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTAATTTGGGATTTGCAGGGATTActct                                                                                                                                                         >  2:905494/1‑150 (MQ=255)
 gtgAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTACCGAAGTTATTTGGGATTTGCAGGGATGACTCTg                                                                                                                                                        >  1:937214/1‑150 (MQ=255)
       ataAAATCTCATACCGGGAAATTAAAAGAAGATGTACCTTGTGCACCTTCCCTACTTAAGTCGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATc                                                                                                                                                  <  2:222608/150‑1 (MQ=255)
       ataAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATc                                                                                                                                                  <  1:715611/150‑1 (MQ=255)
        taAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATcc                                                                                                                                                 <  1:1763848/150‑1 (MQ=255)
           aaTCTCATACCGTGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAATCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCtg                                                                                                                                              <  1:1483062/150‑2 (MQ=255)
           aaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCtg                                                                                                                                              <  1:1796441/150‑2 (MQ=255)
           aaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCtg                                                                                                                                              <  1:1770755/150‑2 (MQ=255)
           aaTCTCATACCGGGAAATTAAAAGAACATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCtt                                                                                                                                              <  1:1337740/150‑1 (MQ=255)
            aTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAAGGATGACCAGAGTCATCCCtgc                                                                                                                                             <  1:813856/150‑3 (MQ=255)
             tctcATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAATTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCTGC‑‑a                                                                                                                                          <  1:891786/150‑1 (MQ=255)
                    ccGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTg                                                                                                                                   >  2:694455/1‑150 (MQ=255)
                    ccGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTg                                                                                                                                   >  1:571121/1‑150 (MQ=255)
                     cGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGta                                                                                                                                 <  1:934217/150‑1 (MQ=255)
                      gggAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGta                                                                                                                                 <  2:1062970/150‑1 (MQ=255)
                      gggAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTAGCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGta                                                                                                                                 <  1:113821/150‑1 (MQ=255)
                        gAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGtata                                                                                                                               <  1:1046646/150‑1 (MQ=255)
                          aaTTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAgg                                                                                                                             <  2:1252493/150‑1 (MQ=255)
                           aTTAAAAGAATATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGt                                                                                                                            <  1:333716/150‑1 (MQ=255)
                           aTTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGt                                                                                                                            <  2:1242596/150‑1 (MQ=255)
                           aTTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGt                                                                                                                            <  2:491962/150‑1 (MQ=255)
                               aaaGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGtaata                                                                                                                        >  2:1096789/1‑150 (MQ=255)
                                 agaagaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAAc                                                                                                                      <  1:1106825/150‑1 (MQ=255)
                                 agaagaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAAc                                                                                                                      <  1:498174/150‑1 (MQ=255)
                                  gaagaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACa                                                                                                                     >  1:1735770/1‑150 (MQ=255)
                                   aagaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACaa                                                                                                                    <  2:1359401/150‑1 (MQ=255)
                                   aagaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACaa                                                                                                                    >  2:1570613/1‑150 (MQ=255)
                                    agaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAAt                                                                                                                   <  2:289957/150‑1 (MQ=255)
                                    agaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAACCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAAt                                                                                                                   <  2:392065/150‑1 (MQ=255)
                                      aTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATAc                                                                                                                 <  2:1607150/150‑1 (MQ=255)
                                      aTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATAc                                                                                                                 <  1:487262/150‑1 (MQ=255)
                                        gTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTg                                                                                                               <  1:1902440/150‑1 (MQ=255)
                                         tACATTGTGCACCTTCCCTACTTTAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGt                                                                                                              <  2:148991/150‑1 (MQ=255)
                                          aCATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGtt                                                                                                             >  1:236375/1‑150 (MQ=255)
                                           cATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGTGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTa                                                                                                            <  2:1446635/150‑1 (MQ=255)
                                                  gtttatCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGcc                                                                                                     >  1:1565447/6‑150 (MQ=255)
                                                      atCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAc                                                                                                 >  2:152562/2‑150 (MQ=255)
                                                       tCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAcc                                                                                                <  2:577967/150‑1 (MQ=255)
                                                       tCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAcc                                                                                                <  2:1590013/150‑1 (MQ=255)
                                                       tCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAcc                                                                                                <  1:987896/150‑1 (MQ=255)
                                                       tCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAcc                                                                                                <  2:760697/150‑1 (MQ=255)
                                                         ccTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACcca                                                                                              >  2:1901165/1‑150 (MQ=255)
                                                         ccTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCACGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACcca                                                                                              <  1:212391/150‑1 (MQ=255)
                                                           tACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatc                                                                                            <  1:1657799/150‑1 (MQ=255)
                                                           tACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatc                                                                                            >  2:1060497/1‑150 (MQ=255)
                                                           tACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatc                                                                                            <  2:164953/150‑1 (MQ=255)
                                                           tACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatc                                                                                            >  1:487731/1‑150 (MQ=255)
                                                           tACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatc                                                                                            >  2:1637596/1‑150 (MQ=255)
                                                                 aGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCAtccatcca                                                                                      <  1:85118/150‑1 (MQ=255)
                                                                 aGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCAtccatcca                                                                                      >  2:1821336/1‑150 (MQ=255)
                                                                  gTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCTGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAg                                                                                     <  1:1313131/150‑1 (MQ=255)
                                                                     gggATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGAtt                                                                                  >  1:1029032/1‑150 (MQ=255)
                                                                        aTAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGTCCCATCCATCCAGATTCTg                                                                               <  1:63669/150‑1 (MQ=255)
                                                                          aaaCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGtt                                                                             <  2:1697929/150‑1 (MQ=255)
                                                                           aaCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGTTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGttt                                                                            >  2:1304247/1‑150 (MQ=255)
                                                                            aCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTc                                                                           <  2:1539779/150‑1 (MQ=255)
                                                                            aCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTc                                                                           <  2:1113461/150‑1 (MQ=255)
                                                                             ccGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCTAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTcc                                                                          <  1:1373427/150‑1 (MQ=255)
                                                                                 tATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTc                                                                      >  2:193669/1‑150 (MQ=255)
                                                                                  aTTGGTCTTAATATCGTCCTTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTcc                                                                     <  2:1875609/150‑1 (MQ=255)
                                                                                   ttGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCa                                                                    >  1:1355253/1‑150 (MQ=255)
                                                                                   ttGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCa                                                                    >  2:1419198/1‑150 (MQ=255)
                                                                                    tGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAg                                                                   >  1:17294/1‑150 (MQ=255)
                                                                                              atCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAAc                                                         <  1:1146265/150‑1 (MQ=255)
                                                                                                cGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCTTCCAGATTCTGTTTCCCTTCCAGCCAAAGTAAccc                                                       >  1:599909/1‑150 (MQ=255)
                                                                                                 gTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCa                                                      <  1:1084501/150‑1 (MQ=255)
                                                                                                 gTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCa                                                      >  1:864016/1‑150 (MQ=255)
                                                                                                  tCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCAt                                                     >  2:1432800/1‑150 (MQ=255)
                                                                                                    aTTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATcg                                                   >  2:1160342/1‑149 (MQ=255)
                                                                                                     ttGATAACAATCATTCCCGAAGTTATTTGTGATTTGCAGGGATGACTCTGGTCATCGCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCtt                                                  <  1:135527/150‑1 (MQ=255)
                                                                                                         tAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATc                                              >  1:1290149/1‑150 (MQ=255)
                                                                                                         aaaCAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACACTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTACATc                                              >  1:336110/2‑150 (MQ=255)
                                                                                                          aaCAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATct                                             >  2:34616/1‑150 (MQ=255)
                                                                                                          aaCAATCATTCCCGAAGGTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATct                                             <  1:1404557/150‑1 (MQ=255)
                                                                                                           aCAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATcta                                            >  1:204087/1‑150 (MQ=255)
                                                                                                           aCAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCATTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTAACCCTGCCAGCCCCATGCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCAATCTTCATcta                                            >  2:51727/1‑150 (MQ=255)
                                                                                                                 aTTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTACCCCATCTTCATCTACGACtt                                      >  2:787366/1‑150 (MQ=255)
                                                                                                                     ccGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGa                                  >  2:654661/1‑150 (MQ=255)
                                                                                                                       gAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCAAGCCAAAGTAACCCATCTTCATCTACTACTTCAGaaa                                >  2:1091525/1‑150 (MQ=255)
                                                                                                                         aGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGaaaaa                              <  1:1096789/150‑1 (MQ=255)
                                                                                                                         aGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGaaaaa                              <  1:157858/150‑1 (MQ=255)
                                                                                                                           ttATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCACTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACt                            >  1:718068/1‑150 (MQ=255)
                                                                                                                             aTTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTAc                          >  2:165749/1‑150 (MQ=255)
                                                                                                                             aTTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTAc                          <  2:659175/150‑1 (MQ=255)
                                                                                                                                     tttGCAGGGATGACTCTTGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAg                  <  1:1907764/150‑1 (MQ=255)
                                                                                                                                       tGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCAGCCAGATTCTGTGTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGGTTCCAGGt                >  2:891696/1‑150 (MQ=255)
                                                                                                                                       tGCAGGGATGACTCTGGTCATCCCGTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCGGGt                <  1:207703/150‑1 (MQ=255)
                                                                                                                                        gaaGGGATGACCAGAGTCATCCCTGCAAATCCCATGTATAGGTATTAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACACATCTTCATCTACTACTTCAGAAATACTACTTTTCCAgg                 >  1:69322/3‑150 (MQ=255)
                                                                                                                                        gCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCCTCTTCATCGACTACTTCAGAAAAACTACGTTTCCAGGTg               >  2:1547552/1‑150 (MQ=255)
                                                                                                                                         cAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTCCCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGa              >  1:1147400/1‑150 (MQ=255)
                                                                                                                                          aGGGATGACCAGAGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGaa             >  1:1732536/1‑150 (MQ=255)
                                                                                                                                           gggATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAg            >  1:199658/1‑150 (MQ=255)
                                                                                                                                           gggATGACTCTGGTCATCCCGTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAg            <  1:905494/150‑1 (MQ=255)
                                                                                                                                            ggATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAgg           <  2:1200215/150‑1 (MQ=255)
                                                                                                                                             gATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAggg          <  2:1355253/150‑1 (MQ=255)
                                                                                                                                             gATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAggg          <  2:461986/150‑1 (MQ=255)
                                                                                                                                               tGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTa        <  1:165922/150‑1 (MQ=255)
                                                                                                                                                gACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTAc       <  2:17294/150‑1 (MQ=255)
                                                                                                                                                  ctctTGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACAt     <  1:650562/150‑1 (MQ=255)
                                                                                                                                                     tGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGa  >  2:439226/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
TGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGA  >  NC_000913/2884430‑2884728

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