Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,880,912 A→G I46T (ATC→ACC)  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,9120AG96.2% 255.8 / 0.6 81I46T (ATC→ACC) casDtype I‑E CRISPR system Cascade subunit CasD
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (40/37);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

CATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGC  >  NC_000913/2880665‑2881158
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
cATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTcc                                                                                                                                                                                                                                                     >  1:149719/1‑251 (MQ=255)
 aTCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTccc                                                                                                                                                                                                                                                    <  2:167724/251‑1 (MQ=255)
    gTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCa                                                                                                                                                                                                                                                   >  1:506290/1‑249 (MQ=255)
             ggTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTGGGGTCCTAAGACAAGcccc                                                                                                                                                                                                                                        <  1:484819/251‑1 (MQ=255)
                              gACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAAcc                                                                                                                                                                                                                       >  2:656434/1‑251 (MQ=255)
                              gACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAAcc                                                                                                                                                                                                                       >  1:234061/1‑251 (MQ=255)
                               aCGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAAccc                                                                                                                                                                                                                      >  1:327647/1‑251 (MQ=255)
                                cGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCg                                                                                                                                                                                                                     <  2:156709/251‑1 (MQ=255)
                                         aggCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTc                                                                                                                                                                                                            <  2:478453/251‑1 (MQ=255)
                                          ggCATCACATCAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAAGAACCCGCTTCGGGTCg                                                                                                                                                                                                           <  1:225378/251‑1 (MQ=255)
                                            cATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCgg                                                                                                                                                                                                          >  2:113697/1‑250 (MQ=255)
                                             aTCACATAAATATTCGCGCCATGTTTGAGTCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGa                                                                                                                                                                                                         >  2:4370/1‑250 (MQ=255)
                                             aTCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTACTCACGAAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGa                                                                                                                                                                                                         >  2:12890/1‑250 (MQ=255)
                                                  aTAAATATTCGCACCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCtt                                                                                                                                                                                                   >  1:488445/1‑251 (MQ=255)
                                                              gccaagcttgcaTCGTTTCATGACTTTTCAAACCACGGCAATCTTCCCGCGCTCCAAGGACTATATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAACTTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTcgc                                                                                                                                                                                       >  2:296704/12‑251 (MQ=255)
                                                               ccATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTcgcg                                                                                                                                                                                      >  1:602312/1‑251 (MQ=255)
                                                               ccATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTcgcg                                                                                                                                                                                      >  1:538170/1‑251 (MQ=255)
                                                               ccATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTcgcg                                                                                                                                                                                      >  2:465811/1‑251 (MQ=255)
                                                                   gTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTcc                                                                                                                                                                                  <  1:564312/251‑1 (MQ=255)
                                                                         aTCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCaa                                                                                                                                                                             >  2:71061/1‑250 (MQ=255)
                                                                            gTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGt                                                                                                                                                                         >  1:215253/1‑251 (MQ=255)
                                                                             tttCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTc                                                                                                                                                                        >  2:560210/1‑251 (MQ=255)
                                                                              ttCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTc                                                                                                                                                                        <  1:577119/250‑1 (MQ=255)
                                                                                cATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCgg                                                                                                                                                                      <  2:602312/250‑1 (MQ=255)
                                                                                   gACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAg                                                                                                                                                                                                                                     >  1:438167/1‑184 (MQ=255)
                                                                                   gACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAg                                                                                                                                                                                                                                     <  2:438167/184‑1 (MQ=255)
                                                                                   gACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGc                                                                                                                                                                  >  2:688628/1‑251 (MQ=255)
                                                                                      ttttCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCtaat                                                                                                                                                                                                                           >  2:370924/1‑191 (MQ=255)
                                                                                      ttttCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCtaat                                                                                                                                                                                                                           <  1:370924/191‑1 (MQ=255)
                                                                                         tCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGccccc                                                                                                                                                              >  1:639864/1‑249 (MQ=255)
                                                                                           aaaCCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGc                                                                                                                                                          >  1:137704/1‑251 (MQ=255)
                                                                                             aCCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCCCGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCtt                                                                                                                                                        <  1:606047/251‑1 (MQ=255)
                                                                                                 cGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCAt                                                                                                                                                    <  1:184803/251‑1 (MQ=255)
                                                                                                   gTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTg                                                                                                                                                  >  2:566856/1‑251 (MQ=255)
                                                                                                     aaTCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGc                                                                                                                                                >  2:135495/1‑251 (MQ=255)
                                                                                                          tctcGCGCGCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGccc                                                                                                                                              <  2:599803/248‑1 (MQ=255)
                                                                                                             cgcgCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAg                                                                                                                                        <  1:465811/251‑1 (MQ=255)
                                                                                                                gcTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCTTTGGCCCAGCAAGCCg                                                                                                                                     >  2:585896/1‑251 (MQ=255)
                                                                                                                        gACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATc                                                                                                                              >  2:651202/1‑250 (MQ=255)
                                                                                                                        gACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCa                                                                                                                             <  1:71061/251‑1 (MQ=255)
                                                                                                                            gTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCaaa                                                                                                                           <  1:173309/249‑1 (MQ=255)
                                                                                                                              aTGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAAAGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGATCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATa                                                                                                                         >  2:365108/1‑249 (MQ=255)
                                                                                                                                    aTCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGAtctc                                                                                                                  >  1:522122/1‑250 (MQ=255)
                                                                                                                                      cACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATg                                                                                                               <  1:688628/251‑1 (MQ=255)
                                                                                                                                       aCGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGt                                                                                                              >  1:637290/1‑251 (MQ=255)
                                                                                                                                        cGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGtt                                                                                                             <  1:678670/251‑1 (MQ=255)
                                                                                                                                           aaCCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATTTCATGTTCAc                                                                                                          <  1:135495/251‑1 (MQ=255)
                                                                                                                                                      aGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCg                                                                                                                                                                                                                <  1:387496/138‑1 (MQ=255)
                                                                                                                                                      aGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCg                                                                                                                                                                                                                >  2:387496/1‑138 (MQ=255)
                                                                                                                                                       gACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCatt                                                                                               >  1:646931/1‑250 (MQ=255)
                                                                                                                                                                   tCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACCCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGTTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCa                                                                                   <  1:631115/250‑1 (MQ=255)
                                                                                                                                                                       aaaGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCCCTGATAACGCCTGTAATGAAGAAGTATCACCACGTTGGGTCCCAAGACAAGCCCCGGGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGAtt                                                                              <  2:639864/250‑1 (MQ=255)
                                                                                                                                                                           aaTGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTa                                                                           >  2:647154/1‑250 (MQ=255)
                                                                                                                                                                           aaTGAGCTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTa                                                                           >  2:573764/1‑250 (MQ=255)
                                                                                                                                                                                gTTCATCGCAGCGCACTGCAAATTGCACCCTCTCGGAAAACGCCTGTAATGAAGAAGTATCATCACGTGGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGtt                                                                     <  2:637290/251‑1 (MQ=255)
                                                                                                                                                                                gTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGAtctc                                                                                                                  <  2:641671/206‑1 (MQ=255)
                                                                                                                                                                                gTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGAtctc                                                                                                                  >  1:641671/1‑206 (MQ=255)
                                                                                                                                                                                gTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGtt                                                                     >  2:441555/1‑251 (MQ=255)
                                                                                                                                                                                gTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTGCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGtt                                                                     <  2:215253/251‑1 (MQ=255)
                                                                                                                                                                                 ttCATCGCAGCGTACTGCAAATTGCACGCTCTCTGATAACCCCTTTAATGAAAAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTGGGAAATCTTCCGGAAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGTTTGCATTGGCACAGCAAGCCGCAAGATAAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGtt                                                                     <  1:104258/250‑1 (MQ=255)
                                                                                                                                                                                  tCATCGCAGCGCACTGCAAATTGCACGCTCGCTGATAACCCCTTGAATGAAGAAGTTTCATCGCTTTGGGTCCCAGGACAAGCCCCGATTGGCCCTAATACCCCGCTTCGGGTCGAAAATCTTCCGGTAGATCGCTTTCCTTCAAAGGTCGGCTGCCCCGAGGCTTGCATTGGCCCGGCAAGCCGCAATATCAAATAAGACCTCGTGTTCCCGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTc                                                                    <  1:296704/250‑1 (MQ=255)
                                                                                                                                                                                        cAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCat                                                                                                >  2:44743/1‑216 (MQ=255)
                                                                                                                                                                                        cAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCat                                                                                                <  1:44743/216‑1 (MQ=255)
                                                                                                                                                                                         aGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGt                                                              >  1:385607/1‑249 (MQ=255)
                                                                                                                                                                                           cgcACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGATCGGTCGCGTTCCTTCAAAGGTCGGCTGACCCCAGGCTTGCATTGACAAAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTCAAGTAGGc                                                          >  1:277805/1‑251 (MQ=255)
                                                                                                                                                                                                  cAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGtt                                                   <  1:4370/251‑1 (MQ=255)
                                                                                                                                                                                                    aatggCACACCCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCTTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATTTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTa                                                 <  2:488445/247‑1 (MQ=255)
                                                                                                                                                                                                     aTTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTa                                                 >  1:337069/1‑250 (MQ=255)
                                                                                                                                                                                                          acactatctgaaaaggaatGTAATGAAGAAGTATCATCACGTTGGGCCCCATGACAAGCCCCGAGTAGCCCTAATAACCAGCTTCGGTTCGGAAATCTTCCGGCATGTCGCGTTCCTTCAACGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCACGCCGCACGCTGAAATAAGATCTCATGTTCACGCCTCGCAATTATTACGACACCAGGATTTCACCTGTTCTAAAGGAGGCATTGGTTTacgttg                                            >  2:277390/19‑245 (MQ=255)
                                                                                                                                                                                                            cctctctGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGc                                         <  1:566856/250‑1 (MQ=255)
                                                                                                                                                                                                            cctctctGAAAACGCCTGTAATGAAAAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGc                                         <  2:506290/250‑1 (MQ=255)
                                                                                                                                                                                                              cctctGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTAAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAg                                       <  1:656434/250‑1 (MQ=255)
                                                                                                                                                                                                                 ctGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTa                                     >  1:262463/1‑250 (MQ=255)
                                                                                                                                                                                                                   gATAACGCCTGAAATGAAGAAGTATCACCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATACCCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAAt                                   <  2:385607/250‑1 (MQ=255)
                                                                                                                                                                                                                             gaaaTGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTan                          <  2:327647/247‑2 (MQ=255)
                                                                                                                                                                                                                               aaTGAAGAAGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTTGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCa                      <  2:234061/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gaagaaGTATCATCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGAt                   <  2:262463/251‑1 (MQ=255)
                                                                                                                                                                                                                                     gaagTATCACCACGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAg                <  1:113697/251‑1 (MQ=255)
                                                                                                                                                                                                                                            caccaCGTTGGTCCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAAt          <  2:149719/247‑1 (MQ=255)
                                                                                                                                                                                                                                             atcaCGTTGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTg        >  2:664770/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   ttGGGTCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGCCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCgcagc  >  1:81902/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
CATCGTTGCATGGGGTGTTAACCAGAGAGCGACGGTAAAGGAGGCATCACATAAATATTCGCGCCATGTTTGAATCGTTTCATGACTTTTCAAACCACGGTAATCTTCTCGCGCTCCAAGGACTGTATGGTAATCACGCAACCCCGTTACAGACACACGACGATCGTCAAGAATGAGTTCATCGCAGCGCACTGCAAATTGCACACTCTCTGATAACGCCTGTAATGAAGAAGTATCATCACGTTGGATCCCAAGACAAGCCCCGAGTAGCCCTAATAACCCGCTTCGGGTCGGAAATCTTCCGGTAGGTCGCGTTCCTTCAAAGGTCGGCTGCCCCCAGGCTTGCATTGGCCCAGCAAGCCGCAAGATCAAATAAGATCTCATGTTCACGCCTCGCCATTATTACGAACCCAGGATTTTAACTGTTCTAAAGTAGGCATTTGTTTAACTTGAGCAGTAATTGGGTCTACATCAGATAAGCTGAATTGCGCAGC  >  NC_000913/2880665‑2881158

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