Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,604,2370CT43.1% 64.1 / 24.7 59S64L (TCG→TTG) lsrDautoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand):  ref base C (10/23);  new base T (25/0);  total (36/23)
Fisher's exact test for biased strand distribution p-value = 1.13e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTT  >  NC_000913/1604088‑1604385
                                                                                                                                                      |                                                                                                                                                      
gCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTCg                                                                                                                                                       <  1:599453/151‑1 (MQ=255)
  tCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTTGtt                                                                                                                                                     >  2:82023/1‑151 (MQ=255)
   cGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTCGt                                                                                                                                                      <  1:579218/149‑1 (MQ=255)
     ttATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTTTTTTgg                                                                                                                                                  >  2:877576/1‑151 (MQ=255)
      tATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATA‑TTTTGttt                                                                                                                                                    >  1:425858/1‑148 (MQ=255)
          gagaTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTA‑TTTTGTTTGGtttgt                                                                                                                                             >  2:449091/1‑148 (MQ=255)
           tgaTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTTGTTTGGt                                                                                                                                                 >  2:309811/2‑146 (MQ=255)
           tgaTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTCGTTTGGt                                                                                                                                                 <  1:309774/145‑1 (MQ=255)
             aTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTCGTTTGGTTC‑GACCa                                                                                                                                         <  2:186043/151‑1 (MQ=255)
               tGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTTGTTTTGGTT‑GACCATc                                                                                                                                       >  2:479764/1‑151 (MQ=255)
                   gCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTCGTTTGGTTC‑GACCATCGGcc                                                                                                                                   <  2:425916/151‑1 (MQ=255)
                     aTTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGGT‑TTTCGCTTGGTTC‑TGCCACCGGcc                                                                                                                                   >  1:794759/1‑149 (MQ=255)
                             cAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCcat                                                                                                                        <  1:506076/151‑1 (MQ=255)
                              aaTTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGGTA‑TTTTGTTTGGTTC‑GGCCATCCGCCTCCT‑CGCcat                                                                                                                        >  1:571802/1‑150 (MQ=255)
                               aTTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCcattg                                                                                                                      <  2:770761/151‑1 (MQ=255)
                                ttAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGTTGTATA‑TTTTGTTTTGGTT‑GGCCCTAAGCCTCTT‑CGcc                                                                                                                          >  2:504697/1‑146 (MQ=255)
                                     ccGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATA‑TTTTGTTTTGGTC‑GACCATCCGCCCCTG‑CGCcatttcatt                                                                                                                  >  2:605082/1‑149 (MQ=255)
                                     ccGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGAGCCATT‑TTTTTTTTCGTTC‑GGACATCAGGCTCTG‑GGcccatggcttt                                                                                                                 >  1:257650/1‑141 (MQ=255)
                                        cgAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTTGATA‑TTTTCTTTTGTTT‑GACCAACTACCCCCG‑CGCCCCTTCATTGGGc                                                                                                              >  2:326315/1‑150 (MQ=255)
                                            tGTTAGATCTCAATATGTTTCTGTTCATCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTTGTTTTGTTC‑GACCCTCCGCCCCTG‑CGCCAATTCATTTGGCGGACt                                                                                                         >  2:813313/1‑151 (MQ=255)
                                                 gATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATA‑ATTTTTTTGGTTT‑GGCCCTCGGCCTCTG‑CGCCCCTGCATTTGGGGGGCTGtttt                                                                                                    >  2:453898/1‑150 (MQ=255)
                                                     tcAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGAGCGATA‑TTTTGTTTTGTTC‑GACCATCGGCCTCTG‑CGCCATTGTATTTGGCGTACTGgttt                                                                                                    >  1:99445/1‑146 (MQ=255)
                                                      cAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAg                                                                                                <  2:771679/150‑1 (MQ=255)
                                                           tGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTTGTTTTGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTTGGCGGACTGGTTCAAAatggtg                                                                                           >  2:212887/1‑150 (MQ=255)
                                                               gCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTTGTTTTGTTT‑GACCCTCGGCCTCTG‑CGCCCTTTCATTTGGGGGACTGTTTTAAAgtgggg                                                                                           >  2:837184/1‑144 (MQ=255)
                                                                     cAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTTGTTTTGTTT‑GACCATCGGCCTCTG‑CGCCCTTTCATTGGGGGGACTGTTTTAAAGTGGGGTGCCGATGcc                                                                                 >  2:622234/1‑150 (MQ=255)
                                                                        cACCTGTGACTTGATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATTGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCg                                                                                <  1:33595/148‑1 (MQ=255)
                                                                         aCCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATTGGTT‑TTTCCTTTGGGTT‑GCCCACCGCCCCCCC‑CGCCATTGTACTGGGGGGTCTTTTTCTCACCGGTGTGTCGCTGcc                                                                                 >  1:348296/1‑146 (MQ=255)
                                                                          ccAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTgg                                                                            <  2:723487/150‑1 (MQ=255)
                                                                           cAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGGTA‑TTTCGTTTGGGTT‑TGCCAACGGGCTCCT‑TGCCCTTTCATTTGGGGGAGTGTTTTAAAGTGGTGTGGCGGTGGCGCTGGCg                                                                          >  2:857441/1‑151 (MQ=255)
                                                                           cAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCg                                                                          <  1:94389/151‑1 (MQ=255)
                                                                                aCTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACt                                                                     <  2:733606/151‑1 (MQ=255)
                                                                                aCTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACt                                                                     <  2:555877/151‑1 (MQ=255)
                                                                                    tATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCt                                                                  <  2:775060/150‑1 (MQ=255)
                                                                                         gcattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTA‑TTTCGTTTGGTTC‑GCCCCTCCGCCCCTC‑CGCCATTGCATTGGGCGGGCTGTTTTAAAAAGGAGAGCCGATGCCGCTGGCGAGACCCCTGACCtt                                                            >  2:127239/1‑151 (MQ=255)
                                                                                          cattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTTGTTTTGTTC‑GCCCATCCGCCTCTG‑CGCCCTTGCATTGTGGGGGCTGTTTCTAAGTGGGGTGCCGCTGCCGCTGGCGCTACTCCTGACCtt                                                            >  1:509235/1‑150 (MQ=255)
                                                                                             tggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGCTTGGTT‑TTTTCTTTGGTTT‑GACCACCCGCCTCCG‑CGGCCTTTCCTTGGGGGGTCTGTTTCATAGTGGGGTGCTGGTGCGGCTGGCGAGTATCCTGGCCTTGCTg                                                        >  1:31696/1‑151 (MQ=255)
                                                                                                cattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTc                                                     <  2:244021/151‑1 (MQ=255)
                                                                                                 attgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGGTA‑TTTTGTTTGGTTC‑TACCCTCGGGCTCTG‑CGCCCTTGCCTTGGGGGGACTGTTTTTAAGGGGGGGGCCCGTGCCGCTGGCGGTACTCCTCGCCTTACTGCTc                                                     >  2:124820/1‑150 (MQ=255)
                                                                                                 attgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATTGATA‑TTTCGTTTTGTTCTCACCA‑CCGCCTCTG‑CGCCCTTGGCTTGGGCGGTCTGGTTTAAAGGTGTGGGCCGATGGCGCTGGCGGTACTCCTGAcccttc                                                          >  1:667150/1‑141 (MQ=255)
                                                                                                     tCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTTGTTTTGGTT‑GACCCTCGGCCCCTG‑CGCCCTTGCCTTTGGGGTACTGTTTTAAAATGGGGTGCCGATGCCGCTGGCGGTACTCCCGACCTTACTGCTCGgcg                                                 >  2:518268/1‑150 (MQ=255)
                                                                                                           tACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTATTTTGTTTC‑GACCA‑CCGCCCCTG‑CGCCCATGCCTTTGGGGTGCTTTTTTTAAGTGGGGTGCCGATGCCGCCGGCGATACTCCTGTACTTTCTGCTCGGCGCact                                             >  1:345439/1‑146 (MQ=255)
                                                                                                             ccGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTTGTTTTGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTTGGCGTACTGTTTTAAAGTGGGGTGTCGATGCCGCTGGCGATAATCCTGACCCTACTGCTCGGCGCCTTGTGCg                                        >  2:277158/1‑151 (MQ=255)
                                                                                                                 tAACGATGGTGATTGTCAGTGGCGGGAGCGGGA‑TTTCCTTTTGTTC‑CGGCACCAGCCCCTGCCGCCTTTCCATGTGCGGTCTTCTTTCAAAG‑GGCGTGCCCGTGCCGCCCGCGAAACCCCTTACCCTACTGCTCCGCGCCtttgtcgggg                                     >  1:131456/1‑142 (MQ=255)
                                                                                                                  aaCGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTTGTTTTGTTC‑GACCATCGGCCCCTG‑CGCCATTGCATTGGGCGGACTGTTTTAAAGTGGTGGGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGGGCATTGTGCggg                                      >  2:508228/1‑148 (MQ=255)
                                                                                                                  aaCGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCt                                    <  2:634934/150‑1 (MQ=255)
                                                                                                                    cGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTTTTTTTGGTC‑GGCCATCGGCCTCTT‑CGCCATTGCATTTGGGGTGCGGTTTTAAAGTGGGGGGCCGATGCCGGTGGGGGTACTCCTGACCTTTCTGCTCGGCGGAGTGTGggggggg                                   >  1:65897/1‑146 (MQ=255)
                                                                                                                      aTGGTGATTGTCAGTGGCGGGATCGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCa                               <  2:137982/151‑1 (MQ=255)
                                                                                                                      aTGGTGATTGTCAGTGGCGGGATCGATA‑TGTCGGTTGGTTC‑GGCCATCGGCCTCTG‑CGCCTTTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCa                               <  1:92601/151‑1 (MQ=255)
                                                                                                                           gATTGTCAGTGGCGGGATCGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGcc                          <  1:721887/151‑1 (MQ=255)
                                                                                                                              tGTCAGTGGCGGGGACTATA‑TTTTGTTTGGTTC‑GGTCAACCGCCTCTG‑CGCCCTTGTCTTGGGGGGACTGTTTTAAAGGGGGGTGGCGATGCCGCTGGCGATACTCCTGACCTTTCTTCTCGGGGCATTGTGTGGGCTGGTCAAACCCgg                        >  1:215418/1‑150 (MQ=255)
                                                                                                                               gTCAGTGGCGGGATCGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGat                      <  2:270586/151‑1 (MQ=255)
                                                                                                                                     ggCGGGATCGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGattaatta                <  2:521262/151‑1 (MQ=255)
                                                                                                                                     ggCGGGATCGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGattaatt                 <  1:3331/150‑1 (MQ=255)
                                                                                                                                        gggATCGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCt             >  1:248385/1‑151 (MQ=255)
                                                                                                                                             cGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATAc         <  2:273195/150‑1 (MQ=255)
                                                                                                                                                ta‑tTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCaa      >  2:463686/1‑150 (MQ=255)
                                                                                                                                                 a‑tTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCaa      >  2:772509/1‑149 (MQ=255)
                                                                                                                                                     tCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGtt  >  2:319091/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GCTCGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATA‑TTTCGTTTGGTTC‑GACCATCGGCCTCTG‑CGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTT  >  NC_000913/1604088‑1604385

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