Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,817,1320TC49.2% 2.2 / 85.2 61G124G (GGA→GGGyqaBfructose‑1‑phosphatase
Reads supporting (aligned to +/- strand):  ref base T (17/14);  new base C (18/12);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 7.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAAC  >  NC_000913/2817055‑2817211
                                                                             |                                                                               
aacaaAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAAcccc                                                                              <  1:507807/81‑1 (MQ=255)
aacaaAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATcc                                                                               >  2:507857/1‑80 (MQ=255)
    aaaaTGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGt                                                                          >  2:507868/1‑81 (MQ=255)
    aaaaTGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGt                                                                          <  2:507808/81‑1 (MQ=255)
      aaTGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTgc                                                                        <  2:507809/81‑1 (MQ=255)
       aTGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTgcg                                                                       <  1:507810/81‑1 (MQ=255)
        tgtgTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTgcgc                                                                      <  1:507811/81‑1 (MQ=255)
         gtgtCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTgcgc                                                                      >  2:507906/1‑80 (MQ=255)
          tgtCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGcc                                                                     <  1:507813/80‑1 (MQ=255)
          tgtCTGGCGCGGGTTTATGGTGTTTGACGTAATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGcc                                                                     <  2:507812/80‑1 (MQ=255)
          tgtCTAGCGCGAGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGcc                                                                     <  2:507814/80‑1 (MQ=255)
           gtCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAg                                                                   >  2:507899/1‑81 (MQ=255)
            tCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAg                                                                   <  2:507815/80‑1 (MQ=255)
                gcgcGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATg                                                              <  2:507816/81‑1 (MQ=255)
                  gcgGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGcc                                                            <  2:507817/81‑1 (MQ=255)
                         tATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCg                                                      <  2:507818/80‑1 (MQ=255)
                         tATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCg                                                      <  2:507819/80‑1 (MQ=255)
                         tATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGa                                                     >  1:507872/1‑81 (MQ=255)
                          aTGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGAt                                                    >  1:507877/1‑81 (MQ=255)
                          aTGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGAt                                                    >  2:507870/1‑81 (MQ=255)
                          aTGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGAt                                                    >  1:507887/1‑81 (MQ=255)
                           tGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGAt                                                    <  1:507820/80‑1 (MQ=255)
                             gtgtTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATgg                                                  <  2:507822/80‑1 (MQ=255)
                             gtgtTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGgc                                                 <  2:507821/81‑1 (MQ=255)
                             gtgtTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGgc                                                 >  2:507883/1‑81 (MQ=255)
                              tgtTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGgcg                                                <  2:507823/81‑1 (MQ=255)
                               gtTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGgcgc                                               >  1:507879/1‑81 (MQ=255)
                               gtTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGgcgc                                               >  1:507896/1‑81 (MQ=255)
                                tttGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGgcgc                                               <  2:507824/80‑1 (MQ=255)
                                tttGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGgcgc                                               <  1:507825/80‑1 (MQ=255)
                                 ttGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCtt                                             >  2:507876/1‑81 (MQ=255)
                                  tGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCttt                                            >  2:507882/1‑81 (MQ=255)
                                   gACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCttt                                            <  2:507826/80‑1 (MQ=255)
                                        gATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTcc                                      >  1:507900/1‑81 (MQ=255)
                                         aTCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTccc                                     >  1:507908/1‑81 (MQ=255)
                                         aTCGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTccc                                     >  1:507892/1‑81 (MQ=255)
                                           cGGCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCg                                   >  1:507907/1‑81 (MQ=255)
                                            ggCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGt                                  >  1:507912/1‑81 (MQ=255)
                                            ggCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGt                                  >  2:507905/1‑81 (MQ=255)
                                            ggCAGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGt                                  <  1:507827/81‑1 (MQ=255)
                                             gCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGGGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGtt                                 <  1:507828/81‑1 (MQ=255)
                                              cAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTc                                <  2:507829/81‑1 (MQ=255)
                                               aGCGACGACGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTcc                               >  1:507904/1‑81 (MQ=255)
                                                 cgacgaCGGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTa                             >  2:507861/1‑81 (MQ=255)
                                                  gacgacGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTa                             >  1:507909/1‑80 (MQ=255)
                                                     gacgGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGc                         >  2:507888/1‑81 (MQ=255)
                                                      acgGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGc                         >  2:507875/1‑80 (MQ=255)
                                                      acgGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGcc                        <  2:507830/81‑1 (MQ=255)
                                                       cgGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCa                       >  2:507915/1‑81 (MQ=255)
                                                       cgGCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCa                       >  1:507894/1‑81 (MQ=255)
                                                        ggCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCAt                      >  2:507874/1‑81 (MQ=255)
                                                        ggCGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCAt                      >  1:507895/1‑81 (MQ=255)
                                                          cGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTg                    <  2:507831/81‑1 (MQ=255)
                                                          cGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTg                    >  1:507931/1‑81 (MQ=255)
                                                          cGTCAAAATAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTg                    >  2:507898/1‑81 (MQ=255)
                                                               aaaTAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGcga               <  2:507832/81‑1 (MQ=255)
                                                               aaaTAATGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGcga               >  1:507889/1‑81 (MQ=255)
                                                                aataatGGCGTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGcgac              >  2:507893/1‑81 (MQ=255)
                                                                  taatGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGcgacga            >  2:507878/1‑81 (MQ=255)
                                                                         gTAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGcc     <  2:507833/81‑1 (MQ=255)
                                                                          tAACCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGcc     >  1:507886/1‑80 (MQ=255)
                                                                            aTCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAAc  >  2:507941/1‑81 (MQ=255)
                                                                             |                                                                               
AACAAAAATGTGTCTGGCGCGGGTTTATGGTGTTTGACGTGATCGGCAGCGACGACGGCGTCAAAATAATGGCGTAATCCCAGGTGCGCCAGCAATGCCTCAGCGATGGCGCTTTCACTCCCCGTTCCTACAGCCATTGGGCGACGACCATGCCAAC  >  NC_000913/2817055‑2817211

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