Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,164,123 C→G A162G (GCA→GGA)  btuB → cobalamin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,164,1230CG97.7% 282.3 / ‑3.1 87A162G (GCA→GGA) btuBcobalamin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base G (38/47);  total (40/47)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTGCGCGGCTATGGCTATG  >  NC_000913/4163879‑4164368
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
ggTCAACTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGATCCCGGAACGGAAATTTCAGGAGGGTg                                                                                                                                                                                                                                                  <  1:362407/251‑1 (MQ=255)
  tCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTggg                                                                                                                                                                                                                                                >  1:228126/1‑251 (MQ=255)
  gcAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTggg                                                                                                                                                                                                                                                <  1:911961/250‑1 (MQ=255)
    aGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAAACCGGAACGGAAATTTCAGGAgggggggg                                                                                                                                                                                                                                               >  1:773431/1‑250 (MQ=255)
        ctcATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGTCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGGGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAg                                                                                                                                                                                                                                            <  2:267093/249‑1 (MQ=255)
           aTCTATTTTTATTCGCGGTACAAAAGCCAGTCATGTGATGGTGTTAATTGATGGAGTACGCCTGAATCTGGCGGGGGTGAGTTGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCCGCGTGTTGAATATATCCGTGGGCCGCGCTCAGCTGTTTATGGGTCCGAGGCAATAGGCGGGGAGGTGAGTATCATCACAAAGCGCGCTGAACTCGGACCAGAAATTTAAGGAGGGTGGGGaaacact                                                                                                                                                                                                                                        >  1:679635/1‑245 (MQ=255)
            tCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGGTCCGGTGCAATAGGCGGGGGTGTGAATATCAGTGCGAAGAGCGATGACAACGGAACGGAAATTTTAGGGAGGGTGGGGagcagt                                                                                                                                                                                                                                         >  2:534381/1‑244 (MQ=255)
             cTATTTTTATTCTCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGGGGGGAAGCAATAGt                                                                                                                                                                                                                                    >  1:169534/1‑251 (MQ=255)
               aTTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCTTGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTTAATATATCCGTGGGCCGCGCCCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATAGCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGtt                                                                                                                                                                                                                                   <  1:593675/250‑1 (MQ=255)
               aTTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAg                                                                                                                                                                                                                                     <  1:822565/248‑1 (MQ=255)
               aTTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGAGGGGAAGCAATAGttt                                                                                                                                                                                                                                  >  1:267490/1‑250 (MQ=255)
                tttttATTCGCGGTACAACTGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAAAATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCCGGAGGGTGGGGAAGCAATAGtttt                                                                                                                                                                                                                                 >  1:59985/1‑249 (MQ=255)
                 ttttATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATc                                                                                                                                                                                                                                >  1:468486/1‑251 (MQ=255)
                   ttATGCGCGGTACAAATGCCAGTTATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTGCCCGATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGTTCCGCTGTGTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGATGCAATAGTTATCa                                                                                                                                                                                                                               <  2:749685/250‑1 (MQ=255)
                      ttCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGGAATTTCAGGAGGGGGGGGAAGCAATAGTTATCAGAAc                                                                                                                                                                                                                           >  1:325827/1‑251 (MQ=255)
                        cgcgGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGGAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACt                                                                                                                                                                                                                          >  1:49908/1‑250 (MQ=255)
                             tACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGGTGAACCCGGAACGGGAATTTCAGGGGGGGGGGGAAGGAATAATTATCACAAATatgggt                                                                                                                                                                                                                    >  2:945693/1‑248 (MQ=255)
                                   tGCCAGACATGTGTTGGAGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCAAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTTTGGTTCCGATGCAATCGGAGGGGGGGGGGAACTTATCACTAGCGGCGATGAACCCGGAACGCGAATTTCAGGTGGGGGGGGAAGCAATCAGTTCTAAAACTATGATGTatctt                                                                                                                                                                                                              >  2:149009/1‑250 (MQ=255)
                                     ccAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACg                                                                                                                                                                                                            <  1:51784/251‑1 (MQ=255)
                                       aGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCCGAACTATGATGTCTCTACgc                                                                                                                                                                                                           >  1:633550/1‑250 (MQ=255)
                                          cATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAAACCGGAACGGAAATTTCAGGGGGGTGGGGAAGCAAAAGGTTTCAGAACTAAGATGtttttacacagc                                                                                                                                                                                                        >  2:165061/1‑239 (MQ=255)
                                            tgtgTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaac                                                                                                                                                                                                     >  1:229370/1‑251 (MQ=255)
                                             gtgtTGGTGGTAATTGATGGCGGCCGCCTGAATCTGGCGGGGGTGAGTGGTTCGGCCGACCTGAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCTATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGcaaca                                                                                                                                                                                                    <  2:768808/251‑1 (MQ=255)
                                                  ggTGTTAATTGATTGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTg                                                                                                                                                                                                <  2:733536/250‑1 (MQ=255)
                                                       tAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGCTGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCCGGAGGGTGGGGAAGCAATAGTTATCAGACCTATGATGTCTCTACGCAACAACAACTGGGGGa                                                                                                                                                                                           >  1:213856/1‑250 (MQ=255)
                                                        aaTTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGAt                                                                                                                                                                                          <  1:634150/250‑1 (MQ=255)
                                                            gATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGac                                                                                                                                                                                     >  1:623385/1‑251 (MQ=255)
                                                               ggCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACg                                                                                                                                                                                  <  2:169534/251‑1 (MQ=255)
                                                                      gCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGGATATAATCACGACGCGCGATGGACCCGGGACGGAAATTTCAGGAGGGTGGGGAAGCAAAAGGTATCAGAACTATGATGTCTCTACGCAGAAAAAACTGGGGGGAAAGACACGGGTAACg                                                                                                                                                                           >  2:349467/1‑251 (MQ=255)
                                                                          gAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCCCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGAGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCt                                                                                                                                                                         <  1:534381/249‑1 (MQ=255)
                                                                           aaTCTGGCTGGGGTGTGGGGTTCTGCCGCCCTTATCCTGCTCCCTTGTGCGCTTGTCCAGCGTGTTGAAGATATCCGGGGGCCGCTCTCAGCTGTTTATGGTGCCGATGCAATAGGCGGGGTGGTGAATATCATCACTACGTGCGATGAACTCGGAACGGAAATTGCAGGAGGGTGGGGAAGCAATAGTAATCATAACTATGATGTCTCTACGCAGCAACAACGTGGGGAGAAGACACGGGTAACGCTGtt                                                                                                                                                                      <  2:558778/251‑1 (MQ=255)
                                                                                ggCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGCTGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCCATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGGAAAGACACGGGTAACGCTGTTGGGCg                                                                                                                                                                 >  2:330393/1‑251 (MQ=255)
                                                                                ggCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTGGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGAAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTggg                                                                                                                                                                   <  1:72567/249‑1 (MQ=255)
                                                                                  ggggggTGAGTGGTTCTGCCGACCTGAGCCAGTTCCCTATTGCGCTTGTCCAGCTTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGAt                                                                                                                                                               <  2:315918/250‑1 (MQ=255)
                                                                                         gAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGcc                                                                                                                                                         <  1:259083/250‑1 (MQ=255)
                                                                                             ggttaggCCGACATGAGCCAGTTCGCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATAc                                                                                                                                                    <  2:151510/245‑1 (MQ=255)
                                                                                                   gCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATgg                                                                                                                                              <  2:371722/251‑1 (MQ=255)
                                                                                                    ccGACCTTAGGCAGTTGCCTATTGGCCTTTTCCAGGGTGTTGAAAATAATCGTGGGACGCGGCACGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGGTGATATCATCACGACGCGGGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGGCATAGTTAACAGAACTATTGTGTCTCTACGCAGCAACCACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGt                                                                                                                                             <  2:193379/251‑1 (MQ=255)
                                                                                                      gACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGtt                                                                                                                                            >  1:409505/1‑250 (MQ=255)
                                                                                                           tAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTatgatg                                                                                                                                      >  1:403461/1‑251 (MQ=255)
                                                                                                           tAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTatgat                                                                                                                                       <  1:346373/250‑1 (MQ=255)
                                                                                                            aGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTTAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGAtgt                                                                                                                                     <  1:733173/251‑1 (MQ=255)
                                                                                                              ccAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGTGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGAtgttt                                                                                                                                   >  1:247312/1‑250 (MQ=255)
                                                                                                               cAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGAtgttgt                                                                                                                                  <  1:47187/251‑1 (MQ=255)
                                                                                                                    cccTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGcc                                                                                                                              <  1:986381/250‑1 (MQ=255)
                                                                                                                    cccTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCt                                                                                                                             <  1:952834/251‑1 (MQ=255)
                                                                                                                    ccatATTGCGCGTGTCCAGCGGGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCt                                                                                                                             <  1:18656/248‑1 (MQ=255)
                                                                                                                     ccTATTGCGCTTGTCCAGCGTGTTGCATATACCCGTGGGCCCCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGACTATCATCACCACGTGCCATGACACCGGAACGGCAATTTCAGGAGGGGGGGGCAGGCATCGTTATTCGAACTATGAACTCTCTACACCACAAACACTGGGGGCTAAGTCCCGGGCATCGCAGTTGGGCGATTATGCCCATACTCATGGTTATTAtgtttgttccc                                                                                                                             >  2:120594/1‑243 (MQ=255)
                                                                                                                     ccTATTGCGCTTGGCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATGGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTa                                                                                                                            <  1:782173/251‑1 (MQ=255)
                                                                                                                        aTTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGAAACGGGTAACGCTGTTGGGCGATTATGCCCATACTGATGGTTATGATGTTGTTGCCTATgg                                                                                                                         >  1:6477/1‑251 (MQ=255)
                                                                                                                            cgcTTGTCCAGCGTGTTGAATATATCCGTTGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGTGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAAt                                                                                                                     <  1:904516/251‑1 (MQ=255)
                                                                                                                                   ccAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGa                                                                                                               <  1:753255/250‑1 (MQ=255)
                                                                                                                                   ccAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGTGAATATCATCACGACGCGCGCTGAACCAGGAACGGAAATTTCAGGAGGGGGGGGAAGAAATAGTTATCAGAACTATGATGTCTCTACGGCGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGGTGCCTATGGTAATACCGGa                                                                                                               >  1:464519/1‑250 (MQ=255)
                                                                                                                                     aGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACg                                                                                                            >  1:399974/1‑251 (MQ=255)
                                                                                                                                     aGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGGGGGGAATATCATCACGACGCGCGCTGCACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACAAGGGTAACGCTGTTGGGCGATTATGCCCCTACTCCTTGTTATGATGTTGTTGCCTATGGTAATACCGGaaac                                                                                                            >  2:398263/1‑249 (MQ=255)
                                                                                                                                        gtgtTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAAACCGGAACGGAAATTTCAGGAGGGGGGGGAAGCAATAGTTATCAGAACTATGATGTCTTTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTTTGATGTTGTTGCCTATGGTAATACCGGAACGCa                                                                                                          >  1:509871/1‑250 (MQ=255)
                                                                                                                                             gAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTa                                                                                                                                                             <  2:51097/194‑1 (MQ=255)
                                                                                                                                             gAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTa                                                                                                                                                             >  1:51097/1‑194 (MQ=255)
                                                                                                                                               atatatCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGcaga                                                                                                  <  1:610619/251‑1 (MQ=255)
                                                                                                                                                     ccGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGAt                                                                                             <  1:667358/250‑1 (MQ=255)
                                                                                                                                                        tGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACg                                                                                         >  1:304415/1‑251 (MQ=255)
                                                                                                                                                         gggCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGa                                                                                        <  1:540839/251‑1 (MQ=255)
                                                                                                                                                           gtcgtgcTCCGCTATTTAGGGGTTCGATGCAATAGGAGGGGTGGTGAATATGATCACGACGCGCGATGAACCCGGAACGTAAATTTCAGGAGGGGGGGGAAGCTATAGTTATCAGAACTATGATGTCGCGACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGTTAATACCGGAACGCAAGCGCAGACAGATAACTATg                                                                                      <  2:696212/246‑1 (MQ=255)
                                                                                                                                                           gCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATg                                                                                      <  1:476030/251‑1 (MQ=255)
                                                                                                                                                               cgctgcgctgTTTATGGGTCCGATGCAAGAGGCGGGGTGGTGAATATCATCACGACGCGCGTTTAACCCGGAACGGAAATTTCAGGAGGGTGGTGAAGCAATAGTTATCAGAACTATGATTTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTAGGCCCATACTCATGGTTATGATGTTGTTGCCTATTGTAATACCGGAACACAAGCGCAGACAGATAACGATGtttt                                                                                  <  2:272447/251‑1 (MQ=255)
                                                                                                                                                               cgcTGCGCTTTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTGCAGGAGGGTGGGGAAGCAATAGTGATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGTGTAACGCTTTTGGGCGATTATGCCCATACTCATGGTTATTATTTTGTTGCCTATGGTAATCCCGGAACGCAAGCGCAGACAGATAACGATGGttt                                                                                  <  1:552582/251‑1 (MQ=255)
                                                                                                                                                                 cTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGttttt                                                                                >  1:316868/1‑251 (MQ=255)
                                                                                                                                                                    cGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTaa                                                                             <  1:370069/251‑1 (MQ=255)
                                                                                                                                                                        gTTTATTGTTCCGATGCAATAGGCGGGGCGGTGAATATCATCACTCCGCGCGATGAACCCGGAACGGATATTTCAGGAGGGTGTGGAAGCACTAGTTATCAGCACTATAGTGTCTTTACGCAGCAACAACTGCGGGATAAGACACGGGTAACGCTGATGGGCGATTATGCCCATAATCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAGCGCACGGGCAGTCAGATAACGATGGTTTTTTAAGTa                                                                          >  2:376330/1‑250 (MQ=255)
                                                                                                                                                                            aTGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAAc                                                                      <  1:699424/250‑1 (MQ=255)
                                                                                                                                                                            aTGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAAc                                                                      <  2:228126/250‑1 (MQ=255)
                                                                                                                                                                                  ccGATGCAATTGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCACATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTAt                                                               <  1:178386/251‑1 (MQ=255)
                                                                                                                                                                                   cGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATg                                                              >  1:980281/1‑251 (MQ=255)
                                                                                                                                                                                      tGCAATAGACGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGTACGGCAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCCACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGTTTATGCCCATACTCATGGGTATGATGTTGTTGCCTATGGGTATACCGGAACCCAAGCGCAGACCACTAACGATGGTTTTTTAAGTAAAAACCTTTATTgcg                                                           >  2:699973/1‑251 (MQ=255)
                                                                                                                                                                                          aTAGGCGGGGGGGTGAATATCATCTCGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCt                                                       <  2:805358/251‑1 (MQ=255)
                                                                                                                                                                                              gCGGTGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAg                                                   <  2:633550/251‑1 (MQ=255)
                                                                                                                                                                                               cGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGc                                                  >  1:467962/1‑251 (MQ=255)
                                                                                                                                                                                                        tgAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCCGCCGGGTGCGGCAGCAATAGTTATCAGAACTATGATGGCCCGACGGAGCAACCACTGGGGGATAAGACACGCGTAACGCTGTTGGGCGATTAGGCCCATCCTCACGGTTATGATGGTGTTGCCTATGGTAATACCGGTCAGCAAGCGCAGGCAGATATCGATGGTTTTTTAAGTAAAACGCTGTATGGCTCGCTGGAGCATAACttt                                          >  1:473233/1‑251 (MQ=255)
                                                                                                                                                                                                          aaTATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTAc                                        >  1:203657/1‑251 (MQ=255)
                                                                                                                                                                                                           atatCAGCACGACGCGCGATGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACt                                       <  1:859528/251‑1 (MQ=255)
                                                                                                                                                                                                                              aTGAACCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCttt                     <  1:828911/251‑1 (MQ=255)
                                                                                                                                                                                                                              aTGAAACCGGAACGGCAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCttt                     >  1:275215/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aaCCCGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTtgtg                  >  2:906006/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ccGGAACGGAAATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTgcgc               >  2:187821/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aacggaaATTTCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTgcgagtct           <  2:623385/251‑6 (MQ=255)
                                                                                                                                                                                                                                                 ttCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTGCGCggctatggcta  <  1:328782/250‑2 (MQ=255)
                                                                                                                                                                                                                                                  tCAGGAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAGGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATTCCGGAACTCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTGCGCGgctatggctatg  <  1:561686/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       
GGTCAGCTCTCATCTATTTTTATTCGCGGTACAAATGCCAGTCATGTGTTGGTGTTAATTGATGGCGTACGCCTGAATCTGGCGGGGGTGAGTGGTTCTGCCGACCTTAGCCAGTTCCCTATTGCGCTTGTCCAGCGTGTTGAATATATCCGTGGGCCGCGCTCCGCTGTTTATGGTTCCGATGCAATAGGCGGGGTGGTGAATATCATCACGACGCGCGATGAACCCGGAACGGAAATTTCAGCAGGGTGGGGAAGCAATAGTTATCAGAACTATGATGTCTCTACGCAGCAACAACTGGGGGATAAGACACGGGTAACGCTGTTGGGCGATTATGCCCATACTCATGGTTATGATGTTGTTGCCTATGGTAATACCGGAACGCAAGCGCAGACAGATAACGATGGTTTTTTAAGTAAAACGCTTTATGGCGCGCTGGAGCATAACTTTACTGATGCCTGGAGCGGCTTTGTGCGCGGCTATGGCTATG  >  NC_000913/4163879‑4164368

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