Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,744,3620GA79.0% 74.1 / 15.5 38H104Y (CAC→TAC) aptadenine phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (5/3);  new base A (17/13);  total (22/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACG  >  CP000730/1744228‑1744500
                                                                                                                                      |                                                                                                                                          
aattATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGt                                                                                                                                      >  1:56035/1‑141 (MQ=255)
  ttATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTTGTTTACTTGCATCTTTGTACATTGTTa                                                                                                                                    >  2:185185/1‑141 (MQ=255)
  ttATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTa                                                                                                                                    <  1:317225/141‑1 (MQ=255)
    ataAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTaaa                                                                                                                                  >  1:9936/1‑141 (MQ=255)
               aCCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTAcc                                                                                                                       <  1:63990/141‑1 (MQ=255)
                ccTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCa                                                                                                                      >  2:98226/1‑141 (MQ=255)
                ccTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCa                                                                                                                      >  2:87636/1‑141 (MQ=255)
                  tactacGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCata                                                                                                                    <  2:311872/141‑1 (MQ=255)
                     tacGATACCGCCTAATTTTTCAACTAATTGTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTCACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTc                                                                                                                 <  2:388997/141‑1 (MQ=255)
                      acGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCt                                                                                                                >  2:117875/1‑141 (MQ=255)
                      acGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCt                                                                                                                <  2:291039/141‑1 (MQ=255)
                      acGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCt                                                                                                                <  1:318639/141‑1 (MQ=255)
                            ccGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCa                                                                                                          >  2:141389/1‑141 (MQ=255)
                                        tCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGa                                                                                              <  1:274215/141‑1 (MQ=255)
                                            cTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGa                                                                                          <  1:26075/141‑1 (MQ=255)
                                              aaTTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACt                                                                                        >  1:152237/1‑141 (MQ=255)
                                                   tATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACg                                                                                   >  1:159585/1‑141 (MQ=255)
                                                     ttGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAg                                                                                 <  1:5777/141‑1 (MQ=255)
                                                     ttGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAg                                                                                 >  1:71577/1‑141 (MQ=255)
                                                       gctgctTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGt                                                                               <  1:352957/141‑1 (MQ=255)
                                                       gctgctTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGt                                                                               >  2:329279/1‑141 (MQ=255)
                                                          gctTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACCCACGTTGACCTGGTTTAATTTCATCTTTGTACATTGTTCAAACATTTGTACCATATTCTAGGTCATACTCATAACAAATGACTTCACGAGGTAAt                                                                            >  2:146632/1‑141 (MQ=255)
                                                           ctTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAAtt                                                                           <  2:52717/141‑1 (MQ=255)
                                                                         ccaGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTa                                                             >  2:407084/1‑141 (MQ=255)
                                                                           aGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAAc                                                           <  1:117875/141‑1 (MQ=255)
                                                                              aGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAgg                                                        >  1:359994/1‑141 (MQ=255)
                                                                                cTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTg                                                      <  2:401701/141‑1 (MQ=255)
                                                                                 taataaATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGc                                                     <  1:38082/141‑1 (MQ=255)
                                                                                  aataaATCATCTGTAATTAACACACGTTGGCCTGGTTTAATTGCATCTTTGTACATTGTTAACAGATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCa                                                    >  1:153968/1‑141 (MQ=255)
                                                                                     aaaTCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTACATTGGTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAg                                                 >  2:355923/1‑141 (MQ=255)
                                                                                                      cacacGTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGTAAAGCCAATCCCCATTGAATa                                >  1:279054/1‑141 (MQ=255)
                                                                                                          cGTTGACCTGGTTTAATTGCATCTTTTTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCAATTGAATAAGCt                            >  1:200693/1‑141 (MQ=255)
                                                                                                           gTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTa                           <  2:217818/141‑1 (MQ=255)
                                                                                                           gTTGACCTGGTTTAATTGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTa                           >  1:398818/1‑141 (MQ=255)
                                                                                                                        aaTTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCTCGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAAt              >  2:162860/1‑141 (MQ=255)
                                                                                                                          ttGCATCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGa            >  1:164815/1‑141 (MQ=255)
                                                                                                                              aTCTTTGTACATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCAGAACGAATGACGTCACGAGGTAATTTCCCTTCGTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATaaa        <  1:185185/141‑1 (MQ=255)
                                                                                                                               tCTTTGTACATTGTTACAACATTTGTACAATAGTCTAGGTCATACTCATGATGAATGACTCCAGGAGGTAATTTCCTTTCTTTTCTGACAGGCGCAAAGCCAATCCCCATTGAATAACCTACAGAACTGCCAATGATaacc       >  2:29995/1‑139 (MQ=255)
                                                                                                                                    gtgCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAcg  <  1:313107/141‑1 (MQ=255)
                                                                                                                                    gtgCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGTCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCAcg  <  1:53592/141‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                          
AATTATAAATGCAATACCTACTACGATACCGCCTAATTTTTCAACTAATTTTATTGCTGCTTCAATCGTACCACCAGTAGCTAATAAATCATCTGTAATTAACACACGTTGACCTGGTTTAATTGCATCTTTGTGCATTGTTAAAACATTTGTACCATATTCTAGGTCATACTCATAACGAATGACTTCACGAGGTAATTTCCCTTCTTTTCTAACAGGTGCAAAGCCAATCCCCATTGAATAAGCTACAGGACAGCCAATGATAAAGCCACG  >  CP000730/1744228‑1744500

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