Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,4320AG25.0% 64.5 / 21.3 44pseudogene (1041/2518 nt)ydbAputative outer membrane protein, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base A (21/12);  new base G (5/6);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAG‑AGC‑‑‑ACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGC  >  NC_000913/1466150‑1466713
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                         
ggACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATTCTAAATCTCTAGGTATCCATATTCAAGTTGCCCCGGCCATCGTTCCTAACGAATGCGAG‑GTT‑‑‑ACCCTGCCCACTGGTGGCCCCGGCACGCAGATTAATGGTGATGACGCCGCGGCAAACCACAACGGCAAAACCACCGTTGACGGCAAGGATTACACCGGTACGGAAATCAACGGTAATAAAGGAAAGGTTATTCAGGACgg                                                                                                                                                                                                                                                                          >  2:82511/1‑300 (MQ=12)
        aTCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAAC‑ACGTATATTATC‑‑‑CCCATCCCCAAGCCTAAAATCCTCACTCGCGTCTATGGTGGTGACGCCACGGCAAACAACAACGGAAAAACAAAGGTTGAAGGAAAGGATTCCAACGGTAAGGAAAGCAACGGTAATAAAGGAAAAGTTATTCAGGACGGCGATCTgg                                                                                                                                                                                                                                                                  >  1:691540/1‑300 (MQ=2)
                                aCGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAG‑AGC‑‑‑ACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTa                                                                                                                                                                                                                                          <  1:351065/300‑1 (MQ=14)
                                        gATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAG‑AGCCCGACCACCACCAACA‑T‑‑CCCCGACCCCCCGATCCCCATTAACGCCTCCAAGGCAAACACCAACGGCAAAACAACCGTTGAAGGAAAGGATTCCAACGGTACGGAAATTGATGGTAATAATGGGAAGGTTATTCAGGACGGCGATCTGGATGTAAGCGGCGGCGGTCACGGTATTGATATc                                                                                                                                                                                                                                  >  1:611121/1‑300 (MQ=9)
                                                           cgACGGTGGATAACAAGTGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAG‑AGC‑‑‑ACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACgg                                                                                                                                                                                                               <  1:81123/300‑1 (MQ=255)
                                                                  ggATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAG‑AGC‑‑‑ACCATCAACAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCCCAGGTACGGAAATTGCTGGCAATAACGGGGAGGTGATTCAGGACGGCGATCTGGATGNCAGCGGCGGCGGTCACGGTATTGATATCAACGGCGACAGCGCAACGGTGGATaa                                                                                                                                                                                                        >  2:505374/1‑300 (MQ=11)
                                                                   gATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAG‑AGC‑‑‑ACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAAc                                                                                                                                                                                                       <  1:335830/300‑1 (MQ=255)
                                                                                 cATGACCGCCACTGATCCGGAGTCCAGGGCTAGCCAGATCGACGTCGACAAATCCATTGTCAACAACGAAGGCGAG‑AAC‑‑‑ACCAGCGATAACCGTGACACCGGCACCCAGATTAATGTTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGAc                                                                                                                                                                                         <  2:367529/300‑1 (MQ=17)
                                                                                            ccGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAG‑AGC‑‑‑ACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTAACGACGCCACTGCGAACAACAACGGAAAAACCACTGTTGATGGAAAGGATTCCACTGGTACGGAAATCAAGGGTAATACCGGAAAGGTTATTCAGGACGGTGATCTGGATGTCAGCGGCGGTGGTCATGGTATTGATATCACTGGCGACAACGCCACGGTGAATAACacgtacaccatacccatccccgatc                                                                                                                                                                              >  1:75513/1‑276 (MQ=255)
                                                                                              gATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAG‑AGC‑‑‑ACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACANCAACGGCAAAACCACCGTTGAAGGCANGGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAGGTGATTCAGGAAGTCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAACGGCGACAGCGCAACGGTGGATAACAAGGGAACCATGAACGTCAACGATCCg                                                                                                                                                                            >  2:465601/1‑300 (MQ=17)
                                                                                               aTCCGGAGTCCATGGGGATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAG‑AGC‑‑‑ACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCgg                                                                                                                                                                           <  1:210233/300‑1 (MQ=255)
                                                                                                      gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAG‑AGC‑‑‑ACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAATAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACAGGTACGGAAATCAACGGTAATAACGGAAGGGTTATTCAGGATGGTGGTCTGGGTGTCGGTGGCGGTGGTCACGGTATTGGTATCATCGGCGACAGCGCAACGCTTGATAACAAGAGCACCAGGAGCGTCACTTATCCGGAGTACAt                                                                                                                                                                    >  2:16542/1‑300 (MQ=18)
                                                                                                      gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAG‑AGC‑‑‑ACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCAt                                                                                                                                                                    >  1:566077/1‑300 (MQ=255)
                                                                                                      gTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAG‑AGC‑‑‑ACCATCAACAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAACGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCAGAGTCCAt                                                                                                                                                                    >  2:64213/1‑300 (MQ=255)
                                                                                                            gggTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAG‑AGC‑‑‑ACCATCACCAACGGTGGCACCGGCACCCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAAGATTCCACCGGTACAGAAATCAAAGGTAATAGCGGAAACGTTCTTCAGACCCATGCTCTGGCTGCCGGCGGCGGTGATCGCGGTATTGATATCACCGGCGACAGTGCAACGGTGGATAACAAGGGCACAATGAACGTCACCGATCCGGAGTCCATCGGTAt                                                                                                                                                              >  1:155811/1‑300 (MQ=17)
                                                                                                                aTCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAG‑AGC‑‑‑ACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGATGGCAAGGATTCCACCGGTACGGAAATCGCTTGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAg                                                                                                                                                          >  1:111694/1‑300 (MQ=11)
                                                                                                                   cAGATCGACGGTTATGAGGTCATCGTCACTGACGAAGTCAAG‑AGC‑‑‑ATCAGAAACAACCGCGGCACCGGCATCCGTATGACTGGTGATAACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGGCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGtt                                                                                                                                                       <  2:674047/300‑1 (MQ=14)
                                                                                                                        cGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAG‑AGC‑‑‑ACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACgg                                                                                                                                                  <  1:563685/300‑1 (MQ=255)
                                                                                                                                          cGTCAATAACGAAGGCGAG‑AGC‑‑‑ACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAGCGGTGGATAACAAGGGCACCATGACCGNCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGccatcgtcaa                                                                                                                                >  2:502115/1‑291 (MQ=255)
                                                                                                                                              aagaaCGAAGGCGAG‑AGC‑‑‑ACTATCACCAACGGTGGCACCGGTACTTAGATTAATGTTGATGACGCCACGGCAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaacaac                                                                                                                            <  1:13270/297‑1 (MQ=255)
                                                                                                                                              aaTAACGAAGGCGAG‑AGC‑‑‑ACCATCACCAACGGTGGCACCGGCACCCAGTTTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACNGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGNCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTaacaac                                                                                                                            <  2:554254/300‑1 (MQ=255)
                                                                                                                                                   cGAAGGCGAG‑AGC‑‑‑ACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAACGGCGACAGCGCAACGGTGGATCACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGaataggcaatcgtcaatcacgaagg                                                                                                                       >  1:290685/1‑276 (MQ=14)
                                                                                                                                                       ggCGAG‑AGC‑‑‑ACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGAACGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCAGTTGTTAACAACGAAGGCgag                                                                                                                   >  1:84171/1‑300 (MQ=255)
                                                                                                                                                         cGAG‑AGC‑‑‑ACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCAATCGTCAACAACGAAGGCGagag                                                                                                                 >  1:553529/1‑300 (MQ=255)
                                                                                                                                                         cGAG‑AGC‑‑‑ACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGAACAGGCGGTTGTTAACAACGAAGGCGagag                                                                                                                 >  1:553513/1‑300 (MQ=255)
                                                                                                                                                                    aCCATCAACAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCAGTCGTCAACAACGAAGGCGAGAGCACCATc                                                                                                          >  1:36971/1‑300 (MQ=255)
                                                                                                                                                                      cATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCAc                                                                                                        <  1:512374/300‑1 (MQ=11)
                                                                                                                                                                       agcactaaCGGTGGCATCGACACCCAGATTAATGGTGATAACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGTAAGGATTCCACCGGTACGGAAATCAACGGTAATAACTGAAAAGTTATTCAGNACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCAcc                                                                                                       <  2:164883/294‑1 (MQ=255)
                                                                                                                                                                          aCCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCaat                                                                                                    >  2:349916/1‑299 (MQ=11)
                                                                                                                                                                               cGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTgg                                                                                               >  2:319416/1‑300 (MQ=255)
                                                                                                                                                                                                 gATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCCCCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGCCGGCGATCTGGCTGTCAGCGGCTGCGGTCACGGTATTGATATCACAGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCAAGGTTGGCGGCGACCAGGTGGTTGTTAACAACGAAGGAGAGAGTGCTATCACCACTGGCGGCAACGGCATTCAGATTaa                                                                             >  2:341096/1‑300 (MQ=17)
                                                                                                                                                                                                 gATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATCGACGGCGACCCGGCGATTGTCAATAACGAAGGCGAGAGCACGATCACCAACAGCGGCACCGGCACTTAGATTaa                                                                             >  2:50560/1‑300 (MQ=11)
                                                                                                                                                                                                    tacgggTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGGACGGAAATCAACGTTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACgg                                                                          <  2:368471/296‑1 (MQ=11)
                                                                                                                                                                                                    tAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACgg                                                                          <  1:60797/300‑1 (MQ=255)
                                                                                                                                                                                                    tAATGGTGATAATGCCACGGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGTAAATCAACGGTAATAACGGGAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTCGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACgg                                                                          <  2:368711/300‑1 (MQ=11)
                                                                                                                                                                                                    agatgtgtataagagacaggtaaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACgg                                                                          <  2:388441/279‑1 (MQ=11)
                                                                                                                                                                                                                       gtaaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACggctgtctcttatacacatct                                                       >  1:394245/3‑281 (MQ=11)
                                                                                                                                                                                                                       gtaaaCAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaac                                                       <  1:246460/298‑1 (MQ=255)
                                                                                                                                                                                                                       gCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAGTAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCAACAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAaac                                                       >  1:147711/1‑300 (MQ=255)
                                                                                                                                                                                                                       gCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCAATCGTCAACAACGAAGGCGAGAGCACCATCACCAACGGCGGCACCGGCACTCAGATCAACGGTAACGACGCCAcagcgaat                                                       >  1:147581/1‑293 (MQ=11)
                                                                                                                                                                                                                          aacaacaatGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGTTGATGACGCCACGGCAaacaac                                                    <  1:652793/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                        cGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGTTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATAACCGACGGGGGCGCCGGTACGCGGAGGAATGTGGATGACGCCGCGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGaaa      <  1:439413/300‑1 (MQ=17)
                                                                                                                                                                                                                                                                          gAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAAtt    >  1:240611/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                            aaTCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCNATCCGGAGTCCATCGGTATCCAGGTTNACGGCGACCAGGCAGTTGTTAACAACGAAGGCGAGAGCACCATCAACAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACAGGTACGAAAAtcga  >  2:455082/1‑297 (MQ=255)
                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                         
GGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAG‑AGC‑‑‑ACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGGTTGACGGCGACCAGGCGGTTGTTAACAACGAAGGCGAGAGCGCCATCACCAACGGTGGCACCGGCACGCAGATTAACGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATTGC  >  NC_000913/1466150‑1466713

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