Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007779556,6530AC40.8% 111.3 / 30.9 71G104G (GGT→GGG) folDbifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase/5,10‑methylene‑tetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base A (34/8);  new base C (0/29);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 5.68e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTC  >  NC_007779/556553‑556751
                                                                                                    |                                                                                                  
gACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATa                                                                                                    <  3:15685/101‑1 (MQ=255)
   gcggcgCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATAccc                                                                                                 >  3:259448/1‑101 (MQ=255)
    cggcgCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACccg                                                                                                <  1:1460599/101‑1 (MQ=255)
     ggcgCCCGCTGGCACAGACGACCGACGATGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACccgc                                                                                               <  1:1262048/101‑1 (MQ=255)
      gcgcgcgcTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACccgcc                                                                                              >  3:941182/1‑101 (MQ=255)
       cgcgcgcTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACccgccg                                                                                             <  1:1562125/101‑1 (MQ=255)
        gcgcgcTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCgg                                                                                            <  3:1472200/101‑1 (MQ=255)
        gcgcgcTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGCCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCgg                                                                                            >  1:172321/1‑101 (MQ=255)
                    aGACGCCCGACTTTGTAAGGATGGAACCCGTCCACGTCTTTGTCCGGATGAATCCGTTCCAGCATTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTg                                                                                <  3:314271/101‑1 (MQ=255)
                       cgaccgacGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAAc                                                                             >  3:1551573/1‑101 (MQ=255)
                            gacGTTGTAGGGATGGAACCCGTCCACGTCTTTGTCCGGATGAATCCGTTCCAGCACTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGaa                                                                        <  1:407561/101‑1 (MQ=255)
                            gacGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGaa                                                                        >  3:871322/1‑101 (MQ=255)
                              cGTTGTAAGGATGGAACCCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCATTTGGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATg                                                                      <  1:1102308/101‑1 (MQ=255)
                              cGTTGTAAGGATGGAAACCGCCCACGTCTTTGTCCGGATGAATACGTTCCAGCATTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATg                                                                      <  1:97521/101‑1 (MQ=255)
                               gTTGTAAGGATGGAACCCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCATTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGc                                                                     <  3:1181955/101‑1 (MQ=255)
                                 tGTAAGGATGGAACCCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCa                                                                   <  3:679377/101‑1 (MQ=255)
                                     aGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCAGTAACGGCAGTTGAACCAGAATGCCATCGa                                                               >  1:243019/1‑101 (MQ=255)
                                       gATGGAACCCGTCCCCGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATg                                                             <  3:472622/101‑1 (MQ=255)
                                       aacggACACCGTCCCTGTCTTTTGCCGGTTAAATACATTCCAGCAGTTTGACGTTAAAAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATg                                                             <  1:1493103/98‑1 (MQ=255)
                                        aTGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATgg                                                            >  3:1097607/1‑101 (MQ=255)
                                             aaCCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGtgttg                                                       >  3:1545958/1‑101 (MQ=255)
                                               ccGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTc                                                     >  3:540767/1‑101 (MQ=255)
                                               ccGTCCACGTCTTTGTCAGGATAAATACGTTCCAGCACTTTGCTGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTc                                                     <  3:1170541/101‑1 (MQ=255)
                                                cGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCCCTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCg                                                    <  1:122304/101‑1 (MQ=255)
                                                 gTCCACGTTTTTGTCCGGATAAATACGTCCCAGCACTTTGATGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCgg                                                   <  3:32755/101‑1 (MQ=255)
                                                 gTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCgg                                                   >  3:1396164/1‑101 (MQ=255)
                                                  cccACGTCTTTGTCCGGATAAATACGTTCCAGCACTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGc                                                  <  1:1010216/100‑1 (MQ=255)
                                                   ccACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCa                                                 >  1:545212/1‑101 (MQ=255)
                                                   ccACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCa                                                 >  1:402449/1‑101 (MQ=255)
                                                    cACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCAt                                                >  2:2462/1‑101 (MQ=255)
                                                    cACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCAt                                                >  1:976699/1‑101 (MQ=255)
                                                      cGTTTTTGTCGGGATAAATACGTCCCAGCTTTTTGACTTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTc                                              <  3:983949/101‑1 (MQ=255)
                                                       gtttttGTCCGGAAGAATGCGTCCCAGCATTTTGATTTTACAAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCa                                             <  3:601810/98‑1 (MQ=255)
                                                       gtttttGCCCGGATAAATACGTTCCACCATTTTGACTTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCa                                             <  3:1088553/98‑1 (MQ=255)
                                                        tCTTTGTCCGGATGAATACTTTCCAGCACTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAg                                            <  3:447290/101‑1 (MQ=255)
                                                        tCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAg                                            >  1:879563/1‑101 (MQ=255)
                                                         cTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGc                                           >  3:1089943/1‑101 (MQ=255)
                                                          tttGTCGGGATAAGTACGTTCCAGCTGTTTACGGTTATAATCACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCg                                          <  3:1196402/101‑1 (MQ=255)
                                                          tttGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCg                                          >  1:116253/1‑101 (MQ=255)
                                                           ttGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGt                                         <  1:628382/101‑1 (MQ=255)
                                                             gTCGGGATAAATACGTTCCAGCATTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTAt                                       <  1:1366460/101‑1 (MQ=255)
                                                               ccGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCg                                     >  3:1078341/1‑101 (MQ=255)
                                                                cGGCTGACTCCGTTCCCGCACTTTGCCGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGa                                    >  4:64241/1‑101 (MQ=255)
                                                                cGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGa                                    >  1:1036896/1‑101 (MQ=255)
                                                                 ggATGAATCCGTTCCAGCATTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGAt                                   <  3:622798/101‑1 (MQ=255)
                                                                 ggATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGAt                                   >  1:115428/1‑101 (MQ=255)
                                                                  gATAAATACGTCCAAGCCCTTTGATTTTATGAACACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATa                                  <  3:909671/101‑1 (MQ=255)
                                                                   aTGAATCCGTTCCAGCACTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATaa                                 <  1:725212/101‑1 (MQ=255)
                                                                   aTGAATACGTTCCAGCACTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATaa                                 <  3:26825/101‑1 (MQ=255)
                                                                   aTGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATaa                                 >  1:1209945/1‑101 (MQ=255)
                                                                         aCGTTCCAGCTCTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCca                           <  3:260861/101‑1 (MQ=255)
                                                                           gtcccAGCACTTTGACTTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCcagc                         <  1:709869/98‑1 (MQ=255)
                                                                           gTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCcagc                         >  1:1040394/1‑101 (MQ=255)
                                                                           gTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCcagc                         >  1:589880/1‑101 (MQ=255)
                                                                            ttCCAGCACTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCcagca                        <  1:1404106/101‑1 (MQ=255)
                                                                            tcccAGCATTTTGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCcagca                        <  1:1097822/99‑1 (MQ=255)
                                                                             cccAGCACTTTGACTTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCcagcag                       <  1:55226/100‑1 (MQ=255)
                                                                                aGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTc                    >  2:110409/1‑101 (MQ=255)
                                                                                 gctttttGACGTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTcc                   <  1:974238/97‑1 (MQ=255)
                                                                                   aCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGc                 >  1:484035/1‑101 (MQ=255)
                                                                                   aCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGc                 >  3:313829/1‑101 (MQ=255)
                                                                                   aCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGc                 >  3:833054/1‑101 (MQ=255)
                                                                                      ttGAGGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTc              <  3:777999/101‑1 (MQ=255)
                                                                                      ttGACTTTATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTc              <  3:542696/101‑1 (MQ=255)
                                                                                       tGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCg             >  2:104232/1‑101 (MQ=255)
                                                                                         actttATCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCt           <  3:588756/98‑1 (MQ=255)
                                                                                          ctttACCAATCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCtg          <  3:685264/99‑1 (MQ=255)
                                                                                           gTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCCTCGCtgg         >  1:13817/1‑101 (MQ=255)
                                                                                               tCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTggtggt     >  3:268577/1‑101 (MQ=255)
                                                                                                cAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCtggtgggt    >  3:970142/1‑99 (MQ=255)
                                                                                                cAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGtt    >  1:1398857/1‑101 (MQ=255)
                                                                                                  aTACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTc  >  1:1006715/1‑101 (MQ=255)
                                                                                                    |                                                                                                  
GACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTC  >  NC_007779/556553‑556751

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑