Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,734,830 C→T Q371* (CAG→TAG)  xylH → xylose ABC transporter membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,734,8300CT80.6% 128.9 / 27.1 73Q371* (CAG→TAG) xylHxylose ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (6/8);  new base T (25/33);  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 = 8.54e-01

CACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGCAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTCACCGCATCACATTACTGTTCAATGCCAATAAAGCCTATGACCGGCAGGTAGTAGAAGGCGTAGGGGAATATTTACAGGC  >  NC_000913/3734580‑3735073
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   
cACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGc                                                                                                                                                                                                                                                     <  1:35039/251‑1 (MQ=255)
   cAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGt                                                                                                                                                                                                                                                  <  1:525655/251‑1 (MQ=255)
   cAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCAGGGCTTCACTGGATAACGTCATGAGTATGATGTATGTACCGACCTTCTGGTAGt                                                                                                                                                                                                                                                  <  2:528768/251‑1 (MQ=255)
               gtTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCAGGGCTTCACTGGATAACGGCATGAGTTTGATGGATGTACCGACCTTCTGGTAGTATATCGTTaa                                                                                                                                                                                                                                        >  1:127987/1‑249 (MQ=255)
               gtTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCGCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGAGTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTAGCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTaaa                                                                                                                                                                                                                                       <  2:259168/250‑1 (MQ=255)
                  cgcgCTTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTTGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGGTAACGGCATGAGTATGATGGATGTACCGGCCTTCTGGTAGTCTATCGTTAAAgggg                                                                                                                                                                                                                                   >  2:125810/1‑249 (MQ=255)
                   gcgATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGc                                                                                                                                                                                                                                  >  1:310245/1‑251 (MQ=255)
                     gATTAACGGATTAATGGTAGCCATCTCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCTTTAAAGGTGCGa                                                                                                                                                                                                                                <  1:526663/251‑1 (MQ=255)
                         aaCGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATAGCGTTAAAGGTGCGATTCt                                                                                                                                                                                                                            <  2:293107/251‑1 (MQ=255)
                           cGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGAGTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCGTCTGGCAGTATATCGTTAAAGGTGCGATTCTGt                                                                                                                                                                                                                          <  2:99275/251‑1 (MQ=255)
                            ggATTAAGGGTTGCCAGCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTTCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTGCTGGTATTATATCGTTAAAGGTGCGATTCTGtt                                                                                                                                                                                                                         <  2:538928/251‑1 (MQ=255)
                             gATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGAGTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTg                                                                                                                                                                                                                        <  2:526201/251‑1 (MQ=255)
                               ttAAAGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCGTATGCGGGCAATATCGCCGAACTCGACGAAATTGCACCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGGTGCCGGAGCAGAAATGGGGGCCTTTATCATGGCTTCACTGGCGAACGGCCTGAGTATGTAGGATGTACCTACCTTCTGGTAGTATCTCGTTAAAGGGGCGATTATGTTGCt                                                                                                                                                                                                                      >  2:490077/1‑251 (MQ=255)
                                      tAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGGTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGGGTATGATGGATGTACCGACCTTCTGGTAGGATGTCGTTAAAGGTGCGATTCTGTTGCTTGCAGtt                                                                                                                                                                                                               >  2:195051/1‑250 (MQ=255)
                                          cATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTatgga                                                                                                                                                                                                           <  2:401092/251‑1 (MQ=255)
                                                   aTTAATCCGTAGTGCGCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCg                                                                                                                                                                                                  <  2:20028/251‑1 (MQ=255)
                                                         ccTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGAGTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCa                                                                                                                                                                                            <  2:481199/251‑1 (MQ=255)
                                                             aGTTCTCGACTTGGCGCTGGTTCCCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGCAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAAc                                                                                                                                                                                         <  1:118251/250‑1 (MQ=255)
                                                             aGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGTTGGGCTCCGCAACCAAACg                                                                                                                                                                                        >  2:338822/1‑251 (MQ=255)
                                                              gTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGTCGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGGAATGGGGGCATTTATCATTGCTTCACGGGATAACGTCATGAGTATGATGGATGTTCCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGc                                                                                                                                                                                       <  1:578365/251‑1 (MQ=255)
                                                                tctcGACTTGGCTCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGTCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGATAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGcc                                                                                                                                                                                      <  1:581645/250‑1 (MQ=255)
                                                                tctcGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGcc                                                                                                                                                                                      <  2:480219/250‑1 (MQ=255)
                                                                gctcGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGcc                                                                                                                                                                                      <  2:253459/249‑1 (MQ=255)
                                                                        tGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGa                                                                                                                                                                              <  1:364577/250‑1 (MQ=255)
                                                                         ggCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGCATTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGAtt                                                                                                                                                                            >  1:348167/1‑251 (MQ=255)
                                                                           cgcTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCATCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGAtttt                                                                                                                                                                          <  1:297663/251‑1 (MQ=255)
                                                                            tcTGGTTCACCTTCTGCGGGAAATCTCGCCGAACTGGACGCGATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGCGGGAAGCGTTGCCGGAGCAGTAAAGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGAAAGAAGGATGTACCGGCCTTCTGGTAGCATATCGTTAAAAGTGCGAGTCTGCTGCAGGCAATATGGTTGGACTCCGCAACAAAAAGCCGGTCTTGGTTTTg                                                                                                                                                                         >  2:436647/2‑251 (MQ=255)
                                                                               ggTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGa                                                                                                                                                                        <  1:360499/249‑1 (MQ=255)
                                                                                  tCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTACGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATaaa                                                                                                                                                                    <  1:5284/250‑1 (MQ=255)
                                                                                  tCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGACTGGCGGCACCAGCCTGCCTGGCGGTGTGGGAAGCGTTGCCGGCGCAGTAATGGGGGAATTTATAATGGCTTCACAGGATAACGGCATGAGCATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGCTTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATaac                                                                                                                                                                    >  1:324077/1‑249 (MQ=255)
                                                                                   cACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATaaaa                                                                                                                                                                   >  1:129592/1‑250 (MQ=255)
                                                                                   cACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATaaaa                                                                                                                                                                   <  1:230745/250‑1 (MQ=255)
                                                                                   cACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGCAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATaaaa                                                                                                                                                                   >  1:200093/1‑250 (MQ=255)
                                                                                    aCCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAAt                                                                                                                                                                 >  1:284751/1‑251 (MQ=255)
                                                                                        tcggCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGCAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAAtttt                                                                                                                                                              <  2:200093/247‑1 (MQ=255)
                                                                                        tCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAAtttt                                                                                                                                                              <  2:51254/250‑1 (MQ=255)
                                                                                             gggAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGCAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCaaa                                                                                                                                                        <  2:532525/251‑1 (MQ=255)
                                                                                                  atatCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCgg                                                                                                                                                   >  1:155612/1‑251 (MQ=255)
                                                                                                  atatCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTAGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGCAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCgg                                                                                                                                                   >  1:129869/1‑251 (MQ=255)
                                                                                                    atCGCCGAACTGGACGCAATTGCAGCATTCGTGATTGGCGGCTCCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGtt                                                                                                                                                 <  1:165418/251‑1 (MQ=255)
                                                                                                      cGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTa                                                                                                                                                <  1:125810/250‑1 (MQ=255)
                                                                                                         cGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGGGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGGCCTTCTGGTAGTATATCGTTAAAGGTGCGGTTCTGTTGCTGGCAGTATGGATGGACTCCGCAAACAAACACCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTTCGt                                                                                                                                             >  1:381330/1‑250 (MQ=255)
                                                                                                                gACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGCAGTATATCGTTAAAGGTGCGAGTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGGATTACCgg                                                                                                                                     <  2:402573/251‑1 (MQ=255)
                                                                                                                         gcagcaTGCGTGATTGGCGGCACCAGCCTGGCTGGCGGAGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGGGGATGATGGATGTACCGACCTTCTGGTAGGATATCGTTAAAGGTGCGGTTCTGTTGCTGGCAGGATGGATGGACTCCGCAACCAAAAGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGGTAAGTATTACCGGTATTGTGtt                                                                                                                            >  2:294102/1‑251 (MQ=255)
                                                                                                                                  gTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGCAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATg                                                                                                                    >  1:412956/1‑250 (MQ=255)
                                                                                                                                    gATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGACGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGAtt                                                                                                                 <  1:224098/251‑1 (MQ=255)
                                                                                                                                                        ctggcGGGGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGt                                                                                             >  1:89605/1‑251 (MQ=255)
                                                                                                                                                             agtgtGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGCTAACGGCATGAGTATGATGGATGTACCGAACTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACt                                                                                        >  2:166321/2‑251 (MQ=255)
                                                                                                                                                                   ggAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGt                                                                                  <  1:338822/251‑1 (MQ=255)
                                                                                                                                                                    gAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTc                                                                                 >  1:8270/1‑251 (MQ=255)
                                                                                                                                                                        ggTTGCCGGAGCAGTAATGTGGGCATGTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTTCCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTGGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTCTACGAAACGTCACCg                                                                             <  2:56579/250‑1 (MQ=255)
                                                                                                                                                                               ggAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTCACCGCATcac                                                                       >  1:612721/1‑250 (MQ=255)
                                                                                                                                                                                gAGCAGTAATGGGGGCAGTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGCAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTCACCGCATCACAt                                                                     <  2:139071/251‑1 (MQ=255)
                                                                                                                                                                                         tGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGACAAGAACCATGTTTACTAAACGTCACCGCATCACATTACTGTTc                                                             >  1:456403/1‑250 (MQ=255)
                                                                                                                                                                                         tGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTaaa                                                                                     >  2:454708/1‑226 (MQ=255)
                                                                                                                                                                                         tGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTaaa                                                                                     <  1:454708/226‑1 (MQ=255)
                                                                                                                                                                                         tGGGGGAATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTACAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTCACCGCATCACATTACTGTTc                                                             >  2:294065/1‑250 (MQ=255)
                                                                                                                                                                                              gCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAGAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTCACCGCATCACATTACTGTTCAATGcc                                                       <  1:21425/251‑1 (MQ=255)
                                                                                                                                                                                                   taacATGGCTTAACTTGATAACGGCATGAGTATGATGGATGTACCGACCTTCGGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAGCCAAACGCCGTTCTTGATTTTGATAAAAAGTTTCTCAAAGCCTGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATATTTACTAAACGTCACCGCATCACATTACTGTTCAATGCCAATa                                                   <  2:549312/247‑1 (MQ=255)
                                                                                                                                                                                                   tATCATGGCTTCACTGGATAACGGCATGCGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCCAACAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGAATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTCACCGAATCAAATTACTGTTCAATGCCAATa                                                   >  2:442364/1‑250 (MQ=255)
                                                                                                                                                                                                        tGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGCTGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTCACCGCATCACATTACTGTTCAATGCCAATAAAGCCt                                             >  1:386734/1‑251 (MQ=255)
                                                                                                                                                                                                           cTTCACTGGATAACGGCATGAGTATGATGGCTGTACCGACCTTCTGGCAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTCACCGCATCACATTACTGTTCAATGCCAATAAAGCCTATg                                          >  2:400045/1‑251 (MQ=255)
                                                                                                                                                                                                               aCTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGCAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTCACCGCATCACATTACTGTTCAATGCCAATAAAGCCTATGACCg                                      <  2:115638/251‑1 (MQ=255)
                                                                                                                                                                                                                          gCATGAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTCACCGCATCACATTACTGTTCAATGCCAATAAAGCCTATGACCGGCAGGtagtag                           <  2:155612/251‑1 (MQ=255)
                                                                                                                                                                                                                              gAGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTCACCGCATCACATTACTGTTCAATGCCAATAAAGCCTATGACCGGCAGGTAGTAGAAg                        <  1:195051/250‑1 (MQ=255)
                                                                                                                                                                                                                               aGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTCACCGCATCACATTACTGTTCAATGCCAATAAAGCCTATGACCGGCAGGTAGTAGAAGGc                      >  1:164379/1‑251 (MQ=255)
                                                                                                                                                                                                                               aGTATGATGGATGTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTCACCGCATCACATTACTGTTCAATGCCAATAAAGCCTATGACCGGCAGGTAGTAGAAGGc                      >  1:528226/1‑251 (MQ=255)
                                                                                                                                                                                                                                     atggatgTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTCACCGCATCACATTACTGTTCAATGCCAATAAAGCCTATGACCGGCAGGTAGTAGAAGGCGTAggg                <  1:463147/251‑1 (MQ=255)
                                                                                                                                                                                                                                           gTACCGACCTTCTTGCAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGAACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTATCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTCTCCGCATCACATTACTGTGCAATGCCAATAAAGCCTATGACCGGCAGGTAGTAGAATGTGTAGGGGAatat          <  1:78579/251‑1 (MQ=255)
                                                                                                                                                                                                                                           gTACCGACCTTCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCACACGACGTTCTTGATTTTGATAAAAATTTTCTCACAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTCACCGCATcaca                                                                      >  1:569591/1‑191 (MQ=255)
                                                                                                                                                                                                                                           gTACCGACCATCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATTGATGGACGCCGCAACCAAACGCCGTGCGGGAGTTTGATAAAAATGTTCTCAAAGCCGGTGACGTTTTACCGGTTTGGAGTTTTTGCATGATGCAGCAGGAAAAGAACCATGGTTACTAAACTTCACCGCATcaca                                                                      <  2:569591/191‑1 (MQ=255)
                                                                                                                                                                                                                                                gACCTTCTGGAAGAATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTAACCGCATCACATTACTGTTCAATGCCAATAATGCCTATGACCGGCAGGTAGTAGGAGGCGTAGGGGAATATTTAc      >  2:592042/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 aCCTTCTGGCAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTCCTAAACGTCACCGCATCACATTACTGTTCAATGCCAATAAAGCCTATGACCGGCAGGTAGTAGAAGGAGTAGGGGAATCTTTAc      >  2:233331/1‑249 (MQ=255)
                                                                                                                                                                                                                                                     tCTGGTAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTCACCGCATCACATTACTGTTCAATGCCAATAAAGCCTATGACCGGCAGGTAGTAGAAGGCGTAGGGGAATATTTACAGGc  <  1:55684/249‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   
CACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGCGTGATTGGCGGCACCAGCCTGGCTGGCGGTGTGGGAAGCGTTGCCGGAGCAGTAATGGGGGCATTTATCATGGCTTCACTGGATAACGGCATGAGTATGATGGATGTACCGACCTTCTGGCAGTATATCGTTAAAGGTGCGATTCTGTTGCTGGCAGTATGGATGGACTCCGCAACCAAACGCCGTTCTTGATTTTGATAAAAATTTTCTCAAAGCCGGTTACGTATTACCGGTTTTGAGTTTTTGCATGATTCAGCAGGAAAAGAACCATGTTTACTAAACGTCACCGCATCACATTACTGTTCAATGCCAATAAAGCCTATGACCGGCAGGTAGTAGAAGGCGTAGGGGAATATTTACAGGC  >  NC_000913/3734580‑3735073

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