Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,466,138 G→A 9.3% pseudogene (747/2513 nt) ydbA → pseudogene, autotransporter homolog, interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,466,1380GA 9.3% 218.4 / 6.5 71pseudogene (747/2513 nt)ydbApseudogene, autotransporter homolog, interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (34/30);  new base A (2/5);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

AACTATTCGCGATTCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTA  >  NC_000913/1465856‑1466409
                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                               
aaCTATTCGCGATTCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAg                                                                                                                                                                                                                                                                             <  1:331765/287‑1 (MQ=255)
aaCTATTCGCGATTCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAg                                                                                                                                                                                                                                                                             >  2:331765/1‑287 (MQ=255)
             tCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCg                                                                                                                                                                                                                                                              >  2:404522/1‑289 (MQ=255)
             tCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCg                                                                                                                                                                                                                                                              <  2:540328/289‑1 (MQ=255)
                  ttttACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGAt                                                                                                                                                                                                                                                            <  2:860615/286‑1 (MQ=255)
                           taCCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACAGGAAAGTGATGCAGGACGGCgcacaggttgccacc                                                                                                                                                                                                                                                >  1:640462/1‑275 (MQ=255)
                           taCCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCGTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGGAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAATGATTCAGGACGGGGATCTGGATGTCAgc                                                                                                                                                                                                                                                >  2:140767/1‑289 (MQ=255)
                                agaATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAgcggcg                                                                                                                                                                                                                                            >  2:107193/1‑288 (MQ=255)
                                       tGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCa                                                                                                                                                                                                                                     <  2:644521/288‑1 (MQ=255)
                                       tGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCa                                                                                                                                                                                                                                     <  2:700409/288‑1 (MQ=255)
                                       tGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCa                                                                                                                                                                                                                                     <  1:606352/288‑1 (MQ=255)
                                                       ggggtaaaaggggcagaaagTAAGGCAACGATTAATATTGGGCAGATGGATGAAGGGAATAACACTGTTGCCCATGAAGGGGTGAGCGCAGATGGCGCAAAGAAGAGGAAATATAATCACAACGGTGGGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAcc                                                                                                                                                                                                                     <  1:543412/269‑1 (MQ=21)
                                                        tGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg                                                                                                                                                                                                                   >  2:211331/1‑289 (MQ=255)
                                                        tGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg                                                                                                                                                                                                                   >  1:871719/1‑289 (MQ=255)
                                                        tGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCgg                                                                                                                                                                                                                   >  1:82231/1‑289 (MQ=255)
                                                          cAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTg                                                                                                                                                                                                                 >  2:492470/1‑289 (MQ=255)
                                                          cAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTg                                                                                                                                                                                                                 >  2:375673/1‑289 (MQ=255)
                                                            aGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGAc                                                                                                                                                                                                               >  2:671909/1‑289 (MQ=255)
                                                                aGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAgc                                                                                                                                                                                                            >  2:614815/1‑288 (MQ=255)
                                                                   aaTGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGa                                                                                                                                                                                                        <  2:871719/289‑1 (MQ=255)
                                                                       gTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGt                                                                                                                                                                                                    >  2:592785/1‑289 (MQ=255)
                                                                            aaGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGAt                                                                                                                                                                                                <  1:671909/288‑1 (MQ=255)
                                                                                 catgatgaaagcgtGGCAGAATGCCGAAGCGAATAAGACGGTTGCCCGTGAAGGCGTGGGCGAAGAAGGCGCAAAGAAGTGGAAAACTAATCACAACGGCGGGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACa                                                                                                                                                                                            <  1:29097/274‑1 (MQ=17)
                                                                                 aaCGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAg                                                                                                                                                                                          <  1:614815/289‑1 (MQ=255)
                                                                                 aaCGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAg                                                                                                                                                                                          <  1:404522/289‑1 (MQ=255)
                                                                                 aaCGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACAGGTGACAGCGCGACGGTGGATAACAAg                                                                                                                                                                                          >  2:538538/1‑289 (MQ=255)
                                                                                  aCGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAgg                                                                                                                                                                                         >  2:543412/1‑289 (MQ=255)
                                                                                              tGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCg                                                                                                                                                                             >  1:672334/1‑289 (MQ=255)
                                                                                               ggCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCg                                                                                                                                                                             >  2:29097/1‑288 (MQ=255)
                                                                                                    aTTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGGGCTTGTTATTACGGGTGATAATACCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAATGATTCCACCGGTACGGAAATTGCTGGCAATAACGGGAAGGTGATTCAGGACGGCGATCTGGATGCCAGCGGCCGCGGGCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACTATGACCGTccacg                                                                                                                                                                       >  1:216509/1‑285 (MQ=18)
                                                                                                         tgaAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                                                                                                                                                  <  2:82231/289‑1 (MQ=255)
                                                                                                         tgaAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCg                                                                                                                                                                  <  1:375673/289‑1 (MQ=255)
                                                                                                              cGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGt                                                                                                                                                              >  1:799632/1‑288 (MQ=255)
                                                                                                               gAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTc                                                                                                                                                             <  2:445825/288‑1 (MQ=255)
                                                                                                                aaTAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAg                                                                                                                                                                                                             <  1:318748/239‑1 (MQ=255)
                                                                                                                aaTAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAg                                                                                                                                                                                                             >  2:318748/1‑239 (MQ=255)
                                                                                                                   aaCACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATgg                                                                                                                                                        <  1:548565/289‑1 (MQ=255)
                                                                                                                        tGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTAt                                                                                                                                                    >  2:586312/1‑288 (MQ=255)
                                                                                                                               cTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCg                                                                                                                                            >  1:241548/1‑289 (MQ=255)
                                                                                                                                ttGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCg                                                                                                                                            >  2:31762/1‑288 (MQ=255)
                                                                                                                                ttGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGa                                                                                                                                           <  1:107193/289‑1 (MQ=255)
                                                                                                                                 tGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGAc                                                                                                                                          <  2:241548/289‑1 (MQ=255)
                                                                                                                                  gAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACg                                                                                                                                         >  2:813537/1‑289 (MQ=255)
                                                                                                                                  gAAGGGGGGAGCGCAGATGGGGAAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGAc                                                                                                                                          <  1:538538/288‑1 (MQ=255)
                                                                                                                                        gTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACa                                                                                                                                   >  2:167632/1‑289 (MQ=255)
                                                                                                                                                aGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGCCGGTGCCAAGGCCATc                                                                                                                           >  2:662717/1‑289 (MQ=255)
                                                                                                                                                  aTGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGAGATAAGGCCATCg                                                                                                                          >  2:844894/1‑288 (MQ=255)
                                                                                                                                                                   caaacaaatgaaaacgacgaaccacgtattaAGGGGGATAAGGCCACAGTAAACAACAAAGACAAAACCACCGTTGAGGGTAAAGATTCCGCAGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCgaga                                                                                                        <  1:813537/263‑1 (MQ=14)
                                                                                                                                                                         aaTCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCa                                                                                                  <  2:799632/289‑1 (MQ=255)
                                                                                                                                                                                      aGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACg                                                                                                                                                                                                      <  1:749581/176‑1 (MQ=255)
                                                                                                                                                                                      aGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACg                                                                                                                                                                                                      >  2:749581/1‑176 (MQ=255)
                                                                                                                                                                                         ttGTTATTACGGGTGATAATGCCACAGTAAACAACAAGGGCAAACCCACCGTGGACGGTAAAGATCCACCCGGTACGGAATTCAACGGTAATAACGGAAAGGGGATTCAGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTATCCAGATTGACGGCGACCAGGCCATCGTCAATAACGAAGGCGAGAGCACTATCACCAATGGCggcac                                                                                  <  2:618556/288‑1 (MQ=11)
                                                                                                                                                                                                taagggTGAAAAAGCCAAAAGAAACCACAATGGCAAAACCCACGGTGACGGTAAAGATTCCACCGGTAACGAAATCAACGGTAAAAACGGGAAAGTGATTCAGGGAGGCGGTCTGGATGTCAGGGGCGGCGGTCAAGGTATTGATAACACCGGTGACAGCGCGACGGGGGATAACAAGGGCAACATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcacc                                                                           <  1:592785/286‑1 (MQ=17)
                                                                                                                                                                                                taCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcacc                                                                           <  1:492470/289‑1 (MQ=255)
                                                                                                                                                                                                taCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcacc                                                                           <  1:31762/289‑1 (MQ=255)
                                                                                                                                                                                                taCGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACAGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcac                                                                            >  2:315892/1‑288 (MQ=255)
                                                                                                                                                                                                     gTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGAt                                                                      >  2:825636/1‑289 (MQ=255)
                                                                                                                                                                                                                              tGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGAt                                                                                                                                                                                                >  2:136416/1‑142 (MQ=255)
                                                                                                                                                                                                                              tGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAgcggcg                                                                                                                                                                                                                                            >  1:82111/1‑98 (MQ=9)
                                                                                                                                                                                                                              tGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGAt                                                                                                                                                                                                <  1:136416/142‑1 (MQ=2)
                                                                                                                                                                                                                                gCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcac                                                                                  >  1:136062/1‑250 (MQ=255)
                                                                                                                                                                                                                                gCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcac                                                                                  <  2:136062/250‑1 (MQ=255)
                                                                                                                                                                                                                                       aaccGGTAACGGTAAAGATTCCACCGGTAAGGAAAGCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATGGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGc                                    <  1:211331/288‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggcaaggATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTCTCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaa                                  <  2:249319/275‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggcaaggATTCCACCGGTACGGAAATCAACGGTAATAACGGAAAAGTTATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGCGACAGCGCAACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATCGGTCTCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCaa                                  >  1:249319/7‑281 (MQ=255)
                                                                                                                                                                                                                                                         ttCCTCCGGTACGGAAATCAACGGTAAGAACAGGAAAGTGATTCAGGCCGGCGATCTGGAGGGCAGCGGCGGCGGTAACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACGAACGGGGGCACCGGCACTCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGc                  <  1:315892/289‑1 (MQ=17)
                                                                                                                                                                                                                                                                         aTCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTTAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTa  <  2:672334/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aTCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTa  <  2:640462/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                            aaCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcacc                                                                           >  2:789104/1‑213 (MQ=255)
                                                                                                                                                                                                                                                                            aaCGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTggcaccggcacc                                                                           <  1:789104/213‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   ataaCGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCAcc                                                                                                                                                                                                                     >  1:541637/1‑68 (MQ=35)
                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                               
AACTATTCGCGATTCAGTTTTTACTTATACCGAGAATGCTGACGGGACTATATCTCTGCAAGATAGCAATGGTCGTAAGGCAACGATTAATCTTTGGCAGATTGATGAAGCGAATAACACTGTTGCCCTTGAAGGGGTGAGCGCAGATGGCGCAACGAAGTGGCAATATAATCACAACGGTGAGCTTGTTATTACGGGTGATAATGCCACAGTAAACAACAATGGCAAAACCACCGTTGACGGTAAAGATTCCACCGGTACGGAAATCAACGGTAATAACGGGAAAGTGATTCAGGACGGCGATCTGGATGTCAGCGGCGGCGGTCACGGTATTGATATCACCGGTGACAGCGCGACGGTGGATAACAAGGGCACCATGACCGTCACCGATCCGGAGTCCATGGGTATCCAGATCGACGGTGACAAGGCCATCGTCAATAACGAAGGCGAGAGCACCATCACCAACGGTGGCACCGGCACCCAGATTAATGGTGATGACGCCACGGCAAACAACAACGGCAAAACCACCGTTGACGGCAAGGATTCCACCGGTA  >  NC_000913/1465856‑1466409

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