Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,135,7810GT24.5% 119.1 / 13.1 50R197R (CGC→CGApitBphosphate transporter
Reads supporting (aligned to +/- strand):  ref base G (6/31);  new base T (12/0);  total (19/31)
Fisher's exact test for biased strand distribution p-value = 2.01e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTG  >  NC_000913/3135548‑3136008
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                    
tGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTGCGGGTAATTTCATAGCCGGACGCATTAATATTGACGACGAAGCCAGCAGGGGCAATCCCCCCCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGacc                                                                                                                                                                                                                                 <  1:387146/239‑3 (MQ=255)
       aGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGCAACAATCACCGCAATACGCGTCCAGAACGCCGGTGTAGCTGTGCCCTTTTTTCTTTTTCTGATCTTcggt                                                                                                                                                                                                                           >  2:39921/1‑236 (MQ=37)
       aGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGACTGGCGGGTTACGGTTGCCCTCTTTCTTTTTTTCTTCTTCCggc                                                                                                                                                                                                                          >  2:740173/1‑238 (MQ=37)
                        aCGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGt                                                                                                                                                                                                          <  1:39921/238‑1 (MQ=255)
                              tCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGGTTGCCTTTTTTTTTTTTGCGATCTTTCGGAATGTGGTGAATAAGGTCAc                                                                                                                                                                                                       >  1:691992/1‑235 (MQ=255)
                               cgcgGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGTCTTTTTTTTTTTTTTGTTCTCTCGGACAGCGGTGAATAAGGTCACCCttt                                                                                                                                                                                                  >  2:99489/1‑239 (MQ=37)
                                                    aGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTTCGTTTTCCTTTTTTTTTTTTTCGATCTTCCTGAAAGCGGTGGATATGAGCACGCTTTTTCTTCTCGCCCCCGTAGcc                                                                                                                                                                              >  2:780053/1‑237 (MQ=37)
                                                    aGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGt                                                                                                                                                                             <  2:255401/239‑1 (MQ=255)
                                                       cGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCgc                                                                                                                                                                          <  1:99489/239‑1 (MQ=255)
                                                            gCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGcagcag                                                                                                                                                                     <  1:12781/239‑1 (MQ=255)
                                                            gCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGcagcag                                                                                                                                                                     <  2:530651/239‑1 (MQ=255)
                                                                            cgaAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTTCCTTTTTTTTTTTTTTTCTATTTCGGAATGCGGGGAATTCTGTTACCCTTTTTCCCCCCTCTCCACTAGCGTCGCCACAGGGATATCATGCCTcc                                                                                                                                                     >  2:306107/1‑238 (MQ=25)
                                                                             gaAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTTCg                                                                                                                                                                                                                                    >  1:560007/1‑159 (MQ=255)
                                                                             gaAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCg                                                                                                                                                                                                                                    <  2:560007/159‑1 (MQ=255)
                                                                                      agGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGTTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGAc                                                                                                                                            <  1:798913/237‑1 (MQ=255)
                                                                                      agGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGAc                                                                                                                                            <  1:740173/237‑1 (MQ=255)
                                                                                         gccAATCCCCGCCCGCACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATGAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGAcc                                                                                                                                           <  2:390074/233‑1 (MQ=255)
                                                                                           cAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGcc                                                                                                                                      <  1:780053/238‑1 (MQ=255)
                                                                                                cccACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTcc                                                                                                                                                                                           <  1:775942/181‑1 (MQ=255)
                                                                                                   accaGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCg                                                                                                                                                                                                                                    >  1:332515/1‑137 (MQ=255)
                                                                                                   accaGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCg                                                                                                                                                                                                                                    <  2:332515/137‑1 (MQ=255)
                                                                                                    ccaGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGc                                                                                                                                                                                                                      <  1:287858/150‑1 (MQ=255)
                                                                                                    ccaGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGcc                                                                                                                                      >  1:383254/1‑229 (MQ=255)
                                                                                                    ccaGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGAAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGcc                                                                                                                                      <  2:383254/229‑1 (MQ=255)
                                                                                                    ccaGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTCCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTc                                                                                                                                                                                                                              <  1:86602/142‑1 (MQ=255)
                                                                                                        tACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTTCCTTTTTTCTTTTTTCGCTCTTCCGGAATGCGGTGAATACGGGCACGCTTTTTCGTCCCCCCCCAGTAGCGTCGCAGCAGGAATATAAGGCGTCCCGCAATGACACGGACGACGCTAGGGGcaa                                                                                                                        >  2:443649/1‑239 (MQ=37)
                                                                                                        tACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGc                                                                                                                                                                                                                      <  1:556242/146‑1 (MQ=255)
                                                                                                         aCCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCGCGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGa                                                                                                                          <  1:306107/236‑1 (MQ=255)
                                                                                                          ccAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATCCGCGTCCAGAATGGCGGATTACGGTTGCCTTTTTTTTTTTTTTGATCCTCCGGAATGCGGTGAATAAGGTCACGCTTTTTTTTCCCGCCTCAGtt                                                                                                                                                                                 >  2:775942/1‑180 (MQ=37)
                                                                                                              cATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGCTTACGTTTGCCCTTTTTTTTTTTTCGGTCTTCCGGAAGGGGGTGAATACg                                                                                                                                                                                                            >  2:287858/1‑150 (MQ=25)
                                                                                                                 tACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTGATCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGa                                                                                                                          <  1:862052/228‑1 (MQ=255)
                                                                                                                        ccGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACCTTTTCCTTTTTTTTTTTTTCGATCTTCCGGGATGCGGTGAATACGGCTACGCTTTTTTGTCCCCCCCCCGGACCGTCCGCGCAGGAATATCAGGCCTCCCGCAATGACCCAGCCCACGGTCCGGGGAACCATCAGCGCGgggac                                                                                                        >  2:380642/1‑238 (MQ=25)
                                                                                                                           aTCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTTTTTTTTCGATCTTTCGGGATGCGGTGAATACGGGCACGCTTTTTTGTCCCCCCCCCGTAGCGTCGCAGCAGGAATATCAAGCCCCCCCCCAAGACCCGGGCGCCGAACGGGGGaacaaacaac                                                                                                                >  2:862052/1‑226 (MQ=25)
                                                                                                                                tttttGTCCGTCGTTCGCGCCGTGCGAAAACGCCCCGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGAt                                                                                                                                <  1:279757/207‑1 (MQ=255)
                                                                                                                                        cgtcgtTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCa                                                                                        <  1:805861/239‑1 (MQ=255)
                                                                                                                                          tcgtTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTTTTTTTTCGATTTTCCGGAAAGCGGGGAAAAAGGGCCCGCTTTTTTTTCCCCCCCCCGTAGGGGCGCAGCGGGAAATATCGGGCCCCCCCCAAAACCCGGGCCACGATAGGGGGAACa                                                                                                                      >  2:279757/1‑207 (MQ=19)
                                                                                                                                           cgtTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGc                                                                                     <  1:443649/239‑1 (MQ=255)
                                                                                                                                                                                     tCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGcgcg                                           <  1:67534/239‑1 (MQ=255)
                                                                                                                                                                                         cgcAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGc                                                                                                                                                                                        <  1:451848/95‑1 (MQ=255)
                                                                                                                                                                                         cgcAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGc                                                                                     <  1:71370/193‑1 (MQ=255)
                                                                                                                                                                                               aCGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAAc                                 <  1:273100/239‑1 (MQ=255)
                                                                                                                                                                                                 gggTCCAGAATGGGGGTTTACGGTTGCCGTTTTTCTTTTTGCGATCTTCCGGAATGCGGTTAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCTCGCAGGTTTAACGCAGCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCg                               <  2:184752/237‑1 (MQ=255)
                                                                                                                                                                                                   gTCCAGAATGGCGGTTTACGTTTTCCTTTTTTTTTTTTTCGATCTTCCCGGATGCGGTGAATACGGTCACCCTTTTTCTTCCCGCTCCCGTAGCg                                                                                                                                                                              >  2:451848/1‑95 (MQ=255)
                                                                                                                                                                                                   gTCCAGAATGGCGGTTTACGTTTGCCTTTTTTTTTTTTGCGATCCTCCCGAAAGGGGTGAATACGGTCACGCTTTTTTGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCCCCATGACCCGGCCGACGATAGGGGAAACAAAAAGGGAGGGGAAAATTTTTGTCACCTCCCGCAGGGTTaaag                                                                           >  2:71370/1‑191 (MQ=255)
                                                                                                                                                                                                    tCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCGCGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATg                                                            <  1:782169/207‑1 (MQ=255)
                                                                                                                                                                                                                            ccTTTTTTCTTTTTGCGCTCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCgcagca                                                                                                                                                                      >  1:687063/1‑79 (MQ=255)
                                                                                                                                                                                                                            ccTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCgcagca                                                                                                                                                                      <  2:687063/79‑1 (MQ=255)
                                                                                                                                                                                                                               tttttCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGtgtg  <  1:12797/239‑1 (MQ=255)
                                                                                                                                                                                                                                 tttCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTgg                                                                                                >  2:63730/1‑143 (MQ=255)
                                                                                                                                                                                                                                 tttCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTgg                                                                                                <  1:63730/143‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                    
TGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGTTTGCCTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATACGGTCACGCTTTTTCGTCCCGCTCCAGTAGCGTCGCAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAACAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTG  >  NC_000913/3135548‑3136008

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