Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,678,0390TG42.9% 86.3 / 16.5 56N108T (AAC→ACC) yphEputative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (7/25);  new base G (24/0);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 8.52e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGC  >  NC_000913/2677890‑2678186
                                                                                                                                                      |                                                                                                                                                   
ccaGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGGGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGt                                                                                                                                                     <  1:72868/150‑1 (MQ=255)
  aGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCTGACC‑GAGGGAGAGGgttt                                                                                                                                                  >  1:1001055/1‑151 (MQ=255)
      gctTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCg                                                                                                                                               <  2:1117451/150‑1 (MQ=255)
                   gCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCt                                                                                                                                 <  1:808016/151‑1 (MQ=255)
                   gCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGGTTTTCGCCCCTGTCAGccc                                                                                                                                  >  1:714149/1‑150 (MQ=255)
                        gggTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCtt                                                                                                                                <  2:672853/145‑1 (MQ=255)
                        gCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCCGACC‑GAGGCGGAGGGTTTTCGCCCCCGTCAGCCCCTcc                                                                                                                              >  2:879046/1‑149 (MQ=255)
                          gTTGAAACAAGTTGTTCAGGACTAACGTCACCGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTccacca                                                                                                                          <  2:43379/151‑1 (MQ=255)
                          gTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCAGTCCGGAGGC‑GAGGGTTTCTCGCACTGTCAGCCCCTccacc                                                                                                                           >  2:275313/1‑150 (MQ=255)
                             gAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCCACC‑GAGGCAGAGGGTTTCCGCCCCCGTCAACCCTTccacc                                                                                                                           >  2:1047492/1‑147 (MQ=255)
                              aaaCAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACC‑CAGGCGGAGGGTTTTCCCCACTGCCCGCCCTTccccc                                                                                                                           >  2:769190/1‑146 (MQ=255)
                               aaCAAGTTGTTCAGGACTAACGTCAACGCCCAGCCCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACtt                                                                                                                     <  1:187634/151‑1 (MQ=255)
                                   aGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACC‑GAGGGGGAGGGTTTTCGCCCCCTTCAGCCCTTCCCCCAGAcctta                                                                                                                   >  2:524650/1‑145 (MQ=255)
                                       gttCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTg                                                                                                             <  1:399925/151‑1 (MQ=255)
                                          cAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGAt                                                                                                           <  2:1032634/150‑1 (MQ=255)
                                           aGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑CAGGGAGAGGGTTTCCCCCCCTGTCAGCCCTTCCACCAGACCTAAATCCCGATa                                                                                                          >  1:9719/1‑150 (MQ=255)
                                              aCTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACC‑GCGGCAGAGGGTTTTCGCCCCTGTCAGCCCCTCCCCCAGACCTAACACCTGTTaaaac                                                                                                      >  2:418016/1‑149 (MQ=255)
                                                 aacgccaacgccCAGCTCCGGTAAGCAACGTGGGGCATCCGGGGCCATTGGCAGGTAATCAATCATGCCGTTGCGGCGGGCCCACCGACC‑GAGGCAGAGTTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACcgc                                                                                                    <  2:940047/150‑1 (MQ=255)
                                                 aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGTAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGGTTTTCGCCACCGTCAGCCCCTTCCCCACACATAACTCCTGATAAACcgc                                                                                                    >  1:323991/1‑150 (MQ=255)
                                                 aaCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCTCTTCCACCAGACTTAACTCCTGATAAACcgc                                                                                                    <  1:114452/150‑1 (MQ=255)
                                                  aCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACTGACC‑CAGGCAGAGGGTTTTCGCCCCTGTCCGCCCTTCCCCCACACTTAACTCCTTATAAAAcccgc                                                                                                  >  1:800046/1‑151 (MQ=255)
                                                  aCGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACC‑CAGGCGGAGGGTTTTCGCCCCCGGCCGCCCCTCCCCCAGACCTAACACCTGATaaa                                                                                                        >  2:464768/1‑145 (MQ=255)
                                                   cGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACcgcgcg                                                                                                 <  1:224684/151‑1 (MQ=255)
                                                   cGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACcgcgcg                                                                                                 <  2:259521/151‑1 (MQ=255)
                                                     tCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACC‑GAGGCAGAGGGTTTTCGCCCCCGTCAGCCCCTCCCCCAGACTTTACTCCTGATAAACcgcgc                                                                                                  >  2:1116955/1‑148 (MQ=255)
                                                      cAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCAc                                                                                              <  1:863084/151‑1 (MQ=255)
                                                       aaCGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCCCTGACC‑GAGGCAGAGGTTTTTCGCCACTGTCAGCCCCTCCCCCCCACTTAACTCCTGATAAAACGCGCGcccc                                                                                             >  2:1146181/1‑151 (MQ=255)
                                                             cAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGt                                                                                       <  2:1054037/151‑1 (MQ=255)
                                                                cgcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACC‑GCGGCCGAGGGTTTTCGCCACTGTCAGCCCCTCCCCCCGACTTAAATCCTGATAAAACGCGCGCACCCCCCGtt                                                                                      >  2:77794/1‑149 (MQ=255)
                                                                 gcCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACC‑GAGGCAGAGGTTTTTCGCCCCTGTCAGCCCCTCCACCCGACTTAACTCCTCTTAAACCGCGCGCACCCCCacatt                                                                                     >  1:216782/1‑145 (MQ=255)
                                                                    tGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGAAC‑GAGGGAGAGGGTTTCCCCCCCTGCCAGCCCCTCCACCCCACTTAACTCCTGATAAAACACGCCCGCCCCCCGATCAGCgg                                                                                >  1:62576/1‑151 (MQ=255)
                                                                       aaGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACC‑GAGGGAGAGGGTTTTCCCCCCTGTCAGCCCCTCCACCAGACTTAAATTCTGATAAAACCCGCGCACCCCCCCGTCAGCGGCa                                                                              >  1:50956/1‑150 (MQ=255)
                                                                           aaCGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCg                                                                         <  1:865100/151‑1 (MQ=255)
                                                                               ttGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGGAGAGGGTTTTCGCCCCTGTCAGCCCCTCCCCCCGACTTAAATCCTGATAAAACGCGCGCACCCCCAGGTCAGCGGCACGGCg                                                                         >  1:680881/1‑147 (MQ=255)
                                                                                tGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCCACC‑CAGGGAGAGGGGTTCCTCCACCGTCACCCCCTCCACCCGAATTAACTCCCGTTAAAACCCCCGCCCCCCCCGTTCACCGGCACGCCggggcc                                                                    >  2:592041/1‑146 (MQ=255)
                                                                                   gCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGGCACCGACC‑GAGGCAGAGGGTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGCTAAAACCCGCGCCCCCCCCGTTCAGCGGCACGGCGa                                                                        >  1:615485/1‑144 (MQ=255)
                                                                                     aTCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACCGACC‑CGGGCGGAGGGTTTTCGCCACTTTCCGCCCCTTCCCCAGCCTTAACTCCTGATAAAACGCGCGCCCCCCCACGTCCGcggcgcggcggggca                                                                    >  2:1118508/1‑141 (MQ=255)
                                                                                      tCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACC‑CAGGCAGAGGGTTTTCGCCACTGTCAGCCCTTCCACCCGACTTAACTCCTGATAAAACCCGCGCACCCCCCGTTCAGCGGCACGGCGGGTCAg                                                                   >  2:561485/1‑146 (MQ=255)
                                                                                          gcgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTCAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTc                                                          <  2:1008469/151‑1 (MQ=255)
                                                                                           cgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTc                                                          <  2:230977/150‑1 (MQ=255)
                                                                                           cgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCg                                                         <  1:690127/151‑1 (MQ=255)
                                                                                           cgcCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCg                                                         <  2:494063/151‑1 (MQ=255)
                                                                                              cATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCAACC‑CAGGGCGAGGGTTTTCGCCCCTGTCCGCCCTTCCCCCCGAATTAACTCCTTATAAAACCCGCGCACCCCCCCCTCAGTGGCCCggggggcacacgga                                                               >  2:773199/1‑130 (MQ=255)
                                                                                              cATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGGTTTTCGCCACTGTCAGCCCTTCCCCCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGCAGCTTCGTc                                                       >  1:410686/1‑150 (MQ=255)
                                                                                                 ttGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCGCTTACC‑GGGGGAGAGGGTTTTCGCCACCGTCAGGCCTTCCCCCAAACTTAAATCCTGCTAAACCCCGCGCACCCCCCCTTCCGCGGCACGCCGCGTCAGCGGAGCTTCGTcccc                                                    >  2:707659/1‑150 (MQ=255)
                                                                                                  tGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCt                                                   <  1:437490/150‑1 (MQ=255)
                                                                                                   gCAGGTAATCAATCATGCCGTTGCGGCGGGGCCCCCGACC‑GAGGCAGAGGGTTTTCCCCCCTGTCAGCCCTTCCCCCCGACTTAACTCCTGTTAAAACGCGCGCACCCCCCCTTTAGCGGCACGGCGAGGCAGCCGAGCTTCGTCACCtt                                                  >  1:95984/1‑150 (MQ=255)
                                                                                                         aatcaatcATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTc                                            <  2:710758/150‑1 (MQ=255)
                                                                                                                          cggcggGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCCTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATcc                             >  1:892613/1‑148 (MQ=255)
                                                                                                                           ggcggGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGa                          <  2:23978/150‑1 (MQ=255)
                                                                                                                            gcggGGCCACTGCCC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGat                         <  2:189384/150‑1 (MQ=255)
                                                                                                                              ggggCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATc                      >  1:874150/1‑151 (MQ=255)
                                                                                                                               gggCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATc                      >  1:5553/1‑150 (MQ=255)
                                                                                                                               gggCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCa                     <  1:1147691/151‑1 (MQ=255)
                                                                                                                                       gACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTa               <  2:723896/149‑1 (MQ=255)
                                                                                                                                       gACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATc             >  1:454274/1‑151 (MQ=255)
                                                                                                                                          c‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCgg          >  1:248448/1‑151 (MQ=255)
                                                                                                                                                   agGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGc  <  1:587487/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CCAGCTGCTTTTGCGCCGGGCTTAGCGTTGAAACAAGTTGTTCAGGACTAACGTCAACGCCCAGCGCCTGTAAGCAACGTTGGGCATCCTGCGCCATTTGCAGGTAATCAATCATGCCGTTGCGGCGGGGCCACTGACC‑GAGGCAGAGGTTTTCCGCCACTGTCAGCCCTTCCACCAGACTTAACTCCTGATAAACCGCGCGCACCCCCAGTTCAGCGGCACGGCGAGTCAGCGTAGCTTCGTCACCTTCCAGTCGCGTCTCGCCAATCCAGATATCACCGCTATCCGGGCGTTCGC  >  NC_000913/2677890‑2678186

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