Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,884,580 G→A 6.0% T880T (ACC→ACT 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,5800GA 6.0% 422.2 / 10.6 134T880T (ACC→ACTygcBCascade complex anti‑viral R‑loop helicase‑annealase Cas3
Reads supporting (aligned to +/- strand):  ref base G (50/74);  new base A (2/6);  total (52/82)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

Only 100 of 139 total aligned reads displayed.

GTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGA  >  NC_000913/2884431‑2884728
                                                                                                                                                     |                                                                                                                                                      
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)
                 tACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATcccc                                                                                                                                     >  2:650562/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:1046646/150‑1 (MQ=255)
                       gAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGtata                                                                                                                               <  1:1776208/150‑1 (MQ=255)
                          aTTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGt                                                                                                                            <  2:491962/150‑1 (MQ=255)
                          aTTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGt                                                                                                                            <  2:1242596/150‑1 (MQ=255)
                             aaaaGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGtaat                                                                                                                         >  2:232284/1‑150 (MQ=255)
                                agaagaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAAc                                                                                                                      <  1:498174/150‑1 (MQ=255)
                                agaagaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAAc                                                                                                                      <  1:1106825/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:1590013/150‑1 (MQ=255)
                                                      tCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAcc                                                                                                <  1:987896/150‑1 (MQ=255)
                                                      tCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAcc                                                                                                <  2:760697/150‑1 (MQ=255)
                                                        ccTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCACGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACcca                                                                                              <  1:212391/150‑1 (MQ=255)
                                                          tACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatc                                                                                            <  1:1657799/150‑1 (MQ=255)
                                                          tACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatc                                                                                            >  2:1637596/1‑150 (MQ=255)
                                                          tACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatc                                                                                            >  1:487731/1‑150 (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)
                                                                 ctAGCGATACACCGTAATTCGTCATCTTATCGTCACTGATAACAATCAATCCCGCCGTTAAATGGGATTTGCAGCGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGATACCCTGCCAGACCCATCCATCCAg                                                                                     <  2:818561/149‑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:1113461/150‑1 (MQ=255)
                                                                           aCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTc                                                                           <  2:1539779/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:1084501/150‑1 (MQ=255)
                                                                                                gTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCa                                                      >  1:864016/1‑150 (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)
                                                                                                                 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‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACt                            <  2:678401/150‑1 (MQ=255)
                                                                                                                          ttATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCACTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACt                            >  1:718068/1‑150 (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: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:477852/150‑1 (MQ=255)
                                                                                                                                               gACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTAc       <  2:1743788/150‑1 (MQ=255)
                                                                                                                                                    tGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGa  >  2:439226/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTGA  >  NC_000913/2884431‑2884728

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