Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,884,580 G→A 5.4% T880T (ACC→ACT 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,5800GA 5.4% 538.6 / 11.4 166T880T (ACC→ACTcas3CRISPR‑associated endonuclease/helicase Cas3
Reads supporting (aligned to +/- strand):  ref base G (72/85);  new base A (3/6);  total (75/91)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

Only 100 of 175 total aligned reads displayed.

GTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTG  >  NC_000913/2884431‑2884727
                                                                                                                                                     |                                                                                                                                                     
gtgAATATAAAATATAATACAGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACAGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTg                                                                                                                                                       >  2:276767/1‑150 (MQ=255)
 tgAATATTAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTgg                                                                                                                                                      <  1:1825498/150‑1 (MQ=255)
 tgAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTgg                                                                                                                                                      <  1:242200/150‑1 (MQ=255)
 tgAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTgg                                                                                                                                                      <  2:532103/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)
       taAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATcc                                                                                                                                                <  1:1813419/150‑1 (MQ=255)
        aaaaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATccc                                                                                                                                               >  2:1844247/1‑150 (MQ=255)
         aaaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTATGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCt                                                                                                                                              <  1:2447154/150‑1 (MQ=255)
         aaaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCt                                                                                                                                              <  1:1758364/150‑1 (MQ=255)
         aaaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCt                                                                                                                                              <  1:108384/150‑1 (MQ=255)
         aaaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCt                                                                                                                                              <  1:491231/150‑1 (MQ=255)
          aaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCtt                                                                                                                                             <  2:541347/150‑1 (MQ=255)
           aTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTc                                                                                                                                            <  1:479382/150‑1 (MQ=255)
           aTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCtgc                                                                                                                                            <  1:2245692/150‑3 (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)
                   ccGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTg                                                                                                                                  >  1:688436/1‑150 (MQ=255)
                    cGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGt                                                                                                                                 >  1:256721/1‑150 (MQ=255)
                      ggAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGtat                                                                                                                               <  2:434303/150‑1 (MQ=255)
                         aaTTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAgg                                                                                                                            <  2:423307/150‑1 (MQ=255)
                               aagaagATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGtaataa                                                                                                                      <  2:941139/150‑1 (MQ=255)
                               aagaagATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGt                                                                                                                           >  2:2067894/1‑145 (MQ=255)
                               aagaagATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGt                                                                                                                           <  1:2067894/145‑1 (MQ=255)
                                agaagaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAAc                                                                                                                     <  2:834241/150‑1 (MQ=255)
                                 gaagaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACa                                                                                                                    >  1:813984/1‑150 (MQ=255)
                                    gaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATa                                                                                                                 <  1:2393668/150‑1 (MQ=255)
                                    gaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATa                                                                                                                 >  2:1079728/1‑150 (MQ=255)
                                      tGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACt                                                                                                               >  2:1959906/1‑150 (MQ=255)
                                       gTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTg                                                                                                              <  2:49312/150‑1 (MQ=255)
                                       gTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTg                                                                                                              <  1:1064386/150‑1 (MQ=255)
                                        tACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCACGAAGTTATTAGGGATTTGCAGGGATGACTCTGGTAATCCCTTC‑‑ATCCCCTGTATAGGTAATAAAAATACTGt                                                                                                             >  1:2223379/1‑150 (MQ=255)
                                         aCATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGtt                                                                                                            >  1:1151853/1‑150 (MQ=255)
                                           aTTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTAc                                                                                                          <  2:1220589/150‑1 (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                                                                                                    <  1:373372/150‑1 (MQ=255)
                                                  tttatCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCa                                                                                                   >  1:1935753/5‑150 (MQ=255)
                                                     ttCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAc                                                                                                <  1:2468822/150‑1 (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                                                                                     <  2:692556/150‑1 (MQ=255)
                                                                 gTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAg                                                                                    <  2:2253500/150‑1 (MQ=255)
                                                                  tAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGa                                                                                   <  1:1654961/150‑1 (MQ=255)
                                                                    gggATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGAtt                                                                                 >  1:2407377/1‑150 (MQ=255)
                                                                     ggATAAACCGTTATTGTTCTTATTATCGTCATTGATAACAATCTTTCACGAAGTTATTTGGGATTTGCATGGATGAATCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAAAAATACTGTTACCCTGCCAGACCCATCCATCCAGATTc                                                                                <  1:1863434/150‑1 (MQ=255)
                                                                       aTAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTg                                                                              >  1:578086/1‑150 (MQ=255)
                                                                          aaCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCCTCCAGATTCTGttt                                                                           >  2:1120296/1‑150 (MQ=255)
                                                                           aCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTc                                                                          <  2:1151853/150‑1 (MQ=255)
                                                                            ccGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTcc                                                                         <  2:761536/150‑1 (MQ=255)
                                                                              gTTATTGGTCTTATTATCGTCATTTATAACAATCATTCCCGAAGTTATTTGGGATTTGCATTGATGACTCTGGTCATCTGTTT‑‑TTTTTATGTATAGGTATTAACAATACTTTTACACTGCAATACACATTTTTCCATATTGTGTTTcttt                                                                       <  2:1065231/150‑4 (MQ=255)
                                                                                tATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTc                                                                     <  2:2302321/150‑1 (MQ=255)
                                                                                tATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTc                                                                     >  2:2297894/1‑150 (MQ=255)
                                                                                 aTTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTcc                                                                    <  2:2126042/150‑1 (MQ=255)
                                                                                  ttGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCa                                                                   >  1:136507/1‑150 (MQ=255)
                                                                                    ggTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGc                                                                 <  2:1029932/150‑1 (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                                                          <  1:2287505/150‑1 (MQ=255)
                                                                                            tatCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTaa                                                         >  2:63349/1‑150 (MQ=255)
                                                                                            tatCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTaa                                                         <  2:153230/150‑1 (MQ=255)
                                                                                            tatCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTaa                                                         >  2:786359/1‑150 (MQ=255)
                                                                                               cGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAAccc                                                      >  2:2026977/1‑150 (MQ=255)
                                                                                                 tCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCAt                                                    >  1:898692/1‑150 (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)
                                                                                                              tCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCtactac                                       >  1:1053684/1‑150 (MQ=255)
                                                                                                                  tCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCa                                   <  1:2119142/150‑1 (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:555561/1‑150 (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)
                                                                                                                                    tttGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAg                 <  2:136507/150‑1 (MQ=255)
                                                                                                                                    tttGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAg                 <  1:435048/150‑1 (MQ=255)
                                                                                                                                      tgaaGGGATGACCAGAGTCATCCCTGCAAATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAg                 >  2:1658625/4‑150 (MQ=255)
                                                                                                                                       gCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTg              >  2:1692603/1‑150 (MQ=255)
                                                                                                                                        cAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGCCCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGa             >  2:251653/1‑150 (MQ=255)
                                                                                                                                        aaGGGATGACCAGAGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACAATGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGGGa             >  2:732651/2‑150 (MQ=255)
                                                                                                                                         aGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGaa            <  1:671270/150‑1 (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)
                                                                                                                                              tGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTa       <  1:606169/150‑1 (MQ=255)
                                                                                                                                               gACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTAc      <  2:301982/150‑1 (MQ=255)
                                                                                                                                                aCTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAGAACTACGTTTCCAGGTGAAGGGTACa     >  2:613003/1‑150 (MQ=255)
                                                                                                                                                aCTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACa     >  1:2330008/1‑150 (MQ=255)
                                                                                                                                                   ctGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACACTGCCAGCCACATCCATCCAGATTATGTTTCCCTTCCACCCAAAGTAACCAATCTTCATCTACTACTTCAGAAAAACAACGTTTCCACCTCAAGGGTACATTg  >  1:2306431/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTC‑‑ATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACATTG  >  NC_000913/2884431‑2884727

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: