Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,604,2370CT46.0% 65.4 / 30.2 64S64L (TCG→TTG) lsrDautoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand):  ref base C (8/26);  new base T (29/0);  total (38/26)
Fisher's exact test for biased strand distribution p-value = 5.31e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAA  >  NC_000913/1604091‑1604387
                                                                                                                                                   |                                                                                                                                                        
cGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGttt                                                                                                                                                      >  1:519483/1‑151 (MQ=255)
 gTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGGCCATATTTTGttt                                                                                                                                                      >  1:74981/1‑150 (MQ=255)
  ttATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGttt                                                                                                                                                      >  1:1080584/1‑149 (MQ=255)
   tATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTgg                                                                                                                                                    <  1:167872/150‑1 (MQ=255)
          aTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCa                                                                                                                                            <  1:797020/151‑1 (MQ=255)
            tGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTGGTTCGACCCTc                                                                                                                                          >  1:22180/1‑151 (MQ=255)
                 cATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTTTTTGGtttgtccctccg                                                                                                                                        >  2:430987/1‑138 (MQ=255)
                 cATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCct                                                                                                                                     <  1:1229584/151‑1 (MQ=255)
                      ggggCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGGGATTGTCAGTGGCGGGTGTGGTTATTTGTTTTGTTCGACCagctggcgcag                                                                                                                                  >  1:765349/4‑139 (MQ=255)
                       gTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTgcgc                                                                                                                               <  1:631026/151‑1 (MQ=255)
                       gTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTCATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGAATGATGTTTCCTTTTGTTCGCCCATCGGCCCCTgcg                                                                                                                                >  2:659423/1‑150 (MQ=255)
                           aaTTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTCGTTTTGGTTGACCCTCCGCCCCTTCGCcctt                                                                                                                           >  2:250470/1‑151 (MQ=255)
                              tAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCcattgc                                                                                                                         <  1:362056/150‑1 (MQ=255)
                                 cccGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCcattgcattg                                                                                                                    <  1:891532/151‑1 (MQ=255)
                                   cgcgAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTCGACCATCGGGCTCTGCGCCCTTGTATTggg                                                                                                                  >  2:1050805/1‑151 (MQ=255)
                                     cgAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGGTCGACCATCGGCCTCTGCGCCCTTGCCTTtgg                                                                                                                  >  2:671267/1‑149 (MQ=255)
                                     cgAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCg                                                                                                                <  1:120610/151‑1 (MQ=255)
                                         tGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTa                                                                                                              <  1:251894/149‑1 (MQ=255)
                                           ttAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGGTTGACCATCGGCCTCTGCGCCATTGCATTtgggggaat                                                                                                            >  2:404310/1‑143 (MQ=255)
                                             aGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGtt                                                                                                         <  1:1182599/150‑1 (MQ=255)
                                               aTCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTCTTTTTTTTGACCCTCGGCCCCCCCGCCATTGCATTGGGGGGACTGttt                                                                                                        >  1:765679/1‑149 (MQ=255)
                                                tctcAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTTGTTTGACCATCGGCCTCTTCGCCCTTGCATTTGGCGTACTGttt                                                                                                        >  2:694469/1‑148 (MQ=255)
                                                   cAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTc                                                                                                       <  1:87601/146‑1 (MQ=255)
                                                        tGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTCGTTTTGTTCGACCCTCGGCCTCTGCGCCACTGCATTGGGCGTGCTGTTTTTAAgtgg                                                                                                >  1:1061048/1‑148 (MQ=255)
                                                          ttGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTCTTTTTTTTCCACCCACCGCCTCTGGGCCCTTGCATTGGGCGCACTGTTTTaaa                                                                                                    >  1:987909/1‑142 (MQ=255)
                                                                 tCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTTGGCCATCCGCCTCTGCGCCCATGCATTGGGGGGACTGTTTTTAAGTGGtgtg                                                                                            >  1:1246293/1‑143 (MQ=255)
                                                                  cAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGGTATTTTGTTTGGTTTGCACATCGGCCTCCGCGCCCTTGCATTGGGGGGGCTGGTTTTAAGTGGTGGGCCGATGcc                                                                                   >  1:1118178/1‑150 (MQ=255)
                                                                   aGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCACCGGCCTCCGCGCCATTGCATTGGGCGTACTGTTTTAAAGTGGGGTGCCGATGCCg                                                                                  >  1:1119078/1‑150 (MQ=255)
                                                                     cacccgttgcTTTATCTGCATTGGCCTTGTCGCCCTTCCGCTAACGAAGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCg                                                                                  <  1:1153132/139‑1 (MQ=255)
                                                                     cACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCCTTGCATTTGGCGGACTGTTTTAAAGTGGTGTGCCGATGCCGc                                                                                 >  2:197644/1‑149 (MQ=255)
                                                                     cACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTCGACCATCCGCCTCTGCGCCATTGCATTGGGGGTACTGTTTCAAAGGGGGGGGCCGATGCCGccg                                                                               >  2:908620/1‑149 (MQ=255)
                                                                         aGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGa                                                                           <  1:272006/151‑1 (MQ=255)
                                                                         aGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGAAATTTTCTTTTGTTCGACCCTCGGCCTCTGCGCCACTGCATTGGGGGTACTGGTTTAAAGTGGTGGGGCGAAGCCGc                                                                                 >  2:1143042/1‑145 (MQ=255)
                                                                          gTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTTGTTTTGTTTGACCCACCGCCCCTTCGCCATTGCATTGGGGGGACTTTTTCTAAAGGGTGGGCCGGTGCCGCTGGCggt                                                                          >  1:828077/1‑149 (MQ=255)
                                                                           tGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATTGATATTTTGTTTTGTTCGACCATCCGCCTCTGCGCCCTTGCATTGGGCGGGCTGTTTTAAAGTGGGGTGCCGATGCCGCCGGCggt                                                                          >  2:359931/1‑148 (MQ=255)
                                                                           tGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTTGACCATCCGCCTCTGCGCCCTTGCCTTTGGCGGACTGTTTTCAAGAGGTGTGCCGATGCCGCTGGCg                                                                            >  2:629945/1‑148 (MQ=255)
                                                                             aCTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCg                                                                            >  2:379624/1‑146 (MQ=255)
                                                                                ttATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTcc                                                                     <  2:519484/150‑1 (MQ=255)
                                                                                ttATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTc                                                                      <  2:696907/149‑1 (MQ=255)
                                                                                  aTCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTTGACATTTTACTTCGCTCGGGCACAGGCCTCCGCCCCGCTGCATTAGGGCGTATGTTTTAAAGGGGTGTCCCGGAGCCGCCGGCGGTAcccc                                                                     >  1:1105245/1‑148 (MQ=255)
                                                                                  aTCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGGGGGCTGTTTTAAAGTGGGGTGCCGATGCCGCTGGCGGTACTcc                                                                     >  1:425934/1‑148 (MQ=255)
                                                                                   tCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTCGACCATCGGCCTCTGCGCCCTTGCATTTGGGGTACTGTTTTAAAGTGGGGTGCCGATGCCGCCGGCGATACTCCTTAc                                                                 >  2:1030412/1‑151 (MQ=255)
                                                                                      gcattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCATTGCATTTGGCGGACTGTTTTAAAGTGGGGTGCCGATGCCCCGGGCGGTACTCCTGAcc                                                                >  2:1177076/1‑149 (MQ=255)
                                                                                      gcattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCtt                                                              <  2:232831/151‑1 (MQ=255)
                                                                                        attggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTAc                                                            <  2:303308/151‑1 (MQ=255)
                                                                                          tggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACt                                                           <  2:332606/150‑1 (MQ=255)
                                                                                            gcattgTCGCCCTACCGCTAACGCTGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGGTCGAACAACCGCCTCCGCGCCCTTGCATTGGGCGGGCTGGTTCCAAGAGGGGTGCCCATGCCGCCGGCGGTACTCCTGGCCTTAc                                                            >  1:1181550/1‑147 (MQ=255)
                                                                                            gcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCt                                                        <  2:286854/151‑1 (MQ=255)
                                                                                              attgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGGTTTGTTCGACCATCGGCCCCTGCGCCCTTTCATTGGGGGGGCTGTTTCTAAGTGGGGTGGCGGTGCCGCCGGCGAGAATCCTGGCCTTACTGc                                                         >  2:1240608/1‑148 (MQ=255)
                                                                                                tgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCgg                                                     <  1:44999/150‑1 (MQ=255)
                                                                                                   cGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGgcgc                                                  <  1:1070819/150‑1 (MQ=255)
                                                                                                   cGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCa                                                 <  1:838233/151‑1 (MQ=255)
                                                                                                    gCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTTGTTCTACCATCCGCCCCTGCCCCCTTTCATTGGGGGTGCTGTTTTAAAGAGGGGTGCCGCTGCCCCCGGCGGTACTCCCGACCTTTCTGCTCcgcggct                                                >  1:1228958/1‑148 (MQ=255)
                                                                                                    gCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTTGTTTTGTTTGACCAACGGCCTCTGGGCCATTTTCATGGGCGGACTGGTTTAAAATGGGGTGCCCATGGCCGCGGCGATTACCCTGGCCCTACTGGTCgg                                                     >  1:825797/1‑146 (MQ=255)
                                                                                                        tACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTCGACCAACGGGCCCCGCGCCATTGCGTTGGGCGTACTTTTTTTAAGAGGGGGGCCCATGCCGCCGGCGGTACTCCTGCCCTTTCTGCTCGGCGCGtt                                               >  1:367442/1‑148 (MQ=255)
                                                                                                         aCCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATtgtg                                            <  1:853814/150‑1 (MQ=255)
                                                                                                            gCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTCGACCATCGGCCTCTTCGGCCTTTCATTTGGCGGACTGTTTTAAAGGGGGGGGCCGGTGCCGCTGGCGGTACTCCTGACCTTACTGCTCGGCGGCTTGTGCggg                                        >  1:1083682/1‑151 (MQ=255)
                                                                                                               aaCGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTg                                     <  1:213600/151‑1 (MQ=255)
                                                                                                                      gTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGGTTGACCAACGGCCTCTGCGCCATTTGATTTGGGGTACTGTTTTAAAGTGGTGTGTCGCTGCCCCTGGCGGTACTCCTGACCTTACTGCGCGGCGCATTTTGgggg                                        >  1:1010987/1‑141 (MQ=255)
                                                                                                                          ttGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCg                           <  1:1237686/150‑1 (MQ=255)
                                                                                                                              cAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGatt                       <  2:748068/150‑1 (MQ=255)
                                                                                                                                     gggATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCt               <  2:237780/151‑1 (MQ=255)
                                                                                                                                      ggATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCt               >  2:52253/1‑150 (MQ=255)
                                                                                                                                             atatTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCa         >  2:568221/1‑150 (MQ=255)
                                                                                                                                                   cGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTaa  <  1:525062/151‑1 (MQ=255)
                                                                                                                                                   cGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTa   >  2:380999/1‑150 (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 < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: