Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,057,8620TC22.2% 74.3 / 22.7 45intergenic (+288/+165)yeeN/asnWUPF0082 family protein/tRNA‑Asn
Reads supporting (aligned to +/- strand):  ref base T (14/21);  new base C (4/6);  total (18/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

AAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGAAATCAAAAAACCCGGTCACTTTTTTACAAGGT  >  NC_000913/2057767‑2057957
                                                                                               |                                                                                               
aaGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGt                                                                                             >  1:34040/1‑100 (MQ=255)
aaGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGCCTGt                                                                                             >  2:886014/1‑100 (MQ=255)
 aGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGtt                                                                                            <  2:1265300/100‑1 (MQ=255)
       gATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTGCATGGAACAGCGTTTTTGACCTTGCGAATCGTGATGTCTGTTGgggag                                                                                      <  2:469186/100‑1 (MQ=255)
           gCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGAc                                                                                  <  1:121023/100‑1 (MQ=255)
                 tCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGCCTGTTGGGGAGGGACAAttt                                                                             >  1:958163/1‑99 (MQ=255)
                  cccGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGCCTGTTGGGGAGGGACAATTTGc                                                                           <  1:886021/100‑1 (MQ=255)
                      cATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCAc                                                                       <  2:214878/100‑1 (MQ=255)
                      cATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCAc                                                                       <  2:657615/100‑1 (MQ=255)
                       aTAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGCCTGTTGGGGAGGGACAATTTGCTCACt                                                                      <  2:1144413/100‑1 (MQ=255)
                            gCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGc                                                                 >  2:720571/1‑100 (MQ=255)
                                 tttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTgag                                                            <  1:729345/100‑1 (MQ=255)
                                   tctcTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGAc                                                          <  2:299956/100‑1 (MQ=255)
                                      ctTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCg                                                       >  1:1009775/1‑100 (MQ=255)
                                        ttATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGCCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGAt                                                     <  1:813192/100‑1 (MQ=255)
                                        ttATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGCCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGAt                                                     <  2:291405/100‑1 (MQ=255)
                                             tttGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAgc                                                >  2:147004/1‑100 (MQ=255)
                                              ttGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAgcg                                               <  2:78669/100‑1 (MQ=255)
                                                 tACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGCCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCAc                                            <  1:400089/100‑1 (MQ=255)
                                                      tACTTCATGGAGCAGAGTTTTTGACCTTGCGAGTCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGTCTCGATTAAGCGCACGaaa                                        <  2:358306/99‑1 (MQ=255)
                                                       aCTTCATGGAACAGAGTTTTTGTTTTTTTTTTTTCTTATGTCTGTTAGGGAGGGACAATTTGCTCACTTAAGCGTGATACTCGCTTCAGCGCAc                                            >  2:459714/1‑94 (MQ=255)
                                                       aCTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGCCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAaca                                      >  1:548185/1‑100 (MQ=255)
                                                        cTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAacac                                     >  1:28919/1‑100 (MQ=255)
                                                        cTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAacac                                     >  2:546703/1‑100 (MQ=255)
                                                         ttCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAAcaca                                    <  1:1139008/100‑1 (MQ=255)
                                                          tCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAacattc                                   >  2:498811/1‑97 (MQ=255)
                                                           cATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGa                                  >  1:233249/1‑100 (MQ=255)
                                                            aTGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGaa                                 >  1:569897/1‑100 (MQ=255)
                                                             tGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGaaa                                <  1:525138/100‑1 (MQ=255)
                                                             tGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGaaa                                <  2:388423/100‑1 (MQ=255)
                                                               gAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGAAATc                              >  1:457734/1‑100 (MQ=255)
                                                               gAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGAAATc                              <  1:261225/100‑1 (MQ=255)
                                                                     agTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGAAATCaaaaaa                        <  1:1041298/100‑1 (MQ=255)
                                                                      gTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGAAATCAAAAAAc                       <  2:980824/100‑1 (MQ=255)
                                                                        ttttGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGAAATCAAAAAAccc                     <  1:979863/100‑1 (MQ=255)
                                                                        ttttGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGAAATCAAAAAAccc                     <  1:1245795/100‑1 (MQ=255)
                                                                                 tGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGAAATCAAAAAACCCGGTCACttt            >  2:77011/1‑100 (MQ=255)
                                                                                  gCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGAAATCAAAAAACCCGGTCACtttt           <  2:507667/100‑1 (MQ=255)
                                                                                  gCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGAAATCAAAAAACCCGGTCACtttt           >  2:1022856/1‑100 (MQ=255)
                                                                                   ngAATCGTGATGCCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGAAATCAAAAAACCCGGTCACttttt          >  2:362858/2‑100 (MQ=255)
                                                                                         gTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGAAATCAAAAAACCCGGTCACTTTTTTACAAg    <  2:34039/100‑1 (MQ=255)
                                                                                         gTGATGTCGGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGAAATCAAAAAACCCGGTCACTTTTTTACAAg    <  2:818552/100‑1 (MQ=255)
                                                                                         gTGATGCCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGAAATCAAAAAACCCGGTCACTTTTTTACAAg    <  2:933608/100‑1 (MQ=255)
                                                                                           gATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGAAATCAAAAAACCCGGTCACTTTTTTACAAGGt  >  2:16798/1‑100 (MQ=255)
                                                                                           gATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGAAATCAAAAAACCCGGTCACTTTTTTACAAGGt  <  1:981600/100‑1 (MQ=255)
                                                                                               |                                                                                               
AAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAATTTGCTCACTGAAGCGTGAGACTCGATTAAGCGCACGAAACACAGAAATCAAAAAACCCGGTCACTTTTTTACAAGGT  >  NC_000913/2057767‑2057957

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