Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,2760TC23.5% 133.7 / 11.7 69pseudogene (885/2518 nt)ydbAputative outer membrane protein, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base T (33/19);  new base C (15/1);  total (49/20)
Fisher's exact test for biased strand distribution p-value = 2.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGG  >  NC_000913/1466028‑1466520
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
cacaACTGTGTGGTTGTTATTACGTGTGAGAATGCCACAGTAAAAAACAATGCCAAAACCACCGTTGACGGTAAATATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTTGATGTCAGCGGCGTCGTTCACGGTATTGATATCGCCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGGACGTTACCGATCCGGTGTCCATGGGTATCCAGATCGACGGTg                                                                                                                                                                                                                                                       <  2:229445/250‑1 (MQ=18)
 acaACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCCAACCCGCGTCAATGGGTATCCAGATTGACGGCGAc                                                                                                                                                                                                                                                     >  1:60637/1‑251 (MQ=18)
        tGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACACTGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGGCACGGTATTGATATCAACGGAGACAGCGCGACGGTGGATCACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGGATCCAGATCGACGGTGCCAAGGccc                                                                                                                                                                                                                                              >  1:290911/1‑250 (MQ=18)
         gAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCCGCGGGAACGGTGTTGAGAGCACGGGTGTCAGCCCCGGGGTGGATAACATGGTCACCAAGACCGTCACCGAACCGGAGTCCNGGGGTATCCAGATCGACGGAAACAAGGCCac                                                                                                                                                                                                                                             >  2:275292/1‑250 (MQ=18)
                  attaCGGGTGAGAATGCCTCATTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCACCGGTAATAACGGGAAAGTGATTCAGGACGGCGAGCTGGATGTCAGCGGCGGCGGTCACGGTATTGAGAGCACCGGTGACATCGCGACGGTGGATAACAAGGGCACCAGGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATa                                                                                                                                                                                                                                     <  2:175896/250‑1 (MQ=18)
                      cGGGTGATAATGCCACAGTAAACAACAATGGAAAAACCACCGTTGACGGTAAAGATTCCACCGGTACTGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGGCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGaa                                                                                                                                                                                                                                <  2:422269/251‑1 (MQ=255)
                        ggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAgg                                                                                                                                                                                                                              <  1:360605/251‑1 (MQ=255)
                            aTAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGCGACAAGGCGGTTGTTAACAACGAAGGCg                                                                                                                                                                                                                            >  1:121258/1‑249 (MQ=255)
                                           acaacaATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAAc                                                                                                                                                                                                          <  1:112113/251‑1 (MQ=255)
                                                   ggCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcac                                                                                                                                                                                                  <  2:505525/251‑1 (MQ=255)
                                                        atccaccGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggca                                                                                                                                                                                             <  1:141205/249‑1 (MQ=255)
                                                               gTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGAGATTCAGGTCGGCGATCTGGATGTCAGCGGCGGAGGTCACGGTCTTGATATCAACGGTGACAGCGCGACGGAGGATACAAAGGGCACCCTGACCGGCACCGAACCGGAGTACATGGGTATCCAGGTCGACGGCGACAAGGCCCTCGTCAATAACCAAGGAGAAAGCACCATCCACAAAGGCggcaccggacccccga                                                                                                                                                                                       >  2:428218/1‑242 (MQ=17)
                                                                  gACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCaatg                                                                                                                                                                                                         <  1:47325/229‑3 (MQ=11)
                                                                  gACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCaatg                                                                                                                                                                                                         >  2:47325/1‑227 (MQ=11)
                                                                  gACGGTAAAGATTCCACCGGTACGGAAATCAACGGCAATAACGGGAAAGAGCTTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGGATTGATATCACCGGCGACAGCGCAACGGTGGGTAACAAGGGCACCATGACCGTCACCGATCCGGCGTCCATAGGTATCCAGATAGACGGCGGCAAGGCCCTCGTCCATACAGGATGCGATAGCACAAACACCAACGGTggcacggagcaccagagca                                                                                                                                                                                    >  2:385080/1‑236 (MQ=9)
                                                                               ccgccGGGACGTAAATCAACGGTAAGAACGGNAAAGTGATTCAGTACGGCGATCTGGATGACAGCGGCGGCGTTCACGGTAGTGATATCACCGGGGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCGCCGAGCCGGAGTCAATGGGTAACCAGATCGACGGTGACAAGGCGAGCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACTGGCACCCGGAGTGAGGGTGATGACGcg                                                                                                                                                                      <  2:249860/248‑2 (MQ=255)
                                                                                        cGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCCACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCacac                                                                                                                                                             >  1:458922/1‑248 (MQ=255)
                                                                                         ggAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGCTATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAaac                                                                                                                                                             >  1:464231/1‑250 (MQ=255)
                                                                                           aaaTCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCCGGTAGACGACGCCAATGCCATCGTCAAAAACGAGGGCGACATAACCATCCCCATCGGTGGCAACGGCACCCAGATTAATGGTGATGACGCCACGGCaacaaac                                                                                                                                                          >  1:437230/1‑246 (MQ=14)
                                                                                                  cGGTAATAACTGGAAAGTGATGCAGGACGGCGATCTGGATGTCAGCGGCGGCGGGCACGGTATGGAGATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCGCCATGACCGTCACCGTGCCGGAGTCCATGGGTAGCCAGATCGACGGTGACAAGGACATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCa                                                                                                                                                  <  2:421301/251‑1 (MQ=18)
                                                                                                    gTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGCGACAAGGCCGTCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCAGAGCCGCCCAGATCAATGGTGATGCCCCCACGCCAAACAACACCGGCaa                                                                                                                                                 >  1:12116/1‑250 (MQ=17)
                                                                                                      aataaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCCCCGATACGGAGCCCATGGGCATCCAGAACGGCGGTGACAAGGCCGTCGTCAATAACTAAGGCGAAAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAac                                                                                                                                              >  1:195794/1‑251 (MQ=18)
                                                                                                          aCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTCGACGGCGACAAGGCCGTCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAGTAAAGGGGATGACGCCACGGCAAACAACAACGGCAAAaccacc                                                                                                                                          >  2:298653/1‑251 (MQ=14)
                                                                                                          aCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAGTGTTGGCACCGGGGACAGCGCGACGGTGGATAACAAGGGCACCATGACCCTCACCTCTCCGGAGTCCATGGGTATCCAGATCGACGGTGACTAGGCCATCGTAAATAATGAAGGCGAGAGAACCATCACCAACGGTGGCCGCGGCACCCAGATTAATGGTGCTGAGGCCACGGCAAACAACACACGCAAAaccacc                                                                                                                                          >  1:100438/1‑251 (MQ=17)
                                                                                                            gggAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGGGGATAACAAGGGCACCATGACCGTCACCGATCCGGGGTCCATGGGTATCCAGGTCGACGGCGACCAGGCCGTTGGCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAAAGGGGATGACGCCCCGGCAAAAAACAACGGCAAAACCACCgg                                                                                                                                        >  2:352370/1‑250 (MQ=14)
                                                                                                                    gATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGCGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGAcaga                                                                                                                                 >  1:54535/1‑247 (MQ=255)
                                                                                                                         aGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCAACAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAt                                                                                                                           >  1:305183/1‑251 (MQ=255)
                                                                                                                            aCGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGCTCCGGAGTCCATGGGTATCCAGGTTGACGGCGACAAGGCCGTCGTTAATAACGGAGGCGAGAGCCCCATCAACAACGGTGGCACCGGCCACCAGATTAACGGGGCTGACGCCCCGGCCAACAACAACGGCCAAACCACCGTTGACGGCAAGGATTcc                                                                                                                        >  2:145140/1‑251 (MQ=11)
                                                                                                                              ggCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCa                                                                                                                                                  <  2:382060/223‑1 (MQ=255)
                                                                                                                              ggCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCa                                                                                                                                                  >  1:382060/1‑223 (MQ=255)
                                                                                                                                          gTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAt                                                                                                          >  1:112287/1‑251 (MQ=255)
                                                                                                                                             aGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGCGACAAGGCCAACGTCGAAAACGAAGGCGAGAGAACCCCCACCAACGGTGGAACCGGCACCCCGATTTATGGTTGTGACAACAAGACAAACAACAACGGCAAAACCACCGTTGACGGCAAGGAATACAACGGAACGGAAATCa                                                                                                        >  1:48800/1‑250 (MQ=17)
                                                                                                                                              gcggcggcgGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAAc                                                                                                      >  1:367027/1‑251 (MQ=255)
                                                                                                                                                      ggTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACACGGCCGTCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCCCCGGCACCCAGATTAACGGTGCTGACGCCACGGCAAACAACAACGGCACAACCACCGTTGACGGCAAGGATTCCCCCGGAACGGAAATAAACGGtaata                                                                                               >  2:280704/1‑250 (MQ=17)
                                                                                                                                                      ggTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACAGCAAAAACACCGTTGACCACAAGAATACCACCGGTACgatactcaccggtaatt                                                                                               >  1:373465/1‑234 (MQ=18)
                                                                                                                                                         cACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGAACCTCGGACCAGCGTTGTGCAGGTAGACGGTGAAGAGCCCGCTGTTCACAGCATAGGCGAAAACACCATCTCCAACGGTGGCACCGGCACCCAGAATAACGGTGATGCAGACACGGCAAACAACAACGGCCAAAACACCGGTGACGGCAAGGATTCCCCCGGTACGGCACTCACCGGTAATAACgg                                                                                           >  2:134534/1‑251 (MQ=2)
                                                                                                                                                          aCGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGGACGGAAATAAACGGTAATAACgg                                                                                           >  1:365899/1‑250 (MQ=255)
                                                                                                                                                          aCGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCAGCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGACATCAACGGTAATAACGGa                                                                                          <  2:458922/251‑1 (MQ=255)
                                                                                                                                                             gTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAACGGTAACGACGCCACGGCGAATAACAGCGGAAAAACCACCGTTGACGGAAAAGATTCCACGGGTACGGAAATCGACGGTAATATCggcatt                                                                                       >  1:332839/1‑247 (MQ=14)
                                                                                                                                                                  gATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGCGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTAt                                                                                  >  1:386487/1‑251 (MQ=255)
                                                                                                                                                                  gATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCGACGGTAATAACGGAACAGTTAt                                                                                  >  1:372542/1‑251 (MQ=255)
                                                                                                                                                                           ggTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACAGGTACGGAAATCAACGGCAATAACGGAAAAGATATTCAGGACgg                                                                         >  1:317136/1‑251 (MQ=255)
                                                                                                                                                                             tGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGCGACAAGGCCATCGTCAATAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACCACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGCAATAACGGAAAGGTTATTCCGGACGGCg                                                                       >  2:461683/1‑251 (MQ=255)
                                                                                                                                                                                cAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGcc                                                                                                                                                                      >  1:41200/1‑154 (MQ=255)
                                                                                                                                                                                cAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGcc                                                                                                                                                                      <  2:41200/154‑1 (MQ=255)
                                                                                                                                                                                cAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGACTAATGTGGATGACGCCCCGGCAAACAACAACGGCAACACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATc                                                                    >  1:50701/1‑251 (MQ=255)
                                                                                                                                                                                 aGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                                                                                                  <  1:165899/188‑1 (MQ=255)
                                                                                                                                                                                 aGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACgg                                                                                                                                  >  2:165899/1‑188 (MQ=255)
                                                                                                                                                                                 aGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACAAACCGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCTAACGGTAAGGAAATCAACAGGAATAACGGAAAGGTTGATCAGGACGGCGATCt                                                                   >  1:424802/1‑251 (MQ=18)
                                                                                                                                                                                  gcgcGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGGGGCACCGGCTGCCAGATGAACGGTGATGAGGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTg                                                                  <  1:239121/251‑1 (MQ=255)
                                                                                                                                                                                   cgcgACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCCCGGCAAACAACAACGGCAAAACCACCGTTGACGGCCAGGATTCCACCGGCACGGAAATCACCGGCAATAACGGGAAGGTGATTCAGGACGGCGAGCTgg                                                                 >  2:80852/1‑251 (MQ=14)
                                                                                                                                                                                      gACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGCGACCAGGCCATCGTCAGGAACGGCGGCGCGTAGGGCATCAGCGAAGGTGTCAGTGGCAGCCGGATTAATGGTGAAGACGACACGGCAAACAACAAGGGCAAAACCACCGTTGACGGCAAGGACTCCACCGGAACGGAAATAAAAGGTAATAACGGGAAGGTGATTGCGGACGGAGATCGGGATg                                                              >  2:416219/1‑251 (MQ=9)
                                                                                                                                                                                          gTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGAACGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAg                                                          <  1:275292/251‑1 (MQ=255)
                                                                                                                                                                                               tAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAGGTGATTCAGGACAGCgcggtggatcacaggtgtg                                                     >  1:52276/1‑233 (MQ=18)
                                                                                                                                                                                                  cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGt                                                                                                                                        >  1:331537/1‑165 (MQ=255)
                                                                                                                                                                                                  cAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGt                                                                                                                                        <  2:331537/165‑1 (MQ=255)
                                                                                                                                                                                                       gCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCAc                                             >  1:528972/1‑251 (MQ=255)
                                                                                                                                                                                                            aTGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCGATCGTCAATAACGAAGGCGAGAGCACCATCTCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTAt                                        <  1:237354/251‑1 (MQ=255)
                                                                                                                                                                                                              gACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACACAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCGACGGCAATAACGGCATAGTAATTCAGGACGGCAGCCTGGATGTCAGCGGCGGCGGGCACGGTGTTg                                      >  1:204057/1‑251 (MQ=17)
                                                                                                                                                                                                               aCCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGACGGATCCGCAGAAAACAACAACGACAAAACCACAGCGGTGGGCACGGAATGCCGCGGTACGCATGTCAACGGAAATAACGGACAAGAGATTCCGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTg                                      >  1:206769/1‑250 (MQ=16)
                                                                                                                                                                                                                    cACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGATCAAGGTATTGATATc                                 >  1:446769/1‑250 (MQ=255)
                                                                                                                                                                                                                             ggAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGCCGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTc                                                            >  2:134572/1‑214 (MQ=255)
                                                                                                                                                                                                                             ggAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTc                                                            <  1:134572/214‑1 (MQ=255)
                                                                                                                                                                                                                               aGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAgc                     >  1:372832/1‑251 (MQ=255)
                                                                                                                                                                                                                                 tCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGCGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAgcgc                   >  1:363541/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGAACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACggg             >  1:46611/1‑250 (MQ=255)
                                                                                                                                                                                                                                            aaagaccGACGGTGAGAAGGCCATCGTCAATAACGAAGGCGAGAACACGAACAACAACGGTGGCACCGGGACCCAGATTACGGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATa        <  1:3237/245‑1 (MQ=18)
                                                                                                                                                                                                                                                aTCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGCAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGAATTGATATCACCGGCGACAGCGCACAGGTGGATAACa     >  1:484671/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  cGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCCACGGTGGATAACAAgg  >  1:219828/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      
CACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGG  >  NC_000913/1466028‑1466520

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