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

CGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAA  >  NC_000913/1604091‑1604387
                                                                                                                                                    |                                                                                                                                                       
cGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGt                                                                                                                                                       <  2:935954/149‑1 (MQ=255)
 gTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGttt                                                                                                                                                     >  1:1074364/1‑150 (MQ=255)
  ttATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGttt                                                                                                                                                     >  2:557464/1‑149 (MQ=255)
    aTTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTttgttttgtt                                                                                                                                                 >  1:641495/1‑151 (MQ=255)
    aTTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGtt                                                                                                                                                 <  2:623218/151‑1 (MQ=255)
     ttGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTttgttttgtt                                                                                                                                                 >  2:647928/1‑150 (MQ=255)
     ttGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTttgttttgtt                                                                                                                                                 >  2:93342/1‑150 (MQ=255)
       gagaTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCg                                                                                                                                               <  2:368309/150‑1 (MQ=255)
        agaTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGAc                                                                                                                                             <  2:326643/151‑1 (MQ=255)
         gaTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGTTTGTTATTTCTTTTGGtttggcc                                                                                                                                            >  1:94462/1‑146 (MQ=255)
            tGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTTGGTCTTTCTTTTTGTTTGACCCTc                                                                                                                                         >  1:804759/1‑151 (MQ=255)
             gTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGGCCATCg                                                                                                                                        <  2:108975/151‑1 (MQ=255)
                        tGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATTGGTATTTTTTTTGGGTCGGCCATCGGGCCCTgcg                                                                                                                               >  1:545233/1‑149 (MQ=255)
                        tGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTCGACCATCGGCCTCTGCGcc                                                                                                                             >  2:377915/1‑151 (MQ=255)
                            aTTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCcattg                                                                                                                         <  2:100364/151‑1 (MQ=255)
                              tAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCcattgc                                                                                                                        <  1:23004/150‑1 (MQ=255)
                               aaCCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTTGTTTTGTTTGACCATCCGCCCCTGCGccccttca                                                                                                                       >  1:664683/1‑144 (MQ=255)
                               aaCCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCcattgcat                                                                                                                      <  2:287366/151‑1 (MQ=255)
                                    gcgAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGCccttgcatttg                                                                                                                   >  2:758514/1‑147 (MQ=255)
                                    gcgAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGc                                                                                                                 <  2:899750/151‑1 (MQ=255)
                                        aTGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTTGATATTTTTTTTTGTTTGACCATCCGCCTCTGCGCCCTTGCATTTGGCg                                                                                                                >  1:1062532/1‑148 (MQ=255)
                                         tGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTTGGCCAACCGCCTCTTGGCCATTTCATTGGGCGTaat                                                                                                            >  2:858703/1‑149 (MQ=255)
                                          gTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACt                                                                                                            <  2:541645/150‑1 (MQ=255)
                                            tAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGt                                                                                                          <  2:763128/150‑1 (MQ=255)
                                                  tcAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGGTATTTTGTTTGGTTTGACCATCGGCCTCCGCGCCCCTGCATTTGGGGGGCTGTTTcca                                                                                                     >  1:960026/1‑147 (MQ=255)
                                                  tcAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTTGACCATCGGCCTCTGCGCCCTTTCCTTGGGCGTACTGTTTTaaa                                                                                                    >  2:458323/1‑150 (MQ=255)
                                                    aaTATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAgtg                                                                                                 <  2:904517/151‑1 (MQ=255)
                                                      tatGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTTGTTTGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAgtgg                                                                                                >  2:484192/1‑150 (MQ=255)
                                                       atGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAgtggtg                                                                                              <  2:1067421/151‑1 (MQ=255)
                                                         gTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTCGGGTCGACCATCCGCCTCTGCGCCCTTTCATTGGGCGGACTGTTTTTAAATGGggtg                                                                                            >  2:653330/1‑151 (MQ=255)
                                                         gTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGtgtg                                                                                            <  2:788172/151‑1 (MQ=255)
                                                         gTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGtgtg                                                                                            <  2:32237/151‑1 (MQ=255)
                                                         gTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGCGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCCTTGGGCGTACTGTTTTAAAGTGGtgtg                                                                                            >  1:169987/1‑151 (MQ=255)
                                                          ttGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGGTTGACCATCGGCCTCTGCGCCCTTGCCTTTGGCGGACTGTTTTAAAATGGtgtg                                                                                            >  2:721963/1‑150 (MQ=255)
                                                             cTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCCTTATTTTGCTTTGTTCGGCCCTCGGCCCCTGCCCCCTTTCATTGGGCGGACTCTTTTaaaaggg                                                                                                >  1:230986/1‑139 (MQ=255)
                                                                 tCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTTGTTTTGTTCGCCCATCCGCCTCTGCGCCATTGCATTGGGCGTACTGTTTTAAAGTGGGGTGGCgcggc                                                                                     >  1:531359/1‑146 (MQ=255)
                                                                  cAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGcc                                                                                   <  2:388795/150‑1 (MQ=255)
                                                                      aCCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTTGTTATTTTGTTTTGTTTGACCATCGGCCCCTGCGCCCCTGCATTGGGGGGACTGGTTCTAAGGGGGGTGcc                                                                                          >  1:582448/1‑140 (MQ=255)
                                                                        cAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCg                                                                            <  2:395810/151‑1 (MQ=255)
                                                                              cTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCTGCCATCCGCCTCTGCGGCCTTGCATTGGGCGGTCTGTTTTTAAGGGGTGTGCCCGTGCCGCTGGGGATAc                                                                        >  1:86866/1‑149 (MQ=255)
                                                                              cTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCCCTGCGCCATTGCATTGGGGGGACTGTTTTCAAGAGGTGTGCCGATGCCGCTGGCGGTACTc                                                                      >  2:91192/1‑151 (MQ=255)
                                                                              cTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCCGCCCCTGCGCCATTGCATTGGGGGTACTGTTTTAAAGTGGTGTGCCGATGCCGCTGGCGGTACt                                                                       >  1:496091/1‑150 (MQ=255)
                                                                                      gcattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGa                                                                  <  1:721919/147‑1 (MQ=255)
                                                                                      gcattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCt                                                               <  2:38487/150‑1 (MQ=255)
                                                                                       cattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTCGTTTTGTTCGACCATCGGCCTCTGCGCCCTTGCATTGGGCGGACTGTTTTTAAGTGGGGTGCCGGTGCCGCTGGGGGTACTCCTGACCtt                                                              >  1:371709/1‑150 (MQ=255)
                                                                                       cattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTa                                                             <  1:624662/151‑1 (MQ=255)
                                                                                          tggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGTTTTTTTCTTTTGGTTGACCCACCGCCCCCTCGCCCTTTCCTTGGGGGGACTTTTTTTAAGGGGTGTGCCGATGCCGcctgcggttatcc                                                                     >  1:814207/1‑128 (MQ=255)
                                                                                           ggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGACGATATTTTGTTTGGTTTGACCAACCGCCCCTTCGCCCTTTCCTTTTGCGGTCTGGTTCCAACGGGGGTGCCGCTGGCCGCTGGGATCCTCCTTGCCCTACtac                                                         >  2:117145/1‑149 (MQ=255)
                                                                                           ggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTg                                                          <  2:265026/150‑1 (MQ=255)
                                                                                              attgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGACGGTAATTTTTTTTGGTTGACCATCCGCCCCTGCGCCACTGCCTTTTGGGGACTGGTTCAAAACGGTGTGGCGATGGCCCGGGGGAGACTCCCGTTCTTACTGc                                                         >  2:260322/1‑148 (MQ=255)
                                                                                                 gTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTCTTTTGTTCGCCCATCGGCCTCTGCGCCCCTTCCTTGGGCGGACTGTTTTTAAGTGGGGGGCCCGTGCCCCCCGCGGGACTCCCGGCCTTTCTGCTCgg                                                     >  1:604903/1‑149 (MQ=255)
                                                                                                      ccTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTCCTTTTGTTTTACCACCGGCCTCTGCGCCCTTGCCTTTGGCGGTCTGTTTTTAAGAGGGGTGCCGCTGCCCCCGGCGATGCTCCTGACCTTACTg                                                          >  1:566674/1‑139 (MQ=255)
                                                                                                        tACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGGACTGTTTTAAAGTGGGGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATtgtg                                            >  2:476880/1‑151 (MQ=255)
                                                                                                              tAACGATGGTGATTGTCAGTGGCGGGAGCGATATTTTGTTTTGTTCGACCATCGGCCCCTTCGCCCTTGCATTTGGGGTTCTGTTTTTAAGTGGGGTGGCCCTGCCGCTGGCGGTACCCCTGGCCCTTCTGCTCGGCGCATTGTGgggg                                        >  1:221506/1‑149 (MQ=255)
                                                                                                               aaCGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTTGTTTGACCATCCGCCCCTGCGCCATTTCATTTGGCGTACTGTTTCAAAGTGGTGTGGCGATGCCGCTGGCGGTACTCCTTACCTTAATGCTCGGCGCCTTTTGCggg                                        >  2:1008115/1‑148 (MQ=255)
                                                                                                                      gTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACg                              <  2:641536/151‑1 (MQ=255)
                                                                                                                      gTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACg                              <  1:911949/151‑1 (MQ=255)
                                                                                                                       tGATTGTCAGTGGCGGGATCGGAATTTTGTTTTGTTTGACCATCGGCCTCTGCGCCCTTGCATTGGGGGGACTGTTTTAAAGGGGGGTGGCGATGCCGCCGGCGATACTCCTGACCTTTCTGCTCGGCGCCTTGTGGGGGGTGGTCAAccc                             >  1:982297/1‑149 (MQ=255)
                                                                                                                                          cGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATAcc          <  2:34966/151‑1 (MQ=255)
                                                                                                                                              tatTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCa         >  1:150043/1‑149 (MQ=255)
                                                                                                                                                atTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCaaa       <  1:607005/150‑1 (MQ=255)
                                                                                                                                                    cGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTaa  <  2:400450/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
CGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAA  >  NC_000913/1604091‑1604387

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