Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,889,365 C→T K245K (AAG→AAA cysI ← sulfite reductase, hemoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,889,3650CT97.3% 204.6 / ‑1.6 73K245K (AAG→AAAcysIsulfite reductase, hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (31/40);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

GACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGAC  >  NC_000913/2889124‑2889611
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            
gACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGttt                                                                                                                                                                                                                                                     >  1:100880/1‑251 (MQ=255)
 aCGCTCCTCTTCCGCTTTCAACGTCTATACCCCCTCGCTCTCCAGCTGGTTTTTGGTTGTGGCATTTTTACGATCGGTTCTGTTACCCCAGTCTCGCTGAGTTGTCCCGACGTCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCTCCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGtttt                                                                                                                                                                                                                                                    <  2:642335/251‑1 (MQ=255)
 aCGCTCCACTTCCCCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTCCGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGtttt                                                                                                                                                                                                                                                    <  1:867323/251‑1 (MQ=255)
  cgtcccACTGCCGCTTTAAACGTCCCAACCCCCACGCGTTCCAGCGTGTATTGGGGTTGGGCATTTTTACGATCGGTATGGTTCCCCCAGTTACGCTGATTTGTCGCGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTc                                                                                                                                                                                                                                                   <  1:857750/247‑1 (MQ=255)
   gcgccACTTCCGCTTTAAACGTCACACCCCCCACTCGCTCAAGCGTGTATTTGGTTTGGGCAGTTGTACGTTCGGGTCGGTTACCCCAGTCCCGTGGAGTTGTCCGAACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCGCCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCg                                                                                                                                                                                                                                                  <  1:283167/248‑1 (MQ=255)
      ccACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCg                                                                                                                                                                                                                                               >  1:165527/1‑251 (MQ=255)
           tCCGCTTTAAACGTCTCAACCCCCCCGCGCTCCGGCGTGTATTTGGTTGTGGCATTTTTACGGTCGGTTCGGTTCCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCTTCAGATAGCCAAACTCGCTCGCGGGGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCg                                                                                                                                                                                                                                           <  1:858910/250‑1 (MQ=255)
           tCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAGGCCCCCCAGTTTGCCGTTTTCGGCGATCGc                                                                                                                                                                                                                                          >  1:724451/1‑251 (MQ=255)
               cTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCACCAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGGTAGCCACg                                                                                                                                                                                                                                      >  2:268677/1‑251 (MQ=255)
                    aaCGTCTAAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGGTGGTTTTCTTGTTGACGTGTTCGATGGAAAGCCCACGGCCCCCCAACAGGTTTAAGACCACAAGTTTGCCGTTTTCGGCGATCGCCACGa                                                                                                                                                                                                                                     >  2:180372/1‑247 (MQ=255)
                       gTCTCCACCCCCAGTCGCTCCAGCGTGTATTGGGTTTGGGTATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCTCCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCAt                                                                                                                                                                                                                             <  1:569637/251‑1 (MQ=255)
                           cAACCCCCACGCGCTCCAGGGTGTATTTGGTTTTGCCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGACACGAAGTTCATGTCg                                                                                                                                                                                                                        <  1:217643/251‑1 (MQ=255)
                                       gcTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAgg                                                                                                                                                                                                           >  1:723549/1‑250 (MQ=255)
                                          cccgCGTGTTTTTGTTTTTGGCATCTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGTCGTGTTCGATGGAAAGCCCTCCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCg                                                                                                                                                                                                        <  1:449517/248‑1 (MQ=255)
                                                   aTTTGGTCTTGGCTTTTTTACGATCGGTTCGGGTAACCCAGTCACGCTGAGTTGTTACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGACCAACTCACTCGCGGTGCGGGCGGAGGTTTTCTTGTTGCCGAGTTCGATTGAAAGCCCACCGCCCACCACCAGGTTAAAGCCCAACAGTTTGCCGTTTTCGGCGATCGCCACGGAGTTAATGTCGTTGGCGTGcgtatagcaatcgttt                                                                                                                                                                                               >  2:368291/1‑236 (MQ=255)
                                                      tggtttggtctttttttCGATGGGTTCGGTTCCCCCATGCCCGCTGAGTTGTCCCGACGGTGTCGGCCGCCGCCAGCGTAGGCTCCTGCGGCATATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCTCCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCtg                                                                                                                                                                                             <  2:734702/239‑1 (MQ=255)
                                                       ggTTTTGTCTTTTTTTCGTTGGGTGCGGTTCCCCCAGTCCCGCTGAGTTGGCACGCGGCCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCtgtg                                                                                                                                                                                           <  2:81670/251‑1 (MQ=255)
                                                       ggTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGATTCGGCCACCGCCAGCGAATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACAAACAGGTTAAAGACCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGATCAGTGCGTTGGCGTGccaataagtatcgttctt                                                                                                                                                                                             >  2:181066/1‑232 (MQ=255)
                                                        gTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCTTTCTGTgg                                                                                                                                                                                          >  2:35057/1‑251 (MQ=255)
                                                          tttGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCTCGACGGCTTCGGACACCGCCAGCGTCTGCTCAAGCGGCAGTTAGCCAAACTACATCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGACCACCGCCCACCAACAGGCTAAAGCCCACCAGTTTGCCGTTTACGGCGCTCGCCACGACAGTCATGTAGTAGGCGCTGACATCGATATAGTTCTGTGGCg                                                                                                                                                                                        >  2:363806/1‑251 (MQ=255)
                                                          tttGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCAGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGCAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTgcggcg                                                                                                                                                                                        >  1:521246/1‑251 (MQ=255)
                                                             ggCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCgggg                                                                                                                                                                                     >  1:309822/1‑250 (MQ=255)
                                                                      aCGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGACGTTCATGTCGTTGGCGTGCAGATCGATATCTTTTTGTGGCGGGagttcacccc                                                                                                                                                                           >  2:437248/1‑242 (MQ=255)
                                                                        gATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGCCGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCCAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCCACAGGTTACAGCCCACCAGTTTGCCGTTTTCGGCGGTCGCCCCGACGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACCACCGTg                                                                                                                                                                        >  2:338220/1‑251 (MQ=255)
                                                                         aTCGGTTCGGTTACCCCAGCCCCGCTGAGTTGTCACGACGGCTTCGGCCTCCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTgg                                                                                                                                                                       <  1:810563/251‑1 (MQ=255)
                                                                          tcggttcggtTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCgggggtcatcccctggg                                                                                                                                                                        >  2:228784/1‑237 (MQ=255)
                                                                            gggtcggttACCCCAGTCACGTGGATTTGCCATGAGGGCTTCGGCACCCGCCAGCGGATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTGCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGttt                                                                                                                                                                   <  2:612494/248‑1 (MQ=255)
                                                                              ttcggttACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTCGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACCGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGCAGGTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCggggagcaaaaccggagtttg                                                                                                                                                                    >  2:557371/1‑233 (MQ=255)
                                                                               tcggttACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGAACACttccctgggtttt                                                                                                                                                                   >  1:60563/1‑238 (MQ=255)
                                                                               ccggttACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCCGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCACACCAACCGGTTCAAGCCCACCCGTTTGCCGATTTCGGCGATCGCCACGCAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTTTGGCggaggcggaacctgggggtttg                                                                                                                                                                   >  2:368773/2‑231 (MQ=255)
                                                                                         cAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACAAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGCGGTTTTGTATTTATGCgg                                                                                                                                                      >  1:270900/1‑251 (MQ=255)
                                                                                              aCGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGAAGCTAGCCAAACTCACTCACGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGACAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCAATCACCGGGGTTTGAATTTAcgccgaaggg                                                                                                                                                 >  1:780192/1‑244 (MQ=255)
                                                                                                    aGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCAATACCCTGGGTTTTGATTTTCGGGGCaggtggggct                                                                                                                                           >  1:468788/1‑245 (MQ=255)
                                                                                                       gtggtcaggccgtctccGGCCACCGCCAGCGTATGCGCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTg                                                                                                                                          <  1:673691/235‑1 (MQ=255)
                                                                                                           tccctccggtttcggcacccGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACGCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGcc                                                                                                                                       <  2:297305/232‑1 (MQ=255)
                                                                                                              cgacggttccgcccACCGCCAGCGTATGTTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTGCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCATCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGAGCACTACCGTGGGGTTGAATTTACGCGGCAGGTAGGTCTGTCCGAgg                                                                                                                                   <  2:549445/239‑1 (MQ=255)
                                                                                                                cgGCTTCGGCCCCCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACGCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATc                                                                                                                                <  2:165527/251‑1 (MQ=255)
                                                                                                                  gCTTCGGCAGCCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTCCCGTGGTTTTGAATTTCCGCGGCAGGTAGGTCTGGCCGAGGATCgg                                                                                                                              <  1:41356/251‑1 (MQ=255)
                                                                                                                   ggtgcGGCCACCGCCAGCGTATGCTCAAGCGGCAGATAGCCAACCTCACGCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGTGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCgg                                                                                                                              <  2:276765/247‑1 (MQ=255)
                                                                                                                     ttCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCAATACCGTGGTTTTGAATTTAAGCGGGCGGTTGGGCTGGCCGGGGAtggggg                                                                                                                            >  1:577264/1‑246 (MQ=255)
                                                                                                                       cGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTAACGGGGGTTTGGATTTACGCGGCAGGGAGGGCTGGGCGAGGGTCgggttt                                                                                                                          >  1:870869/1‑247 (MQ=255)
                                                                                                                         gacgcgcgcggctttgtcTCCAACGGCTGATAGCCCAACGCACTTGCGGTGCTGGCGTAGGGTTTCTTGTTGCCGTGCGCGATGGAGAGCCCACCGCCCGCAAACAGGTTAAAGCCCACCCGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCATGTAGGTCTGGCCGAGGATCGGttcttc                                                                                                                        <  2:595214/234‑1 (MQ=255)
                                                                                                                          cccccGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTtca                                                                                                                       <  2:734652/248‑1 (MQ=255)
                                                                                                                           ccccgacAGCGTATGCTCCAGCGGCCGATAGCCAAACTCCCTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCCCCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTtcat                                                                                                                      <  1:446851/245‑1 (MQ=255)
                                                                                                                              cGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGGGGTTTTGAATTTACGCGGGCGGGAGGGATGGCCGGGGAGCGGTTCTtcttca                                                                                                                    >  1:591177/1‑250 (MQ=255)
                                                                                                                                 cAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTTACTACCGTGGGTTTTAATTTTTGCGGCAGGGAGGGCTGGCCgggggtagggtcttctttcgt                                                                                                                  >  1:443166/1‑229 (MQ=255)
                                                                                                                                         gCTCCAGCGGCAGATATCCAAACGCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCag                                                                                                                   <  2:599951/240‑1 (MQ=255)
                                                                                                                                         gCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGGTTTTAATTTACGCGGCAGGGAGGGCTGGCCGCGGATCGGTTCTTTATCag                                                                                                                   >  1:599951/1‑240 (MQ=255)
                                                                                                                                            ccAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGGGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACttt                                                                                                      <  2:121693/250‑1 (MQ=255)
                                                                                                                                               gCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTcc                                                                                                  <  1:228784/251‑1 (MQ=255)
                                                                                                                                               gCGGCAGAGAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTcc                                                                                                  <  1:353713/251‑1 (MQ=255)
                                                                                                                                                       aaGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAg                                                                                          <  1:724296/250‑1 (MQ=255)
                                                                                                                                                         ggccaaacccaaccgcgGTGCGGGCGGTGGTTTGCTTGTTGCCGTGTTCGATGGAAAGGCCTCCGCCCACCATCAGGTTAAAGCCCACCAGTGTGCCGTTTTCGTCGATCGCCACTAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGc                                                                                         <  2:267499/238‑1 (MQ=255)
                                                                                                                                                         gCCAAACTCACGCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGcc                                                                                        <  1:35057/251‑1 (MQ=255)
                                                                                                                                                           caacctcactcGCGGTGGGGGCGTAGGTTTTATTGTTGACGTGTTCGATGGAAAGCACGCCGACCGCCAGCAGGTTGAAGCTCACCAGTTTGCCGTTATCGGCGAGCGCCGCGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTCCCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTCCCTGGTCGAGCCAg                                                                                      <  1:145042/247‑1 (MQ=255)
                                                                                                                                                              actcactcGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGAtc                                                                                   <  1:593214/251‑1 (MQ=255)
                                                                                                                                                                 cactcGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTTAATTTTCGCGGCAGGTAGGGCTGGCCGAGGGTCGGTTCTTCATCTGTAGTGGCGACTTTTTTCTTGGCGAGCCAGATCTcc                                                                                >  1:433492/1‑251 (MQ=255)
                                                                                                                                                                  cctcGCGGTGCGGGCGTAGGTTTGCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGTCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCg                                                                               <  1:557371/250‑1 (MQ=255)
                                                                                                                                                                      cgcgGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCCATACCGTGGGTTTGAATTTATGCGGCCGGGAGGGCTGGGCGAGGGTAGGGTTTTCCACCGAAGTGGCGGCTTTTTTCTGGGCGGGCCACATCTCCGCac                                                                            >  1:255625/1‑250 (MQ=255)
                                                                                                                                                                      cgcgGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGCTCACTACCGTGGTTTTTTATTTACGCGGCAGGGAGGGCTGGCCGCGGGTCGGGTCTTCATTAGTAGGGGGGGCATTTTCCTTGTTGGGCCAGATCTCCGCAt                                                                            >  1:72015/1‑251 (MQ=255)
                                                                                                                                                                                gCGTAGGTTTTCTTGTTTCCGTGTTCGATGCAAACCCCACCGCCCACCAACAGGTTACAGCCCACCAGTTTTCCGTTTCCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCTCATACGCGCGGGt                                                                  <  2:657461/251‑1 (MQ=255)
                                                                                                                                                                                   tatgTTTTCTTGTTGCCGTGTTCGATGGACAGCCCACCGCCCACCAACAGGTTAAAGCCCTCCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGGTCTGTGGCGGGAGCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGGTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTAc                                                                <  2:784314/247‑1 (MQ=255)
                                                                                                                                                                                                  ccGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGATACACCCGGGGTTTTGAATTTTCGCGGGAGGGAGGGCTGGGCGGGGATCGGGTTTTCATCTGAAGGGGGGGCTTTTTTCTGGTCGAGCCCGATATCCGCCTACGCCCGGGGAGGAGGCcacaaaagcgc                                                >  2:265041/1‑244 (MQ=255)
                                                                                                                                                                                                      gtTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCCGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTAACGTGGGTTTGGATTTACTCGGGAGGGAGGGCTGGCCGGGGGGCGGGTATTTATTAGGAGAGGGGGCTTTTTCCTGGTCGGGCCAGATCTCCGCATAAGGGCGGGGACGGGGACACAGATggccaga                                             >  1:718718/1‑244 (MQ=255)
                                                                                                                                                                                                         cGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTGTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGTCAGGTAGGTCTGGCCGAGTATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAAt                                          <  1:406523/250‑1 (MQ=255)
                                                                                                                                                                                                             ggAAAGCCCCCCGCACACCTACAGGTTAAAGCCCACCAGTTTGCCGTGTGCGGCGATCGACACGAAGTTCTGGTCGTTGGCGTGCAGATCGATATCGGTCTGGGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGTCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCATCAGATGCTCAGAAAtcttc                                      <  2:100880/250‑1 (MQ=255)
                                                                                                                                                                                                             ggAAAGCCCACCGCCCACCAACAGGTTAAAGCCCCCCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCTCTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTTGTGGCGACTTTTTCCTGTTCGAGCCAGATTTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAAtcttc                                      <  2:309822/250‑1 (MQ=255)
                                                                                                                                                                                                                    ccaccgcccaccAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTTAATTTTCGCGGCAGGTAGGTCTGGCCGGGGATCGGGTCTTTATCAGGAGGGGCGACTTTTTTCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCCGAAATCTTCTTCGcccc                              >  1:348374/1‑250 (MQ=255)
                                                                                                                                                                                                                         gcccaccAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGt                         <  1:237734/251‑1 (MQ=255)
                                                                                                                                                                                                                         gcccaccAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGt                         <  1:217702/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTg            <  1:814668/250‑1 (MQ=255)
                                                                                                                                                                                                                                           gCCCACCAGTTTGCCTGTTGCGGCGATCCCCACGATGTTCTTGTCGATGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAAGTTTCGCGGCAGGTAGTTCGGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCATCAGATGCGCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCt       <  2:333631/251‑1 (MQ=255)
                                                                                                                                                                                                                                                ccaGTTTGCCGTTGTCTGCGATCTCCACTATGTTCATGTCGTTGGCGTGCCGCCCGGTATCGCTCTGTGTCGGGATCACAACCGTGGTTTTGAATTTGCGCGGCAGGTAGGTCTGGCCGAGGATCGTTTCTTCGTCAGTAGAGGCGACTTTTTCCTGGGCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCGCAGAAATCTTCCTCGCCCACTCGTACGCTTCCGCGTGCAgttgtgac  <  2:206561/251‑8 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            
GACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGAC  >  NC_000913/2889124‑2889611

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