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

CGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAA  >  NC_000913/1604091‑1604387
                                                                                                                                                   |                                                                                                                                                        
cGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTTGttt                                                                                                                                                      >  1:519483/1‑151 (MQ=255)
 gTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑GGCCATATTTTGttt                                                                                                                                                      >  1:74981/1‑150 (MQ=255)
  ttATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTTGttt                                                                                                                                                      >  1:1080584/1‑149 (MQ=255)
   tATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTgg                                                                                                                                                    <  1:167872/150‑1 (MQ=255)
          aTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCa                                                                                                                                            <  1:797020/151‑1 (MQ=255)
            tGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑GTCGATATTTTGTTTGGTTCGACCCTc                                                                                                                                          >  1:22180/1‑151 (MQ=255)
                 cATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTTTTTTGGtttgtccctccg                                                                                                                                        >  2:430987/1‑138 (MQ=255)
                 cATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCct                                                                                                                                     <  1:1229584/151‑1 (MQ=255)
                      ggggCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGGGATTGTCAGTGGCGGGTGTGGTTATTT‑GTTTTGTTCGACCagctggcgcag                                                                                                                                  >  1:765349/4‑139 (MQ=255)
                       gTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTgcgc                                                                                                                               <  1:631026/151‑1 (MQ=255)
                       gTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTCATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑AATGATGTTTCCTTTTGTTCGCCCATCGGCCCCTgcg                                                                                                                                >  2:659423/1‑150 (MQ=255)
                           aaTTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGGTATTTCGTTTTGGTTGACCCTCCGCCCCTTCGCcctt                                                                                                                           >  2:250470/1‑151 (MQ=255)
                              tAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCcattgc                                                                                                                         <  1:362056/150‑1 (MQ=255)
                                 cccGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCcattgcatt‑g                                                                                                                    <  1:891532/151‑1 (MQ=255)
                                   cgcgAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTTGTTTGGTTCGACCATCGGGCTCTGCGCCCTTGTATT‑ggg                                                                                                                  >  2:1050805/1‑151 (MQ=255)
                                     cgAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTTGTTTTGGTCGACCATCGGCCTCTGCGCCCTTGCCTT‑tgg                                                                                                                  >  2:671267/1‑149 (MQ=255)
                                     cgAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCg                                                                                                                <  1:120610/151‑1 (MQ=255)
                                         tGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTa                                                                                                              <  1:251894/149‑1 (MQ=255)
                                           ttAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTTGTTTTGGTTGACCATCGGCCTCTGCGCCATTGCATT‑tgggggaat                                                                                                            >  2:404310/1‑143 (MQ=255)
                                             aGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGtt                                                                                                         <  1:1182599/150‑1 (MQ=255)
                                               aTCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTTCTTTTTTTTGACCCTCGGCCCCCCCGCCATTGCATT‑GGGGGGACTGttt                                                                                                        >  1:765679/1‑149 (MQ=255)
                                                tctcAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTTGTTTGACCATCGGCCTCTTCGCCCTTGCATT‑TGGCGTACTGttt                                                                                                        >  2:694469/1‑148 (MQ=255)
                                                   cAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTc                                                                                                       <  1:87601/146‑1 (MQ=255)
                                                        tGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑GTCGATATTTCGTTTTGTTCGACCCTCGGCCTCTGCGCCACTGCATT‑GGGCGTGCTGTTTTTAAgtgg                                                                                                >  1:1061048/1‑148 (MQ=255)
                                                          ttGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGGTATTTCTTTTTTTTCCACCCACCGCCTCTGGGCCCTTGCATT‑GGGCGCACTGTTTTaaa                                                                                                    >  1:987909/1‑142 (MQ=255)
                                                                 tCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTTGTTTTGTTTGGCCATCCGCCTCTGCGCCCATGCATT‑GGGGGGACTGTTTTTAAGTGGtgtg                                                                                            >  1:1246293/1‑143 (MQ=255)
                                                                  cAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑GTCGGTATTTTGTTTGGTTTGCACATCGGCCTCCGCGCCCTTGCATT‑GGGGGGGCTGGTTTTAAGTGGTGGGCCGAT‑Gcc                                                                                   >  1:1118178/1‑150 (MQ=255)
                                                                   aGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTTGTTTTGTTCGACCACCGGCCTCCGCGCCATTGCATT‑GGGCGTACTGTTTTAAAGTGGGGTGCCGAT‑GCCg                                                                                  >  1:1119078/1‑150 (MQ=255)
                                                                     cacccgttgcTTTATCTGCATTGGCCTTGTCGCCCTTCCGCTAACGAAGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCg                                                                                  <  1:1153132/139‑1 (MQ=255)
                                                                     cACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCCTTGCATT‑TGGCGGACTGTTTTAAAGTGGTGTGCCGAT‑GCCGc                                                                                 >  2:197644/1‑149 (MQ=255)
                                                                     cACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTTGTTTGGTTCGACCATCCGCCTCTGCGCCATTGCATT‑GGGGGTACTGTTTCAAAGGGGGGGGCCGAT‑GCCGccg                                                                               >  2:908620/1‑149 (MQ=255)
                                                                         aGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGa                                                                           <  1:272006/151‑1 (MQ=255)
                                                                         aGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGAAATTTTCTTTTGTTCGACCCTCGGCCTCTGCGCCACTGCATT‑GGGGGTACTGGTTTAAAGTGGTGGGGCGAA‑GCCGc                                                                                 >  2:1143042/1‑145 (MQ=255)
                                                                          gTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGGTATTTTGTTTTGTTTGACCCACCGCCCCTTCGCCATTGCATT‑GGGGGGACTTTTTCTAAAGGGTGGGCCGGT‑GCCGCTGGCggt                                                                          >  1:828077/1‑149 (MQ=255)
                                                                           tGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATTGATATTTTGTTTTGTTCGACCATCCGCCTCTGCGCCCTTGCATT‑GGGCGGGCTGTTTTAAAGTGGGGTGCCGAT‑GCCGCCGGCggt                                                                          >  2:359931/1‑148 (MQ=255)
                                                                           tGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTTGTTTTGTTTGACCATCCGCCTCTGCGCCCTTGCCTT‑TGGCGGACTGTTTTCAAGAGGTGTGCCGAT‑GCCGCTGGCg                                                                            >  2:629945/1‑148 (MQ=255)
                                                                             aCTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCg                                                                            >  2:379624/1‑146 (MQ=255)
                                                                                ttATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTcc                                                                     <  2:519484/150‑1 (MQ=255)
                                                                                ttATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTc                                                                      <  2:696907/149‑1 (MQ=255)
                                                                                  aTCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑GTTGACATTTTACTTCGCTCGGGCACAGGCCTCCGCCCCGCTGCATTAGGGCGTA‑TGTTTTAAAGGGGTGTCCCGGA‑GCCGCCGGCGGTAcccc                                                                     >  1:1105245/1‑148 (MQ=255)
                                                                                  aTCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGGGGGCTGTTTTAAAGTGGGGTGCCGAT‑GCCGCTGGCGGTACTcc                                                                     >  1:425934/1‑148 (MQ=255)
                                                                                   tCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTTGTTTGGTTCGACCATCGGCCTCTGCGCCCTTGCATT‑TGGGGTACTGTTTTAAAGTGGGGTGCCGAT‑GCCGCCGGCGATACTCCTTAc                                                                 >  2:1030412/1‑151 (MQ=255)
                                                                                      gcattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCATTGCATT‑TGGCGGACTGTTTTAAAGTGGGGTGCCGAT‑GCCCCGGGCGGTACTCCTGAcc                                                                >  2:1177076/1‑149 (MQ=255)
                                                                                      gcattggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTCCTGACCtt                                                              <  2:232831/151‑1 (MQ=255)
                                                                                        attggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTCCTGACCTTAc                                                            <  2:303308/151‑1 (MQ=255)
                                                                                          tggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTCCTGACCTTACt                                                           <  2:332606/150‑1 (MQ=255)
                                                                                            gcattgTCGCCCTACCGCTAACGCTGGTGATTGTCAGTGGCGGG‑ATCGATATTTTGTTTGGGTCGAACAACCGCCTCCGCGCCCTTGCATT‑GGGCGGGCTGGTTCCAAGAGGGGTGCCCAT‑GCCGCCGGCGGTACTCCTGGCCTTAc                                                            >  1:1181550/1‑147 (MQ=255)
                                                                                            gcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTCCTGACCTTACTGCt                                                        <  2:286854/151‑1 (MQ=255)
                                                                                              attgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTTGGTTTGTTCGACCATCGGCCCCTGCGCCCTTTCATT‑GGGGGGGCTGTTTCTAAGTGGGGTGGCGGT‑GCCGCCGGCGAGAATCCTGGCCTTACTGc                                                         >  2:1240608/1‑148 (MQ=255)
                                                                                                tgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTCCTGACCTTACTGCTCgg                                                     <  1:44999/150‑1 (MQ=255)
                                                                                                   cGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTCCTGACCTTACTGCTCGgcgc                                                  <  1:1070819/150‑1 (MQ=255)
                                                                                                   cGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCa                                                 <  1:838233/151‑1 (MQ=255)
                                                                                                    gCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑GTCGATATTTTGTTTTGTTCTACCATCCGCCCCTGCCCCCTTTCATT‑GGGGGTGCTGTTTTAAAGAGGGGTGCCGCT‑GCCCCCGGCGGTACTCCCGACCTTTCTGCTCcgcggct                                                >  1:1228958/1‑148 (MQ=255)
                                                                                                    gCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGGTATTTTGTTTTGTTTGACCAACGGCCTCTGGGCCATTTTCAT‑GGGCGGACTGGTTTAAAATGGGGTGCCCATGGCCGC‑GGCGATTACCCTGGCCCTACTGGTCgg                                                     >  1:825797/1‑146 (MQ=255)
                                                                                                        tACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTTGTTTGGTTCGACCAACGGGCCCCGCGCCATTGCGTT‑GGGCGTACTTTTTTTAAGAGGGGGGCCCAT‑GCCGCCGGCGGTACTCCTGCCCTTTCTGCTCGGCGCGtt                                               >  1:367442/1‑148 (MQ=255)
                                                                                                         aCCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATtgtg                                            <  1:853814/150‑1 (MQ=255)
                                                                                                            gCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTTGTTTGGTTCGACCATCGGCCTCTTCGGCCTTTCATT‑TGGCGGACTGTTTTAAAGGGGGGGGCCGGT‑GCCGCTGGCGGTACTCCTGACCTTACTGCTCGGCGGCTTGTGCggg                                        >  1:1083682/1‑151 (MQ=255)
                                                                                                               aaCGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTg                                     <  1:213600/151‑1 (MQ=255)
                                                                                                                      gTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGGTTGACCAACGGCCTCTGCGCCATTTGATT‑TGGGGTACTGTTTTAAAGTGGTGTGTCGCT‑GCCCCTGGCGGTACTCCTGACCTTACTGCGCGGCGCATTTTGgggg                                        >  1:1010987/1‑141 (MQ=255)
                                                                                                                          ttGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCg                           <  1:1237686/150‑1 (MQ=255)
                                                                                                                              cAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGatt                       <  2:748068/150‑1 (MQ=255)
                                                                                                                                     ggg‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCt               <  2:237780/151‑1 (MQ=255)
                                                                                                                                      gg‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCt               >  2:52253/1‑150 (MQ=255)
                                                                                                                                             atatTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCa         >  2:568221/1‑150 (MQ=255)
                                                                                                                                                   cGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTaa  <  1:525062/151‑1 (MQ=255)
                                                                                                                                                   cGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTa   >  2:380999/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                        
CGTTATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGG‑ATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATT‑GGGCGTACTGTTTCAAAGTGGTGTGCCGAT‑GCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTTAA  >  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: ‑