Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,884,577 T→C 6.0% 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 6.0% 562.9 / 16.7 182R881R (AGA→AGG) ‡cas3CRISPR‑associated endonuclease/helicase Cas3
Reads supporting (aligned to +/- strand):  ref base T (77/94);  new base C (4/7);  total (81/101)
Fisher's exact test for biased strand distribution p-value = 7.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 189 total aligned reads displayed.

ACTGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACAT  >  NC_000913/2884428‑2884725
                                                                                                                                                     |                                                                                                                                                    
aCTGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGAct                                                                                                                                                      >  2:826206/1‑150 (MQ=255)
       atatAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCa                                                                                                                                               <  2:785289/150‑1 (MQ=255)
            aaaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCt                                                                                                                                          >  1:572634/1‑150 (MQ=255)
             aaTCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCtt                                                                                                                                         <  1:897572/150‑1 (MQ=255)
               tctcATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCTGCa                                                                                                                                       <  1:1840497/150‑1 (MQ=255)
               tctcATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACCAGAGTCATCCCTGCa                                                                                                                                       <  1:1227444/150‑1 (MQ=255)
                    tatcGGTAAATTAAAATATTATGTACATTGTTCACCTTCTCTACTTAAGTAGGGATAAACCTTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATcccc                                                                                                                                  <  2:845488/147‑1 (MQ=255)
                    tACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATcccc                                                                                                                                  >  2:137194/1‑150 (MQ=255)
                      ccGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTg                                                                                                                                <  2:1987908/150‑1 (MQ=255)
                      ccGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTg                                                                                                                                <  2:268300/150‑1 (MQ=255)
                        gggAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGta                                                                                                                              >  2:1214452/1‑150 (MQ=255)
                          gAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGtata                                                                                                                            <  1:547350/150‑1 (MQ=255)
                              ttAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGta                                                                                                                        <  1:836972/150‑1 (MQ=255)
                              ttAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGta                                                                                                                        >  2:389158/1‑150 (MQ=255)
                                aaaaGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGtaat                                                                                                                      >  2:1058433/1‑150 (MQ=255)
                                   agaagaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGtata                                                                                                                            <  1:697870/141‑1 (MQ=255)
                                     aagaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCGTCCCTTCATCCCCTGTATAGGTAATAACaa                                                                                                                 >  2:1524486/1‑150 (MQ=255)
                                       gaTGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATa                                                                                                               >  1:1122472/1‑150 (MQ=255)
                                         tGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACt                                                                                                             >  1:138288/1‑150 (MQ=255)
                                         tGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACt                                                                                                             >  1:781126/1‑150 (MQ=255)
                                          gTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTg                                                                                                            >  2:530655/1‑150 (MQ=255)
                                          gTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTg                                                                                                            <  1:1585638/150‑1 (MQ=255)
                                           tACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGt                                                                                                           >  1:39055/1‑150 (MQ=255)
                                             cATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTa                                                                                                         >  2:1270451/1‑150 (MQ=255)
                                              aTTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTAc                                                                                                        <  1:1505007/150‑1 (MQ=255)
                                                tgtgCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTAccc                                                                                                      <  2:179458/150‑1 (MQ=255)
                                                  tgCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTg                                                                                                    <  1:1371199/150‑1 (MQ=255)
                                                  tgCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTg                                                                                                    >  2:1666259/1‑150 (MQ=255)
                                                   ggtttatCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTAACGAAGATATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGc                                                                                                   >  2:824778/7‑150 (MQ=255)
                                                    gtttatCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGcc                                                                                                  >  2:777378/6‑150 (MQ=255)
                                                    cACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGcc                                                                                                  >  2:433984/1‑150 (MQ=255)
                                                        ttCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAc                                                                                              <  1:604086/150‑1 (MQ=255)
                                                        atCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTATCCTGCCAgcc                                                                                              >  2:848602/2‑148 (MQ=255)
                                                        atCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAc                                                                                              >  1:320376/2‑150 (MQ=255)
                                                        atCCCTACTTAAGTAGCGATAAAAATTTATTGGTCTTATTATCGAAATAGATAACAAACATTCAAGAAGTTATTTGCGATTTGCAGGGATGACTCTGGTCATACATTCATCCCATGTATAGGTAATAACAATACTGTTACCCTGCCAGAc                                                                                              >  2:332282/2‑150 (MQ=255)
                                                         tCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAcc                                                                                             <  1:956936/150‑1 (MQ=255)
                                                          cccTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAccc                                                                                            >  2:629374/1‑150 (MQ=255)
                                                          cccTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGAccc                                                                                            >  2:1403330/1‑150 (MQ=255)
                                                           ccTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACcca                                                                                           >  1:1693888/1‑150 (MQ=255)
                                                           ccTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACcca                                                                                           >  2:1013956/1‑150 (MQ=255)
                                                            cTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccat                                                                                          <  2:299631/150‑1 (MQ=255)
                                                            cTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccat                                                                                          >  2:1321132/1‑150 (MQ=255)
                                                              aCTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatcc                                                                                        <  1:1223005/150‑1 (MQ=255)
                                                              aCTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACccatcc                                                                                        >  2:495892/1‑150 (MQ=255)
                                                                 tAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCcatccatc                                                                                     >  2:405306/1‑150 (MQ=255)
                                                                 tAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATAATTCCAGAAGTTATTTGGGATTTACAGTGATGACTCTGGTCATCCCTACATCCCCTGTATAGGTAATAACAATAATGTTACACTGCCAGACCaatccatc                                                                                     >  2:1941650/1‑150 (MQ=255)
                                                                   aGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCAtccatcca                                                                                   <  1:1862369/150‑1 (MQ=255)
                                                                   aGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCAtccatcca                                                                                   >  1:241431/1‑150 (MQ=255)
                                                                   aGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCAtccatcca                                                                                   >  1:1157762/1‑150 (MQ=255)
                                                                     tAGGGATAATCCGTTATTGGTATTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGa                                                                                 <  2:1007247/150‑1 (MQ=255)
                                                                     tAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGa                                                                                 <  1:1270451/150‑1 (MQ=255)
                                                                        ggATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTc                                                                              >  2:821489/1‑150 (MQ=255)
                                                                          aTAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTg                                                                            >  1:823883/1‑150 (MQ=255)
                                                                           tAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCAGCCAGACCCATCCATCCAGATTCTGt                                                                           <  2:1480456/150‑1 (MQ=255)
                                                                            aaaCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGtt                                                                          <  2:613807/150‑1 (MQ=255)
                                                                             aaacGTTATTGGTCTTATTATCGTCATTGATATCAATGATTCCCGAAGTTATTTGTTATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGttt                                                                         <  2:711976/147‑1 (MQ=255)
                                                                             aaCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGttt                                                                         >  1:1668582/1‑150 (MQ=255)
                                                                             aaCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGttt                                                                         <  1:716021/150‑1 (MQ=255)
                                                                              aCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTc                                                                        <  1:1658826/150‑1 (MQ=255)
                                                                              aCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTc                                                                        <  1:1554859/150‑1 (MQ=255)
                                                                              aCCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTc                                                                        <  1:1362069/150‑1 (MQ=255)
                                                                               ccGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTcc                                                                       <  1:1333387/150‑1 (MQ=255)
                                                                                  ttATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCtt                                                                    >  2:1874058/1‑150 (MQ=255)
                                                                                       ggTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGc                                                               >  1:1930767/1‑150 (MQ=255)
                                                                                          cTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGt                                                                                                                         <  2:356894/89‑1 (MQ=255)
                                                                                          cTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGt                                                                                                                         >  1:356894/1‑89 (MQ=255)
                                                                                            tattatCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAg                                                          <  1:466670/150‑1 (MQ=255)
                                                                                                 tCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAAcc                                                     <  2:1563503/150‑1 (MQ=255)
                                                                                                        tGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTc                                                                   >  2:480459/1‑129 (MQ=255)
                                                                                                         gATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCa                                             >  2:1462921/1‑150 (MQ=255)
                                                                                                         gATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCa                                             <  1:495452/150‑1 (MQ=255)
                                                                                                         gATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCa                                             <  1:777378/150‑1 (MQ=255)
                                                                                                         gATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCa                                             <  1:137194/150‑1 (MQ=255)
                                                                                                           tAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATc                                           <  1:629374/150‑1 (MQ=255)
                                                                                                           tAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATc                                           >  1:153980/1‑150 (MQ=255)
                                                                                                            aaCAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATct                                          <  1:1536474/150‑1 (MQ=255)
                                                                                                                  cATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTactact                                    <  1:1710076/150‑1 (MQ=255)
                                                                                                                     tCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCa                                 <  2:699090/150‑1 (MQ=255)
                                                                                                                      cccGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAg                                >  2:470279/1‑150 (MQ=255)
                                                                                                                        cGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGaa                              >  1:1021313/1‑150 (MQ=255)
                                                                                                                            gTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAAc                          <  2:930945/150‑1 (MQ=255)
                                                                                                                               aTTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTAc                       >  1:1413860/1‑150 (MQ=255)
                                                                                                                               aTTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTAc                       <  1:405306/150‑1 (MQ=255)
                                                                                                                                 ttGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGt                     >  1:301150/1‑150 (MQ=255)
                                                                                                                                  tGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGtt                    <  1:1739087/150‑1 (MQ=255)
                                                                                                                                     gATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTcc                 <  1:1013956/150‑1 (MQ=255)
                                                                                                                                     gATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTcc                 <  1:495895/150‑1 (MQ=255)
                                                                                                                                        ttGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACcca                                                                                           <  2:418281/73‑1 (MQ=255)
                                                                                                                                        ttGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACcca                                                                                           >  1:418281/1‑73 (MQ=255)
                                                                                                                                        ttGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAgg              <  1:704718/150‑1 (MQ=255)
                                                                                                                                        ttGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAgg              <  1:208059/150‑1 (MQ=255)
                                                                                                                                        ttGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAgg              <  1:1649840/150‑1 (MQ=255)
                                                                                                                                         tGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACTCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCa                                 >  2:1840426/1‑130 (MQ=255)
                                                                                                                                         tGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGt             >  1:195001/1‑150 (MQ=255)
                                                                                                                                         tGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGt             >  1:187996/1‑150 (MQ=255)
                                                                                                                                          gCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTg            <  2:714601/150‑1 (MQ=255)
                                                                                                                                               gATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTTAAggg       <  1:258732/150‑1 (MQ=255)
                                                                                                                                               gATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAACCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAggg       <  2:1599809/150‑1 (MQ=255)
                                                                                                                                                    ctctGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACAt  >  2:103094/1‑150 (MQ=255)
                                                                                                                                                    ctctGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACAt  <  2:320376/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ACTGTGAATATAAAATCTCATACCGGGAAATTAAAAGAAGATGTACATTGTGCACCTTCCCTACTTAAGTAGGGATAAACCGTTATTGGTCTTATTATCGTCATTGATAACAATCATTCCCGAAGTTATTTGGGATTTGCAGGGATGACTCTGGTCATCCCTTCATCCCCTGTATAGGTAATAACAATACTGTTACCCTGCCAGACCCATCCATCCAGATTCTGTTTCCCTTCCAGCCAAAGTAACCCATCTTCATCTACTACTTCAGAAAAACTACGTTTCCAGGTGAAGGGTACAT  >  NC_000913/2884428‑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: