Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,880,635 T→C L138L (TTA→TTG casD ← type I‑E CRISPR system Cascade subunit CasD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,880,6350TC95.7% 134.2 / ‑4.0 47L138L (TTA→TTGcasDtype I‑E CRISPR system Cascade subunit CasD
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (22/22);  total (24/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

ACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGAT  >  NC_000913/2880392‑2880878
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
aCCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACt                                                                                                                                                                                                                                               <  1:546343/250‑1 (MQ=255)
aCCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACt                                                                                                                                                                                                                                               <  1:224340/250‑1 (MQ=255)
 cctcctTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAAATGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGAAAGGTGTATACCGAGGCTTCAATACTGc                                                                                                                                                                                                                                             >  1:322247/1‑251 (MQ=255)
    cctTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGGGGCTTCAATACTGCttt                                                                                                                                                                                                                                          >  2:142774/1‑251 (MQ=255)
       ttAAGCACATACCAGTCGCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTGCAATACTGCtttt                                                                                                                                                                                                                                         <  2:637996/249‑1 (MQ=255)
          aTCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAACAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCagg                                                                                                                                                                                                                                    >  2:149570/1‑249 (MQ=255)
              cattccATTCTCGGGAAGCAAATTGGCGAGGCAAGGTTAACATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAgggct                                                                                                                                                                                                                                <  2:14877/247‑5 (MQ=255)
                       ctcGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAAcc                                                                                                                                                                                                                       >  1:298309/1‑251 (MQ=255)
                         cggtaaGCAAATTGTCGAGGCAAGGTTATCATCGGTTCGTCGCGTGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAAATTAATAGCGCCTTCTGATGATCCGATGCCCGACATGTCCCCAAAAAAATCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAATGTGTATACCGAGGCTTCAATACTGCTTTTTCAATTTCTGAGATAACCa                                                                                                                                                                                                                      <  2:628957/246‑1 (MQ=255)
                              aGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTg                                                                                                                                                                                                                >  1:263091/1‑251 (MQ=255)
                                     tGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATgggg                                                                                                                                                                                                         >  1:38076/1‑251 (MQ=255)
                                      gTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGgt                                                                                                                                                                                                        <  1:610734/251‑1 (MQ=255)
                                             caaggtgttctgcgtttCGTTGCGTGCCTTAAATTTTAAATGATGCCCTGTAACTGATTCCACACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGAAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGCGGTGTTAAca                                                                                                                                                                                                 <  2:525280/236‑2 (MQ=255)
                                                  tGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCag                                                                                                                                                                                               >  1:16123/1‑248 (MQ=255)
                                                   ggtcatcGGTTCTTCGCGCGCTGTAAATTTTAAATGATGCCCTGTAACTAATTCCTCTCTATATATATCGCCGCCAACGGGCTAATAATTTAATAGCGCCTTCTGAGGTTCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTTTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCagaga                                                                                                                                                                                            <  2:411909/248‑1 (MQ=255)
                                                   gATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCagagag                                                                                                                                                                                           <  1:215217/251‑1 (MQ=255)
                                                    atcatcGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGc                                                                                                                                                                                          <  1:39723/251‑1 (MQ=255)
                                                       atcGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTTTATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTGCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTGTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGAc                                                                                                                                                                                       <  2:580841/251‑1 (MQ=255)
                                                               gtcgCGCGCCGTAAATTTTACATGATGCCACGTAACTGAGTCCTCACTATATATATAGACCCCAACGGGCTCATAATTTCATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCACCAAAAAAAGCGGGTGTGATAGGGGACAGATTCTCCGCCCCAGGAAAGGTGTATACCGGGGCTTCACTACTGCTTTTTCAAGTTCTAATATACCCATCGTTGTAGGAGGGGTTAACCAGAGGGCGACGGGAAagg                                                                                                                                                                               >  2:542064/1‑251 (MQ=255)
                                                               gtcgCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAaggg                                                                                                                                                                              >  2:71245/1‑250 (MQ=255)
                                                                 cgcgcgCCGTAAATTTTAAATGATGCCATGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACAGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAagg                                                                                                                                                                               >  2:34313/1‑249 (MQ=255)
                                                                   cgcgcCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGGGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGc                                                                                                                                                                           <  1:294303/251‑1 (MQ=255)
                                                                      gcCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAATGGAGGCATc                                                                                                                                                                        >  2:570838/1‑251 (MQ=255)
                                                                            aaTTTTAAATGATGCCCTTTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGATATAGCCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATaa                                                                                                                                                                  <  1:665708/251‑1 (MQ=255)
                                                                                                        ccTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTc                                                                                                                                       <  1:149570/250‑1 (MQ=255)
                                                                                                                                 gggttCATAAGTTAATAGCGCCTTCTGAGGAGCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAAtctt                                                                                                             <  2:106481/247‑1 (MQ=255)
                                                                                                                                   gcccacaaTTTACTACCGCCTTCCGAAGAGCAGATGCCTAACATGTCACCAAAAAAAGCGGGACTAGTAATGGGCAACATCTACGACCCAGGAAAGGTGTATACCGAGGCTTCCATAATGCTTTTTAAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAAACAGAGAGCGACGGAAAGGGAGGCATCCAATAAACATTCGCGCCATGTTTGAATCGTTTCAGGACTTTTAAAACCAAGGTCAtcttct                                                                                                           >  2:628652/7‑251 (MQ=255)
                                                                                                                                     tCATAATTTAATCGCGCCTTCTGAGGGTCCGATGCCTAACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGTGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTcg                                                                                                         >  2:244296/1‑251 (MQ=255)
                                                                                                                                      cATAATTTCATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGTGTAACTTCTCCGCCCCCGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGTCTGTAAAGGAGGCATCACAAAACTATTCGCGCCATGTTTGATTCGTTTCTTGACTTTTCCAACAACGGTAATCTTCCcgc                                                                                                        >  2:206590/1‑251 (MQ=255)
                                                                                                                                          aTTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTATCAAACCACGGTAATCTTCTCGCGCTc                                                                                                    >  2:467923/1‑251 (MQ=255)
                                                                                                                                                                cccGATGCCTGACATGCCCCCAAAAAAAGCTGGTGTGTTAGTGTGCTACTTCTCCGCCCCAGGGAAGGTGTATACCTAGGATTCAACACTGCTTTTTCAAGTACTGAGATAACCATCGTTGCATTGGGTGTTAAACAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCACGTTTGAATCGGTTCATTACTTTTCAAACGACGGTAATCTTCTCGCGCTCCTAGCACTGTATGTTAATCAct                                                                              <  2:76301/250‑2 (MQ=255)
                                                                                                                                                                       ccTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAAcccc                                                                       >  1:300293/1‑251 (MQ=255)
                                                                                                                                                                                           aGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATc                                                   <  1:441694/251‑1 (MQ=255)
                                                                                                                                                                                           aGCGGGTGTGGTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGTTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGGGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTCACAGACACACGACGATc                                                   <  2:298309/251‑1 (MQ=255)
                                                                                                                                                                                              gggTGTGTTAGTGGGCGTCTTCACCGCCCAAGGTAAAGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAGCCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGta                                                <  2:12440/251‑2 (MQ=255)
                                                                                                                                                                                               ggTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTc                                                >  1:75507/1‑250 (MQ=255)
                                                                                                                                                                                                  gtgtTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGa                                            <  1:467923/251‑1 (MQ=255)
                                                                                                                                                                                                             caccTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACAACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATACCGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTcctc                                 >  2:531629/4‑248 (MQ=255)
                                                                                                                                                                                                                     ccGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAAtct                                                                                                              <  2:14481/166‑1 (MQ=255)
                                                                                                                                                                                                                     ccGCCCCAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAAtct                                                                                                              >  1:14481/1‑166 (MQ=255)
                                                                                                                                                                                                                        acccAGGAAAGATGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATAGTTGCATGGGGTGTTAACCAGATAGCGACGGTAAAAGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGATAATCTTCTCTCGCTCCAAGGACTGTATGGTAATCACGAAACCCCGTTACAGACACACGACAATCGTCAAGATTGAGGTCATCGCAGCGCaatg                      >  2:74709/2‑248 (MQ=255)
                                                                                                                                                                                                                         cccAGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTGTCATGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACAAGAGAGCGCCGGTAAAGAGGGCATCACATAAATATTCGCGCCATGTTTACATCGTTTCATGACTTTTCAAACCACGATAAGATTCTCGCTCTGAAAGGACTGTATGGTAATAAATCAGCCCCGTTACAGACACACGACGATCTCCAAGAATAAGCTATTCGCAGCCCACTGc                     >  2:21073/1‑251 (MQ=255)
                                                                                                                                                                                                                            aGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCa                                                                                                                            >  1:96217/1‑145 (MQ=255)
                                                                                                                                                                                                                            aGGTAAGGTGTATACCGAGGCTTCAATACTGCTTTTTCAAGGTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCa                                                                                                                            <  2:96217/145‑1 (MQ=255)
                                                                                                                                                                                                                                             aGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGTTAATCACGCAACCCCGTTa                                                                   <  2:385542/185‑1 (MQ=255)
                                                                                                                                                                                                                                             aGGCTTCAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTa                                                                   >  1:385542/1‑185 (MQ=255)
                                                                                                                                                                                                                                             aGGCTTAAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGAATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGAAGCGCACTGCAAATTGACCACTCTCTGAt  >  2:513146/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
ACCTCCTTTAATCACATACCATTCTCGGGAAGCAAATTGTCGAGGCAAGGTGATCATCGGTTCGTCGCGCGCCGTAAATTTTAAATGATGCCCTGTAACTGATTCCTCACTATATATATCGCCGCCAACGGGCTCATAATTTAATAGCGCCTTCTGAGGATCCGATGCCTGACATGTCCCCAAAAAAAGCGGGTGTGTTAGTGGGCAACTTCTCCGCCCCAGGTAAGGTGTATACCGAGGCTTTAATACTGCTTTTTCAAGTTCTGAGATAACCATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGAT  >  NC_000913/2880392‑2880878

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