Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,222,5960GA41.2% 32.4 / 13.4 36F13F (TTC→TTTymgDuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (14/6);  new base A (0/14);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 4.05e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAG  >  NC_000913/1222448‑1222741
                                                                                                                                                     |                                                                                                                                                  
cAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGa                                                                                                                                                   >  1:89797/1‑150 (MQ=255)
        gCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAaaaca                                                                                                                                             <  2:204522/147‑1 (MQ=255)
        gCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAaacaaaca                                                                                                                                         <  1:313192/151‑1 (MQ=255)
         cattacccattaccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAg                                                                                                                                        <  1:609930/151‑1 (MQ=255)
                     ccGCAACATTTTCTGCAACCTGTGCGCTGTCAGTTTTTCCCTTTTCGTTGGTTTCTGCGAAAGTAATATGGGTTTCTGACGGAAAACCTTTTTTAATATTGTAGCCCTGCCCCCAGGTCATTGGTGCGAAAAAAAACAGGCCCGCCAGTa                                                                                                                             <  2:266490/150‑1 (MQ=255)
                     ccGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTa                                                                                                                             >  2:619754/1‑150 (MQ=255)
                      cGCAACATTTTCTGCACCCTGTGCTCTGTCATTTTTTGCCATTTCGTTGGTTTCTGCGCAAGTAATAGAGGTTTCTGACGGAAAACCGTTTTTAATTTTGTAGTCCTGCGCCCAGGTCATTGGGGCAAAAAAAAACAGGCCCGCCAGTaaa                                                                                                                           <  2:378459/151‑1 (MQ=255)
                       gCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTaa                                                                                                                            >  1:583530/1‑149 (MQ=255)
                            attttCTGCACCCTGTGCGCTGTCATTTTTTGCCATTTCGTTGGCTTTGGCGCAAGTAATATAGTTTTCTGACGGCAACCCTTTTTAAATATTGTAGTCCGGCGCCCGGTTCATTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGc                                                                                                                         <  2:353028/145‑1 (MQ=255)
                              tttCTGCACCCTGTGCGTTGTCATTTTTTGCCTTTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGAAACCCTTTTTTAATATTGTAGTCCTGCGCCCAGTTCATTGGTGAAAAAAAAAACAGGCCCGCCAGTAAAGCAAAtt                                                                                                                    <  1:485278/150‑1 (MQ=255)
                                     aaCCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatc                                                                                                            >  1:253165/1‑151 (MQ=255)
                                          gtgCGCTGTCATTTTTTGCCTTTTCTTTGGTTTCTGCGCAATTAATATAGTTTTCTGACGGAAAACCTTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatcatc                                                                                                         <  1:640970/149‑1 (MQ=255)
                                          gtgCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTcatcatcat                                                                                                       >  2:552112/1‑151 (MQ=255)
                                            gcgtTGTCTTTTTTGGCCTTTTCTTTGGTTTCTGCGCAAGTAATATAGGTTTCTGACGGCAACCCGTTTTTAATATTGTAGCCCTGCCCCCAGTTTATTGGTCCAAAAACAACAAGGCCCGCCAGTAAAGCAAATTTTTtcatcatca                                                                                                        <  2:455594/148‑1 (MQ=255)
                                              gcTGTCATTTTTTGCCTTTTCTTTTGTTTCTGCTAAAGTAATATGGTTTTCTGACGGCAAACCTTTTTAAATATTGTAGCCCTGCCCCCAGGTATTTGGTGAAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTtcatcatca                                                                                                        <  2:630516/146‑1 (MQ=255)
                                                   ctttttttGCCTTTTCGTTGGCTTCTGCGCAATTAATATAGTTTTCTGACGGAAAACCTTTTTTAATATTGTAGTCCTGCGCCCAGTTCTTTGCTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCt                                                                                                   <  2:56120/143‑1 (MQ=255)
                                                      tttttGCCTTTTCGTTGTTTTCTGCGCAAGTAATATAGTTTTCTGACGGAAACCCTTTTTAAATATTGTAGTCCTGCGCCCAGGTCTTTGGTGCAAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTAtt                                                                                               <  1:576442/147‑1 (MQ=255)
                                                               tttCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTAccc                                                                                    >  1:572928/1‑149 (MQ=255)
                                                                   gTTGGCTTCTGCGCAAGTAATATAGGTTTCTGCCGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAAtt                                                                              >  2:140540/1‑151 (MQ=255)
                                                                         ttCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGcacc                                                                         >  1:148184/1‑150 (MQ=255)
                                                                                    tAATATAGTTTTCTGCCGCCAACCTTTTTTTATTATTGTAGCCCCGCCCCCAGTTCATTGGGGCGAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGc                                                              <  1:95597/150‑1 (MQ=255)
                                                                                        ataGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATc                                                         >  1:82380/1‑151 (MQ=255)
                                                                                                  gACGGCACACCGTTTTTAATATTGTAGTCCTGCGCCCAGGACATTGGTGGGAAAACCAACAGACCCGCCAGTAAAGCATATTTTTTCATCATCATCCCTTATTTCTTTTTACCCAGACTTGCACCTCCCGTACCGCCCATCACTGCCCCtt                                               >  2:388304/1‑151 (MQ=255)
                                                                                                     cgCAAACCCTTTTTTAATTTTTTAGTCCTCCCCCCAGTTCTTTGGTGAAAAAACAAACGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCttt                                              <  2:67857/149‑1 (MQ=255)
                                                                                                                 tttAATATTGTAGCCCGGCCCCCAGTGCTTTGTCGAAAAAACAAACGGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTc                                  <  1:401638/150‑1 (MQ=255)
                                                                                                                        ttgtgttCCTGCCCCAAGGTATTGGGTGCAAAAAAAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCg                                 <  1:635048/138‑1 (MQ=255)
                                                                                                                        ttGTATCCCTGCCCCAAGGTATTGGGTGAAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCa                           <  2:89780/150‑1 (MQ=255)
                                                                                                                        ttGTAGCCCTGCCCCAAGTTCATTGGTGCAAAAAAAAACAGCCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGa                                <  1:672875/145‑1 (MQ=255)
                                                                                                                        ttGTAGCCCGGCCCCCGGGTCTTTGGTGCAAAAAAAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGcc                            <  1:71366/149‑1 (MQ=255)
                                                                                                                         tgtatgcccggcccccGGTCCTTGGTGCCAAAACCAACCGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCa                           <  2:395159/138‑1 (MQ=255)
                                                                                                                                cTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTc                  >  1:644009/1‑151 (MQ=255)
                                                                                                                                     cccAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAAc             >  2:149635/1‑151 (MQ=255)
                                                                                                                                          gTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGcc          >  2:411515/1‑149 (MQ=255)
                                                                                                                                             aTTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTg      >  1:527680/1‑150 (MQ=255)
                                                                                                                                              ttGGTGCAAAAAAAAACGGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGt       <  1:113110/148‑1 (MQ=255)
                                                                                                                                               aggggagAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCATATTTCATTTTACCCAGAATTGCACCACCCTTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAGCGCCAg        <  2:318022/140‑1 (MQ=255)
                                                                                                                                                  tGCTAAAACAAACAGGCCAGCCAGTAAAGCAATTTTTTTCATCAACATTACTTATTACATTTTACCCAGAATTGCACCACACGTACCGCCAATCATGGCACCTTTAATCGCTCCATAGATTCTGTTGCCGGTCCGAACCCCAGTGACAg  >  2:360372/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CAACGCTGGCATTACCCATTACCGCAACAATTTCTGCAACCTGTGCGCTGTCAGTTTTTGCCATTTCGTTGGCTTCTGCGCAAGTAATATAGGTTTCTGACGGCAAACCGTTTTTAATATTGTAGTCCTGCGCCCAGGTCATTGGTGCGAAAACAAACAGGCCCGCCAGTAAAGCAAATTTTTTCATCATCATTCCTTATTTCATTTTACCCAGAATTGCACCACCCGTACCGCCAATCACGGCACCTTTAATCGCCCCTTCGAGGCCATTGCCGGTCAGAACGCCAGTGACAG  >  NC_000913/1222448‑1222741

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