Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,884,578 C→A 5.6% R881I (AGA→ATA) ‡ 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,5780CA 5.6% 431.5 / 10.7 137R881I (AGA→ATA) ‡ygcBCascade complex anti‑viral R‑loop helicase‑annealase Cas3
Reads supporting (aligned to +/- strand):  ref base C (52/77);  new base A (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 = 9.87e-01

Only 100 of 139 total aligned reads displayed.

TGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACAT  >  NC_000913/2884430‑2884725
                                                                                                                                                    |                                                                                                                                                     
tgtgAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTAATTTGGGATTTGCAGGGATTActct                                                                                                                                                      >  2:905494/1‑150 (MQ=255)
 gtgAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTg                                                                                                                                                     >  2:1459068/1‑150 (MQ=255)
 gtgAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTg                                                                                                                                                     >  2:1530768/1‑150 (MQ=255)
  tgAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTgg                                                                                                                                                    <  2:418056/150‑1 (MQ=255)
         aaaaTCTCATACCGTGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATccc                                                                                                                                             <  1:36240/150‑1 (MQ=255)
          aaaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCt                                                                                                                                            <  1:1124042/150‑1 (MQ=255)
           aaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCtg                                                                                                                                           <  1:1770755/150‑2 (MQ=255)
            aTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAAGGATGACCAGAGTCATCCCtgc                                                                                                                                          <  1:813856/150‑3 (MQ=255)
                cATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATcc                                                                                                                                    <  1:603799/150‑1 (MQ=255)
                 aTACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGAAGGAATGACCAGAGTCATCCCTGC‑‑aaatc                                                                                                                                   <  1:332267/150‑5 (MQ=255)
                  tACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATcccc                                                                                                                                  >  2:1907764/1‑150 (MQ=255)
                  tACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATcccc                                                                                                                                  >  2:1115060/1‑150 (MQ=255)
                    ccGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTg                                                                                                                                >  2:694455/1‑150 (MQ=255)
                     cGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGta                                                                                                                              <  1:934217/150‑1 (MQ=255)
                      gggAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGta                                                                                                                              <  2:1062970/150‑1 (MQ=255)
                        gAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGtata                                                                                                                            <  1:1776208/150‑1 (MQ=255)
                        gAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGtata                                                                                                                            <  1:1046646/150‑1 (MQ=255)
                           aTTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGt                                                                                                                         <  2:1242596/150‑1 (MQ=255)
                           aTTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGt                                                                                                                         <  2:491962/150‑1 (MQ=255)
                              aaaaGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGtaat                                                                                                                      >  2:232284/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:1570613/1‑150 (MQ=255)
                                    agaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAAt                                                                                                                <  2:289957/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)
                                           cATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTa                                                                                                         >  2:16084/1‑150 (MQ=255)
                                                  gtttatCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGcc                                                                                                  >  1:1565447/6‑150 (MQ=255)
                                                    ccTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAg                                                                                                <  2:168536/150‑1 (MQ=255)
                                                      atCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAc                                                                                              >  1:913372/2‑150 (MQ=255)
                                                       tCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAcc                                                                                             <  2:577967/150‑1 (MQ=255)
                                                       tCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAcc                                                                                             <  1:987896/150‑1 (MQ=255)
                                                       tCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAcc                                                                                             <  2:760697/150‑1 (MQ=255)
                                                       tCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAcc                                                                                             <  2:1590013/150‑1 (MQ=255)
                                                         ccTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCACGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACcca                                                                                           <  1:212391/150‑1 (MQ=255)
                                                           tACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatc                                                                                         >  1:487731/1‑150 (MQ=255)
                                                           tACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatc                                                                                         >  2:1637596/1‑150 (MQ=255)
                                                           tACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatc                                                                                         <  1:1657799/150‑1 (MQ=255)
                                                                 aGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCAtccatcca                                                                                   >  2:1821336/1‑150 (MQ=255)
                                                                  gTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCTGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAg                                                                                  <  1:1313131/150‑1 (MQ=255)
                                                                  gTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAg                                                                                  <  2:1068014/150‑1 (MQ=255)
                                                                     gggATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGAtt                                                                               >  1:1029032/1‑150 (MQ=255)
                                                                      ggATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTc                                                                              <  1:1750713/150‑1 (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)
                                                                                 ttttGGTCTTATTATCGTCATTGATAAGAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATATCAATACTGTTACCCTGCCAGACCCATCCATCCAGCTTCTGTTTCCCTTc                                                                   <  1:204104/148‑1 (MQ=255)
                                                                                 tATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTc                                                                   >  2:193669/1‑150 (MQ=255)
                                                                                  aTTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTcc                                                                  <  1:994143/150‑1 (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)
                                                                                        cTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCaa                                                            <  1:890576/150‑1 (MQ=255)
                                                                                              atCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAAc                                                      <  1:1146265/150‑1 (MQ=255)
                                                                                                cGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCTTCCAGATTCTGTTTCCCTTCCAGCCAAAGTAAccc                                                    >  1:599909/1‑150 (MQ=255)
                                                                                                 gTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCa                                                   >  1:864016/1‑150 (MQ=255)
                                                                                                 gTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCa                                                   <  1:1084501/150‑1 (MQ=255)
                                                                                                  tCATTGATAACAATCCTTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATTACAATACTGTTACCATGCCAGACCGCTCGATCCAGATTCTGTTTCCCTTCCAGCCCAAGTAACCCAt                                                  >  2:1309970/1‑150 (MQ=255)
                                                                                                     ttGATAACAATCATTCCCGAAGTTATTTGTGATTTGCAGGGATGACTCTGGTCATCGCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCtt                                               <  1:135527/150‑1 (MQ=255)
                                                                                                         tAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATc                                           <  1:1665822/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)
                                                                                                            cAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctac                                        <  2:1381923/150‑1 (MQ=255)
                                                                                                                 aTTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTACCCCATCTTCATCTACGACtt                                   >  2:787366/1‑150 (MQ=255)
                                                                                                                  ttCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTc                                  <  1:694455/150‑1 (MQ=255)
                                                                                                                    cccGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCCTCCCTTC‑‑ATCACCTGTATAGGTAATAACAATAATGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTAATCTAATACTTCat                                >  1:775212/1‑149 (MQ=255)
                                                                                                                    cccGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAg                                >  2:1584242/1‑150 (MQ=255)
                                                                                                                     ccGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGa                               >  2:654661/1‑150 (MQ=255)
                                                                                                                       gAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCAAGCCAAAGTAACCCATCTTCATCTACTACTTCAGaaa                             >  2:1091525/1‑150 (MQ=255)
                                                                                                                         aGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGaaaaa                           <  1:157858/150‑1 (MQ=255)
                                                                                                                         aGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGaaaaa                           <  1:1096789/150‑1 (MQ=255)
                                                                                                                          gTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAAc                          <  1:540299/150‑1 (MQ=255)
                                                                                                                           ttATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCACTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACt                         >  1:718068/1‑150 (MQ=255)
                                                                                                                             aTTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTAc                       <  2:659175/150‑1 (MQ=255)
                                                                                                                             aTTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTAc                       >  2:165749/1‑150 (MQ=255)
                                                                                                                                     tttGCAGGGATGACTCTTGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAg               <  1:1907764/150‑1 (MQ=255)
                                                                                                                                      ttGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCACTGTATAGGTTATAACAATACTGATACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAgg              <  2:1203730/150‑1 (MQ=255)
                                                                                                                                       tGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCAGCCAGATTCTGTGTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGGTTCCAGGt             >  2:891696/1‑150 (MQ=255)
                                                                                                                                        gaaGGGATGACCAGAGTCATCCCTGCAAATCCCATGTATAGGTATTAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACACATCTTCATCTACTACTTCAGAAATACTACTTTTCCAgg              >  1:69322/3‑150 (MQ=255)
                                                                                                                                         cAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTCCCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGa           >  1:1147400/1‑150 (MQ=255)
                                                                                                                                          aGGGATGACCAGAGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGaa          >  1:1732536/1‑150 (MQ=255)
                                                                                                                                           gggATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAg         <  1:1459068/150‑1 (MQ=255)
                                                                                                                                           gggATGACTCTGGTCATCCCGTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAg         <  1:905494/150‑1 (MQ=255)
                                                                                                                                            ggATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAgg        <  2:1200215/150‑1 (MQ=255)
                                                                                                                                             gATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAggg       <  2:461986/150‑1 (MQ=255)
                                                                                                                                             gATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAggg       <  2:1355253/150‑1 (MQ=255)
                                                                                                                                             gATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATAATGTTACCCTGCCAGACCCCTCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCATCATCTACTACTTCAGCAAAACTACGATTCCAGGTGAAggg       >  1:1331678/1‑150 (MQ=255)
                                                                                                                                               tGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTa     <  1:194046/150‑1 (MQ=255)
                                                                                                                                               tGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTa     <  1:165922/150‑1 (MQ=255)
                                                                                                                                                gACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTAc    <  2:1743788/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACAT  >  NC_000913/2884430‑2884725

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: