Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,782,9080TG33.9% 78.5 / 14.0 67T252P (ACC→CCC) gpsAglycerol‑3‑phosphate dehydrogenase (NAD+)
Reads supporting (aligned to +/- strand):  ref base T (8/33);  new base G (21/0);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 2.33e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCACGACCTAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATA‑‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCA  >  NC_000913/3782670‑3783145
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
gCACGACCTAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGt                                                                                                                                                                                                                                                     <  1:1922440/239‑1 (MQ=255)
                tCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTTCACGCCTGGTTGTTGGTACAGGTAAAGACCaaaatcccc                                                                                                                                                                                                                                     >  2:1618458/1‑231 (MQ=255)
                          cTCGCGCGTGTTTTTTCCCAATTGTATATTTTGAAAACTTCCCTCGTTTATTGCTATTCAACCGCCCAAGCGGGCCGCATTTGCGCGGTTTTCTTTAGTATTGCGTTGCCCTCCCCCCCCCTCCCCACTTTCCCCTGGCGCGCTTTGTCAATCCAGCCCCTGACCGCGCATCATCCAAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGATAAGCAACAGATCGCAAAGCCCCGCCa                                                                                                                                                                                                                          <  2:1963321/239‑1 (MQ=255)
                           tCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCAt                                                                                                                                                                                                                         <  1:2311298/239‑1 (MQ=255)
                           ccgcgcgctTTTTTTCCGCAATATAATACTTGATAAATTCCCTCGTTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGCCTTCTTTCGTATTGCGGTAGCCTTCCCCCACCGGACCAATCTTCTCCTGCGCGCTTTGTCCATCCATGCCCGGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGTTTGTGGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCAt                                                                                                                                                                                                                         <  2:42172/238‑1 (MQ=255)
                              cgcgcgTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGGTACGCGACTGGGTGTCGGGACAGGGAAACACCACAACGCCAAGCCCCGCCACGc                                                                                                                                                                                                                       >  2:2058266/1‑238 (MQ=255)
                                            tcaaattaataCTTCATAAGTTTCCTCAGTTATTGCCATTCCAACGCCTGAGGGCTCGGTCATTTCGCGATCTCCTTCCGAATGTGGGAAGCCATCCTCCTCCTGACCATTCTACACCTGGACGCTTTGTCCACCCGTGCCTTGACCAAGAGCCACGCCAACAAGGCGGTTACGCGCTTGGTTGCCGGTGCGGGAAACCGCCACATCACCACGCCCCGCCATGCCCATAAAGGTGGCAg                                                                                                                                                                                                       <  2:65693/233‑1 (MQ=255)
                                                ataaagaCTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGGGATCCCCCATTTCGGGGACTTCTTTCGTATTGCGGTAGCCTTCCCCCCCCTCACCAATCTTCCCCTGCCCGCTTTGTACATCCATGCCCTGACCGAGCATCATCCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGcc                                                                                                                                                                                                                           <  2:34760/210‑1 (MQ=255)
                                                           gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACACGACTGGTTGTCGGTACAGGGAAGCACCAGATCGCCAAACCCCGCCATGCCCCTAAAAGGGGCAGGGTCGGGACCCCgc                                                                                                                                                                                         >  1:371647/1‑239 (MQ=255)
                                                           gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTAAGCGACTGGGTGTCGGGGCAGGGAAGCACCCGAAAGCCCAACCCCGCCCAGcccccaa                                                                                                                                                                                                                 >  1:34760/1‑212 (MQ=255)
                                                           gATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGGTACACGACCGGGTGTTGGGACAAGGAAGCACCAAAACGCCAAGCCCCGCCCCGCCCCTCAAAGGGGGCAGGGGGGGCCCCCac                                                                                                                                                                                          >  2:82930/1‑238 (MQ=255)
                                                                           ttgtGGAATTCCAACGCCGAGGAGTTGCCCCAGTTGGCGGACTTCTTGCGTAATGGGGTGTCCTCCCCCCACCTGCCCAATCTTCCCCTGCGCGCTTTGAACACCACTGCCCGGACCGAGCATCATGCCAAAACGGCGGTTACGCGTCTGGTTGTCGGTACAGGTAAGCGCCAGACCCCCAACCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGAcgcg                                                                                                                                                                         <  2:1143899/236‑1 (MQ=255)
                                                                                     tCACCGCCGAAGCGAGGCGCCAGTTCGCGTACTTTTTTCGTATTGCGGTAGCCTTCCCCCACCTGCCCAATCTTCCCCTGCGCGCTTTGTACATCCATCCCCTGACCGAGCATCATGCCAAACCGGCGGTTACGCGACTGGTTGTCGTTACAGGTAAGCACCAGATCGCAAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc                                                                                                                                                              <  2:1639335/239‑1 (MQ=255)
                                                                                     tCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCCCCACCTGACCAATCTTCTCCTGCGCGTTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc                                                                                                                                                              <  1:1618458/239‑1 (MQ=255)
                                                                                     ccaccccggaggCGGTGCGCCGTTCCGCAGATTTCTTTCGTATTGCGGTACCCTTCCCCCCCCGACCCATTCTTTCCTAGCGGGCTTTAAACATCAATGCCCTGCCGGAGCACCATGCAAAACCGGGGGTTACGGGATGGGTTGTCGGTAGAGGTAGGCCCCAGCCCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc                                                                                                                                                              <  2:397921/228‑1 (MQ=255)
                                                                                           ccGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGGACAGGGAAGCCCCAAAACGCCCAGCCCCGCCCTGCCCACAAAAGAGGCAGGGGGGGGACCCCGCCCCGCACCCAAACACGCACATTCAGCCAGccc                                                                                                                                                        >  2:920177/1‑239 (MQ=255)
                                                                                                                      tcttGCGTATTGCGGTAGCCTCCCACCCCCTGACCATTCTTCCCCTGCGCGCTTGGTACCTCCATGCCCTGACCGAGCATCATGCCAAACCGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAACCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCAtt                                                                                                                             <  2:654680/239‑1 (MQ=255)
                                                                                                                                  cGGTAGCCTTCCCCCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACGGGTAAGCACCAGATCGCCAAGCCCCGcc                                                                                                                                                                                                                           <  1:1725777/135‑1 (MQ=255)
                                                                                                                                        ccTTCCCCCCCCTGCCCAATCTTCCCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGAt                                                                                                                                               <  1:20264/203‑1 (MQ=255)
                                                                                                                                        ccTTCCCCCACCTGCCAAATCTTCCCCGGCGCGTTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCAcc                                                                                                                                                                                             <  2:1141504/158‑1 (MQ=255)
                                                                                                                                        ccTTCCACCCCCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc                                                                                                                                                              <  1:70659/188‑1 (MQ=255)
                                                                                                                                         cTTCCCCCCCCTGACCACTCTTCTCCTGCGCGCTTTGTACATCCATGCCCGGCCCGAGCATCATGCCAAAACGGCGGTTACGCGACTGTTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCg                                                                                                         <  1:1558310/239‑1 (MQ=255)
                                                                                                                                         cTTCCCCCACCTGACCAATCTTCCCCTGCGCGTTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGCCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCg                                                                                                         <  1:2058266/239‑1 (MQ=255)
                                                                                                                                            cccccaccTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTCCGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAg                                                                                                                                                                                                                                  <  1:1208155/115‑1 (MQ=255)
                                                                                                                                              ccccccTGACAAATTTCCTCCTGCCGCTTTTGCACCTCCGTCCCCGGCCCGGGCATCATCCCAAAAGGGGGGTTCCGCGACGGGTTGTCGGTACAGGAAAGCCCCAGACCCCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGa                                                                                                    <  2:2039411/235‑1 (MQ=255)
                                                                                                                                               cccccTGCCCAATCTTCCCCTGCGCGCTTTGTACTTCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGTTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGAc                                                                                                   <  1:1251528/236‑1 (MQ=255)
                                                                                                                                                  ccTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGGACAGGGAAGCACCAAAACGCCAAACCCCGCCCCGCCCCTAAAAGGGGGAGGGGCGGCACCCCGCGCCGCACCAAAACGCGACACTTTTgcccgccccccg                                                                                                                                                    >  2:70659/1‑184 (MQ=255)
                                                                                                                                                  ccTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGGTGTCGGGACAGGGAAGCACCAGATCGCCAAACCCCCCCACGCCCCTAAAAGGGGGAGGGGCGGCCCCCCAcgccccc                                                                                                                                                                                   >  1:1141504/1‑155 (MQ=255)
                                                                                                                                                  ccTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCGACCGGGTTTCGGTGCAGGGAAACACCAAAAAGCCAAACCCCCCCCCGCCCCCAAAGGTGGCAGGGGCGGGACCCCCCGCCGCACCCAGACGCCGCAATTCACCCAGCCCACCGGGGGCGAGCGCCgac                                                                                                                                     >  2:20264/1‑201 (MQ=255)
                                                                                                                                                     gACCATCCTCCTCCTGCGGGTTTTGTACACCCATCCCCTGACCGAGCACCATGCCAAACCGGCGGTTCCGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATcccc                                                                                             <  2:1055794/239‑1 (MQ=255)
                                                                                                                                                      aCCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCGACTGGTTGTTGGGGCAGGGAAGCACCCGATCGCCAAGCCCCGCCCTg                                                                                                                                                                                                                        >  2:1208155/1‑118 (MQ=255)
                                                                                                                                                      aCCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCCACTGGGTTTCGGGACAGGGAAGCACCAAAACGCCCAGCCCCGCCCCGCCCATAAAAGGGGCAGGGGCGGCCCCCAGCGCCCCACCAACACGCCACACTTTAACCCACCCCCGGGGGATCCGCCCCGTACGCGCACTCGCCCCCAAAACCCGACCGTCGGACGAcccc                                                                                             >  2:293968/1‑238 (MQ=255)
                                                                                                                                                                  cTGCGCTCTTTGTACACCCATGCCCTGACCGAGCATCATGCAAAAACGGCGGTTACGCGACGGGTTGGCGGTAAAGGAAAGCCCCAGACCGCCAACCCCCGCCATGCCCATAAAGGTGGCAGGGt                                                                                                                                                                                                    <  2:1234756/125‑1 (MQ=255)
                                                                                                                                                                     gcgcgtTTTGTACACCCATGCCCTGCCCGAGCATCATGCCAAACCGGCGGTTCCGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATa                                                                              <  1:920177/239‑1 (MQ=255)
                                                                                                                                                                      cgcgcTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCCCCTGGTTTTCTGTACATGTTAGCACCCAAGAGCCCAGACCCGCCCTccccctcaa                                                                                                                                                                                                                >  1:65699/1‑106 (MQ=255)
                                                                                                                                                                            tttGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcag                                                                                                                                                               <  1:2310934/153‑1 (MQ=255)
                                                                                                                                                                              tGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAg                                                                                                                                                                                                              <  1:1564414/105‑1 (MQ=255)
                                                                                                                                                                              tGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagc                                                                                                                                                              <  1:1749905/152‑1 (MQ=255)
                                                                                                                                                                              tGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTAcgcg                                                                                                                                 <  1:2129244/181‑1 (MQ=255)
                                                                                                                                                                              tGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCCACTGGGTTTCGGGACAGGGAAGCACCAAAACGCCAAACCCCCCCCTGCCCATAAAAGGGGCAGGGGCGGCAccccg                                                                                                                                                                                          >  1:1234756/1‑123 (MQ=255)
                                                                                                                                                                              tGTACACCCAGCCCCTGCCCGAGCATCATGCCAAAACGGCGTTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCCCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGAc                                                                                                   <  1:1636052/210‑1 (MQ=255)
                                                                                                                                                                                tACATCCATGCCCTGACCGAGCATCATGCAAAAACGGCGGTTCCGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCAt                                                                                                                                                                                                                         <  1:1523891/93‑1 (MQ=255)
                                                                                                                                                                                     ccagccccTGACCGAGCATCATCCCAAACCGGCGTTTCCGCGACTGTTTGTCGGTACAGGAAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAg                                                                                                                                                                                                       <  2:225798/100‑1 (MQ=255)
                                                                                                                                                                                     ccATGCCCTGACCGAGCATCATGCCAAACCGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGcc                                                                                                                                                                                                                           <  1:1150712/86‑1 (MQ=255)
                                                                                                                                                                                     ccATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTCCGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCAcc                                                                                                                                                                                             <  1:71044/115‑1 (MQ=255)
                                                                                                                                                                                     ccATGCCCGGACCGAGCATCAGGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCAAGACCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTcagccag                                                                                                                                                           <  1:1331345/148‑1 (MQ=255)
                                                                                                                                                                                      cATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTAAGCGACCGGTTTTCGGGACAAGGAAACACCAAAACGCCCAACCCCCCCCCGCCCCTAAAAGGGGGAGGGGCGGCACCCACCCCCCCCCCCCAAAGCGGCCTTTTAGCCAGCCCCCCGGGGCTCAGCGCCGCACCCGCCATTGCCCCACAAACAAAACAGTCGGACATCCCCCCACCCAACGCAATAACACGTTTTTCCCcg                                                               >  2:2075368/1‑238 (MQ=255)
                                                                                                                                                                                        tGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGGTGTCGGGACAGGGAAGCACCAGAACCCCCAGCCCCCCCCCGCCCCTAAAAGTGGCCgggg                                                                                                                                                                                                    >  2:1564414/1‑104 (MQ=255)
                                                                                                                                                                                        tGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACCCCCCACGTTTTCGGTAAAGGGAAGGACCACAACGCCAAAACCCCCCCACTCCCCAAAAGGGGCAGGGGGGGGACCCCGCGACGCACCCAGAAGCCAACATTCACGCAGCCCACCGGGGGGAAGCGCCGGCCGCCCATTCGCCCCCAAAACCGAACCCGCGGGCAACCCCCCCCCCCTGGCAAAA‑‑ACGcttttttcccccgc                                                           >  2:1141797/1‑229 (MQ=255)
                                                                                                                                                                                        tGCCCTGACCGAGCATCATGCCAAAACGGCGGTTAACCCACTGGTTGTTGGGACAGGGAAGCACCAAAACGCCAAACCCCCCCCCGCCCCTAAAAGGGGGAGGGGCGGCACCCCCCGCCCCCCCCACACGCGACATTTTAGCCAGccccccg                                                                                                                                                    >  2:1749905/1‑148 (MQ=255)
                                                                                                                                                                                        tGCCCTGACCGAGCATCATGCCAAAACGGCGGGTACGCGACTGGGTGTGGGGTCAGGTAAACa                                                                                                                                                                                                                                                >  2:759142/1‑63 (MQ=255)
                                                                                                                                                                                          cccTGACCGGGCATCATGCAAAAAGGGAGGTTACGCGACGGTTTGTGGGTACGGGTAAGCACCAGATCGCCAGGACCCGCCATGCCCATAAAGGTGGCAGGGTCg                                                                                                                                                                                                  <  1:90584/105‑1 (MQ=255)
                                                                                                                                                                                          cccTGACCGAGCATCATGCCAAAACGGCGGTTAAGCGGCCGGGTGTGGGGAAAGGGAAg                                                                                                                                                                                                                                                  >  2:594762/1‑59 (MQ=255)
                                                                                                                                                                                          cccTGACCGAGCATCATGCCAAAACGGCGGTTAAGCGACCGGGTTTCGGGACAGGGAAGCAACAAAACGCCCAGCCCCCCCCTGCCCCTaaaa                                                                                                                                                                                                              >  2:1523891/1‑92 (MQ=255)
                                                                                                                                                                                           ccTGACCGAGCATCATGCCAAACCGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCCCCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAg                                                                                                                                                                                                       <  2:1511480/99‑1 (MQ=255)
                                                                                                                                                                                               aCCGAGCATCATGCCAAAACGGCGGTTTCGCGACTGGTTGTCGGTACAGGGAAGCACCCGACCGCCAAGACCCCCCCTGCCCCTAAAGGAGGGCGGGGCGGCAcc                                                                                                                                                                                             >  1:225798/1‑105 (MQ=255)
                                                                                                                                                                                               aCCGAGCATCATGCCAAAACGGCGGTTACGCGACTGTTTGTCGGTAAAGTTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACgg                                                                                                                                                    <  1:1414438/145‑1 (MQ=255)
                                                                                                                                                                                               aCCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCAccc                                                                                                                                                                                            <  1:607072/106‑1 (MQ=255)
                                                                                                                                                                                               aCCGAGCATCATGCCAAAACGGCGGTTACGCGACCGGGTTTCGGGACAGGGAAAGACCAAAACGCCAAACCCCCCCCTGCCCCTAAAAGGGGCAGGGGCGGCACCCCGcgccgcc                                                                                                                                                                                   >  2:71044/1‑114 (MQ=255)
                                                                                                                                                                                                     catcatGCCAAAACGGCGGTTACGCGGCTGGTTGTTGGGACAGGGAAGCACCCGAACGCCCAGCCCCCCCCCGCCCCTAAAAGGGGGAGGGGCGGCAcc                                                                                                                                                                                             >  1:1511480/1‑99 (MQ=255)
                                                                                                                                                                                                         atGCCAAAACGGCGGTTACGCGCCTGGTTGTCGGTACAGGGAAGCACCAAAACGCCCAACCCCCCCCTGCCCCTAAAAGGGGCCGGGGGGGCACCCCGCCCCCCAc                                                                                                                                                                                  >  2:607072/1‑106 (MQ=255)
                                                                                                                                                                                                                      ggTTACGGGACTGGTTGTCGGTACGGGAAAGCCCCAGTTCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATA‑‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCgg                                <  1:82930/235‑1 (MQ=255)
                                                                                                                                                                                                                                      gTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATA‑‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGaa              >  1:938852/1‑237 (MQ=255)
                                                                                                                                                                                                                                              aGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCAtt                                                                                                                             <  2:538476/121‑1 (MQ=255)
                                                                                                                                                                                                                                              aGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCAtt                                                                                                                             >  1:538476/1‑121 (MQ=255)
                                                                                                                                                                                                                                              aGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATA‑‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGcc     >  2:159869/1‑238 (MQ=255)
                                                                                                                                                                                                                                              aGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATA‑‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAAACCGGATTGCTGTAAACGCGGAAACTTTTGcc     >  1:2105771/1‑238 (MQ=255)
                                                                                                                                                                                                                                                gTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATA‑‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCa  >  1:504108/1‑239 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                   
GCACGACCTAGTAAAGTCAATGCTGCCTCGCGCGCGTTTTTTCCGCAATATAATACTTGATAAATTTCCTCGGTTATTGGCATTTCAACGCCGAAGCGATGCGCCAGTTCGCGGACTTCTTTCGTATTGCGGTAGCCTTCCACCACCTGACCAATCTTCTCCTGCGCGCTTTGTACATCCATGCCCTGACCGAGCATCATGCCAAAACGGCGGTTACGCGACTGGTTGTCGGTACAGGTAAGCACCAGATCGCCAAGCCCCGCCATGCCCATAAAGGTGGCAGGGTCGGCACCCAGCGCCGCACCAAGACGCGACATTTCAGCCAGCCCACGGGTGATCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATA‑‑ACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCA  >  NC_000913/3782670‑3783145

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