Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0009131,594,1241.→C12.9% 13.9 0.0 ‑0.331pseudogene (2/1461 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  new base (2/2):  ref base (15/12):  total (17/14)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.61e-01

TGGGTATCGCCTTTCACGTTCATCTTATCGGTGGGCGAATTATCGCCGCCCAGCGTCGCGTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCAC  >  NC_000913/1593923‑1594330
                                                                                                                                                                                                          |                                                                                                                                                                                                              
tGGGTATCGCCTTTCACGTTCATCTTATCGGTGGGCGAATTATCGCCGCCCAGCGTCGCGTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑cccccccc                                                                                                                                                                                                        >  1:273685/1‑210 (MQ=255)
 gggTATCGCCTTTCACGTTCATCTTATCGGTGGGCGAATTATCGCCGCCCAGCGTCGCGTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCa                                                                                                                                                                                                       >  1:136111/1‑210 (MQ=255)
                             gctgggcggcgataattcgccCAGCGTCGCGTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAgtgg                                                                                                                                                                           >  2:424849/18‑210 (MQ=255)
                                               gccCAGCGTCGCGTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATca                                                                                                                                                         >  1:211173/1‑210 (MQ=255)
                                               ccccAGCGCCGCGTTGACCTGCATCAGGCCATTTTTTCCGGTATAGTTTCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATca                                                                                                                                                         <  2:211166/209‑1 (MQ=255)
                                                                     ctccggccatgtgttCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGGTTCCCCCCGTCGCTATTCATTTGAACCGGGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAAt                                                                                                                                   <  2:157861/197‑1 (MQ=255)
                                                                                gTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTg                                                                                                                        <  1:424876/210‑1 (MQ=255)
                                                                                      ggTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGAACTCACTCCGCTGTTATGTATATATATGTTTTGACTTGCAatat                                                                                                                  >  2:407557/1‑210 (MQ=255)
                                                                                            gTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCCGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATgg                                                                                                            >  1:127238/1‑210 (MQ=255)
                                                                                              tGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGATTTGCCATATTAATGGCa                                                                                                          >  2:173550/1‑210 (MQ=255)
                                                                                               gCCGTTAACGGTCAGTACCTTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCag                                                                                                         <  2:371620/210‑1 (MQ=255)
                                                                                                  gTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTTACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTCCCATATTAATGGCagatt                                                                                                      >  2:280780/1‑208 (MQ=255)
                                                                                                       cGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGaa                                                                                                 <  2:32160/210‑1 (MQ=255)
                                                                                                         gTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAAcc                                                                                               <  2:136109/210‑1 (MQ=255)
                                                                                                           cAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGa                                                                                             <  1:173556/210‑1 (MQ=255)
                                                                                                              tACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAg                                                                                          <  2:148036/210‑1 (MQ=255)
                                                                                                                   tacccGGTTTCCCCCCTTCACTGTACCTTTGAGCCAAGCGGCCAATATCCGGGTTGCCGACGATAGTAGTTTTAGCGACTCGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCAGCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAAtt                                                                                     <  2:400919/208‑1 (MQ=255)
                                                                                                                             ccccccTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCa                                                                           >  2:358869/1‑210 (MQ=255)
                                                                                                                                  cTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTTACTTGCCATATTAATGCCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACgtaa                                                                      >  2:335060/1‑207 (MQ=255)
                                                                                                                                                          cgccATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCcgcg                                              >  1:205710/1‑210 (MQ=255)
                                                                                                                                                           tacacgacgctcttccgatctgATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCcgcg                                              <  2:222012/189‑1 (MQ=255)
                                                                                                                                                                   ccAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTCCCCTGATCCGCCGTGGAATCggtg                                     >  2:246324/1‑210 (MQ=255)
                                                                                                                                                                      ggTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTa                                  <  1:427153/210‑1 (MQ=255)
                                                                                                                                                                                gATAGTGGTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGAGATCggaagagcacacgtct                         >  1:222019/1‑195 (MQ=255)
                                                                                                                                                                                   aGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAAGTTCCCCCCCCACGCTATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATcc                     >  1:221289/1‑210 (MQ=255)
                                                                                                                                                                                       gTTTTAGCGACACGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAACCGTTGTTAAGAGTCAGGATCCCgg                  >  2:462215/1‑210 (MQ=255)
                                                                                                                                                                                       gTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGTGTCAGGATCCCGGt                 >  2:120823/1‑210 (MQ=255)
                                                                                                                                                                                              cGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTCATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGt          >  2:166437/1‑210 (MQ=255)
                                                                                                                                                                                                 cacGCAGTGCCCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGttt        <  1:462247/210‑1 (MQ=255)
                                                                                                                                                                                                   cGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGCTTTGt     >  2:437751/1‑210 (MQ=255)
                                                                                                                                                                                                    gCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTc    <  1:457354/210‑1 (MQ=255)
                                                                                                                                                                                                      aGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTcac  <  1:392056/210‑1 (MQ=255)
                                                                                                                                                                                                          |                                                                                                                                                                                                              
TGGGTATCGCCTTTCACGTTCATCTTATCGGTGGGCGAATTATCGCCGCCCAGCGTCGCGTTGAACGTCATCAGGCCATTGTTTCCGGTATAGTTGCCGTTAACGGTCAGTACATTCCCCGGTTTCCCCCCTTCGCTATTCATTTGAACCGTGCCGCCATTCTCCAGGTTGCCGCCGATAGTGGTTTTAGCGACACGCAGTG‑CCCCCCCCAGGCATAACGTTCATATCACCTGCAGTGGTGACATTTCCCGACATCACACCGCTGTTATGGATATTAATGTGTTGACTTGCCATATTAATGGCAGAGTCGGAACCGAAAGCAATTTCCCCCCCCACGATATCCGTCGTACCCTGATACGCCGCGGAATCGTTGTTAAGAGTCAGGATCCCGGTACCCAGTTTTGTCAC  >  NC_000913/1593923‑1594330

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.