Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,884,577 T→C 5.2% R881R (AGA→AGG) ‡ cas3 ← CRISPR‑associated endonuclease/helicase Cas3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,884,5770TC 5.2% 546.4 / 12.7 173R881R (AGA→AGG) ‡cas3CRISPR‑associated endonuclease/helicase Cas3
Reads supporting (aligned to +/- strand):  ref base T (72/92);  new base C (3/6);  total (75/98)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

Only 100 of 175 total aligned reads displayed.

CTGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACAT  >  NC_000913/2884429‑2884725
                                                                                                                                                    |                                                                                                                                                      
  gtgAATATAAAATATAATACAGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACAGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTg                                                                                                                                                     >  2:276767/1‑150 (MQ=255)
   tgAATATTAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTgg                                                                                                                                                    <  1:1825498/150‑1 (MQ=255)
   tgAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCGTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTgg                                                                                                                                                    <  1:2311905/150‑1 (MQ=255)
   tgAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTgg                                                                                                                                                    <  1:233644/150‑1 (MQ=255)
   tgAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTgg                                                                                                                                                    <  2:532103/150‑1 (MQ=255)
   tgAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTgg                                                                                                                                                    <  1:242200/150‑1 (MQ=255)
   tgAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTgg                                                                                                                                                    >  2:2461539/1‑150 (MQ=255)
     aaTATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTc                                                                                                                                                  >  1:2303457/1‑150 (MQ=255)
     aaTATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTc                                                                                                                                                  <  1:1209962/150‑1 (MQ=255)
        ataAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATc                                                                                                                                               >  1:430000/1‑150 (MQ=255)
         taAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATcc                                                                                                                                              <  1:1813419/150‑1 (MQ=255)
          aaaaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATccc                                                                                                                                             >  2:1844247/1‑150 (MQ=255)
           aaaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCt                                                                                                                                            <  1:1758364/150‑1 (MQ=255)
           aaaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCt                                                                                                                                            <  1:491231/150‑1 (MQ=255)
           aaaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCt                                                                                                                                            <  1:108384/150‑1 (MQ=255)
            aaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCtt                                                                                                                                           <  2:541347/150‑1 (MQ=255)
            aaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCtg                                                                                                                                           <  1:537626/150‑2 (MQ=255)
             aTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCtgc                                                                                                                                          <  1:2245692/150‑3 (MQ=255)
              tctcATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTTCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑a                                                                                                                                       >  1:301982/1‑150 (MQ=255)
                  aTACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATccc                                                                                                                                   <  2:345783/150‑1 (MQ=255)
                   tACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATcccc                                                                                                                                  <  1:276421/150‑1 (MQ=255)
                     ccGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTg                                                                                                                                <  2:905757/150‑1 (MQ=255)
                      cGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGt                                                                                                                               <  2:1142601/150‑1 (MQ=255)
                        ggAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGtat                                                                                                                             <  2:434303/150‑1 (MQ=255)
                        ggAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGtat                                                                                                                             <  1:652922/150‑1 (MQ=255)
                           aaTTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAgg                                                                                                                          <  2:423307/150‑1 (MQ=255)
                              tAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGtaa                                                                                                                       >  2:1410807/1‑150 (MQ=255)
                                 aagaagATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGtaataa                                                                                                                    <  2:941139/150‑1 (MQ=255)
                                 aagaagATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGt                                                                                                                         <  1:2067894/145‑1 (MQ=255)
                                 aagaagATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGt                                                                                                                         >  2:2067894/1‑145 (MQ=255)
                                  agaagaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAAc                                                                                                                   <  2:834241/150‑1 (MQ=255)
                                      gaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATa                                                                                                               >  2:1079728/1‑150 (MQ=255)
                                        tGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACt                                                                                                             >  2:1959906/1‑150 (MQ=255)
                                         gTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTg                                                                                                            <  1:1064386/150‑1 (MQ=255)
                                         gTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTg                                                                                                            <  2:49312/150‑1 (MQ=255)
                                          tACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCACGAAGTTATTAGGGATTTGCAGGGATGACTCTGGTAATCCCTTC‑‑ATCCCCTGTATAGGTAATAAAAATACTGt                                                                                                           >  1:2223379/1‑150 (MQ=255)
                                           aCATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGtt                                                                                                          >  1:1151853/1‑150 (MQ=255)
                                              ttGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTAcc                                                                                                       >  1:1315324/1‑150 (MQ=255)
                                                gtgCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCt                                                                                                     >  1:1868363/1‑150 (MQ=255)
                                                 tgCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTg                                                                                                    >  1:2253500/1‑150 (MQ=255)
                                                  gCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGc                                                                                                   >  2:1572244/1‑150 (MQ=255)
                                                   cACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGcc                                                                                                  >  2:124776/1‑150 (MQ=255)
                                                    tttatCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCa                                                                                                 >  1:1935753/5‑150 (MQ=255)
                                                            tACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGTATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatc                                                                                         <  1:1654987/150‑1 (MQ=255)
                                                              cTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatcca                                                                                       >  2:1519737/1‑150 (MQ=255)
                                                              cTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatcca                                                                                       >  2:419710/1‑150 (MQ=255)
                                                                 aaGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCatccatcc                                                                                    <  1:1140277/150‑1 (MQ=255)
                                                                  aGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCAtccatcca                                                                                   <  1:2422212/150‑1 (MQ=255)
                                                                  aGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCAtccatcca                                                                                   <  2:692556/150‑1 (MQ=255)
                                                                   gTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAg                                                                                  <  2:2253500/150‑1 (MQ=255)
                                                                   gTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAg                                                                                  <  2:1101874/150‑1 (MQ=255)
                                                                    tAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGa                                                                                 <  1:1654961/150‑1 (MQ=255)
                                                                     aGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGAt                                                                                <  1:2460893/150‑1 (MQ=255)
                                                                       ggATAAACCGTTATTGTTCTTATTATCGTCATTGATAACAATCTTTCACGAAGTTATTTGGGATTTGCATGGATGAATCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAAAAATACTGTTACCCTGCCAGACCCATCCATCCAGATTc                                                                              <  1:1863434/150‑1 (MQ=255)
                                                                            aaCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGttt                                                                         >  1:885148/1‑150 (MQ=255)
                                                                              ccGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTcc                                                                       <  2:761536/150‑1 (MQ=255)
                                                                                  tATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTc                                                                   >  2:2297894/1‑150 (MQ=255)
                                                                                  tATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTc                                                                   <  2:2302321/150‑1 (MQ=255)
                                                                                  tATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTc                                                                   <  1:68806/150‑1 (MQ=255)
                                                                                   aTTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTcc                                                                  <  2:2126042/150‑1 (MQ=255)
                                                                                    ttGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCa                                                                 >  1:136507/1‑150 (MQ=255)
                                                                                       gTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGcc                                                              >  1:815380/1‑150 (MQ=255)
                                                                                           tattatCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAg                                                          <  1:1016950/150‑1 (MQ=255)
                                                                                            attatCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGt                                                         <  1:1959906/150‑1 (MQ=255)
                                                                                             ttatCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTa                                                        <  2:1169035/150‑1 (MQ=255)
                                                                                             ttatCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTa                                                        <  1:2287505/150‑1 (MQ=255)
                                                                                              tatCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTaa                                                       <  2:153230/150‑1 (MQ=255)
                                                                                              tatCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTaa                                                       >  2:63349/1‑150 (MQ=255)
                                                                                              tatCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTaa                                                       >  2:786359/1‑150 (MQ=255)
                                                                                                 cGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAAccc                                                    >  2:2026977/1‑150 (MQ=255)
                                                                                                     aTTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCt                                                <  1:447555/150‑1 (MQ=255)
                                                                                                        gATAACAATCATTCCTGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAAACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCa                                             <  2:2185591/150‑1 (MQ=255)
                                                                                                        gATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCa                                             <  1:1712019/150‑1 (MQ=255)
                                                                                                           aaCAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATct                                          <  1:1410184/150‑1 (MQ=255)
                                                                                                             cAATCATTCCCGAGGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctac                                        >  2:704544/1‑150 (MQ=255)
                                                                                                               aTCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctacta                                      <  1:1519737/150‑1 (MQ=255)
                                                                                                               aTCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTCTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctacta                                      <  1:716678/150‑1 (MQ=255)
                                                                                                                tCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCtactac                                     >  1:1053684/1‑150 (MQ=255)
                                                                                                                    tCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCa                                 <  1:2119142/150‑1 (MQ=255)
                                                                                                                     cccGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCGCTTCCAGCCAAAGTAACCGATGTGCCTCGCCTACTTCAg                                >  1:555811/1‑150 (MQ=255)
                                                                                                                     cccGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAg                                >  2:1239297/1‑150 (MQ=255)
                                                                                                                          aGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGaaaaa                           <  1:1230681/150‑1 (MQ=255)
                                                                                                                           gTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAAc                          <  1:1127042/150‑1 (MQ=255)
                                                                                                                              aTTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTAc                       <  1:1239297/150‑1 (MQ=255)
                                                                                                                                 tGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGtt                    <  1:1410807/150‑1 (MQ=255)
                                                                                                                                     aTTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACACTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCa                <  2:1849141/150‑1 (MQ=255)
                                                                                                                                        tgaaGGGATGACCAGAGTCATCCCTGCAAATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAg               >  2:1658625/4‑150 (MQ=255)
                                                                                                                                         gCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATctacta                                      <  1:255566/124‑1 (MQ=255)
                                                                                                                                         gCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTg            >  2:1692603/1‑150 (MQ=255)
                                                                                                                                          cAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGCCCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGa           >  2:251653/1‑150 (MQ=255)
                                                                                                                                          aaGGGATGACCAGAGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGa           >  1:907619/2‑150 (MQ=255)
                                                                                                                                          aaGGGATGACCAGAGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACAATGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGGGa           >  2:732651/2‑150 (MQ=255)
                                                                                                                                            gggATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCACTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGATTCAAGTTGAAg         >  2:1041078/1‑150 (MQ=255)
                                                                                                                                              gATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAggg       <  1:1572244/150‑1 (MQ=255)
                                                                                                                                               aTGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGt      <  2:256721/150‑1 (MQ=255)
                                                                                                                                                tGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTa     <  1:63349/150‑1 (MQ=255)
                                                                                                                                                 gACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTAc    <  2:301982/150‑1 (MQ=255)
                                                                                                                                                 gACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTAc    <  1:600699/150‑1 (MQ=255)
                                                                                                                                                  aCTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAGAACTACGTTTCCAGGTGAAGGGTACa   >  2:613003/1‑150 (MQ=255)
                                                                                                                                                   ctctGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCAACTACTACTTCAGAAAAACAACGTTTCCAGGTGAAGGGTACAt  >  1:2398730/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CTGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACAT  >  NC_000913/2884429‑2884725

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