Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,640,9810CT48.2% 46.0 / 35.6 56F40F (TTC→TTTdtpBdipeptide and tripeptide permease B
Reads supporting (aligned to +/- strand):  ref base C (3/26);  new base T (27/0);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 5.54e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTT  >  NC_000913/3640832‑3641128
                                                                                                                                                     |                                                                                                                                                       
tttCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTttc                                                                                                                                                         <  2:1062826/150‑1 (MQ=255)
tttCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTttc                                                                                                                                                         <  1:642444/150‑1 (MQ=255)
   cAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGttttttt                                                                                                                                                     >  2:195616/1‑149 (MQ=255)
    aGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTttctt                                                                                                                                                      >  1:850997/1‑148 (MQ=255)
       aTGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGttt                                                                                                                                                >  1:1047161/1‑149 (MQ=255)
            aTTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGtttttttttttaa                                                                                                                                               >  2:847777/1‑140 (MQ=255)
              tGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCt                                                                                                                                        <  1:440916/150‑1 (MQ=255)
               ggAATGTTAATAAATATGAATACAACAACACCCATGGGTATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCtt                                                                                                                                       <  1:772494/150‑1 (MQ=255)
               ggAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTtttttttttttaaaaacttg                                                                                                                                        >  2:431403/1‑137 (MQ=255)
                gAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTgg                                                                                                                                     <  2:367096/151‑1 (MQ=255)
                        ataaataTGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTAAAAAGCTTTGAtttt                                                                                                                                >  1:916949/1‑146 (MQ=255)
                         taaataTGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGATTTTTTTTTTTAAAATCTTGGCtt                                                                                                                                  >  2:967765/1‑145 (MQ=255)
                          aaataTGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCaa                                                                                                                           <  1:235211/151‑1 (MQ=255)
                          aaataTGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCa                                                                                                                            <  2:692147/150‑1 (MQ=255)
                           aataTGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTTTTTAAAAGCTTGGGTTTTCGCaa                                                                                                                           >  2:568255/1‑150 (MQ=255)
                           aataTGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCGTTTAACAGCTTGGGTTTTTGCaa                                                                                                                           >  1:699470/1‑150 (MQ=255)
                            atatGAATACAACAACACCCATGGGGATGCTGTAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAaga                                                                                                                         <  2:401919/151‑1 (MQ=255)
                             tatGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTAAACAACTTGGGTTCTCGcca                                                                                                                           >  1:642543/1‑146 (MQ=255)
                               tGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAAAAGCCTGGGTTCTCGCAagag                                                                                                                        >  1:1102659/1‑149 (MQ=255)
                                 aaTACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAgg                                                                                                                    <  1:17399/151‑1 (MQ=255)
                                  aTACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTAGGCTACTACGGCGTGCAGGGCGTAATGGCGGTTTTTTTCTTTAAAAAGCTTAGGTTATCGCGAGAGCCgg                                                                                                                    >  1:395188/1‑150 (MQ=255)
                                     caacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGTGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCtt                                                                                                                 <  2:736829/150‑1 (MQ=255)
                                      aacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCGTTTAAAAATGTTGATTCTCGCCAGAGGAGGGtttt                                                                                                               >  2:1153251/1‑151 (MQ=255)
                                        caacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTCTTTCCTTAAAAACTTGGATTCTCGCCACAGAAGGCgttttt                                                                                                             >  2:597764/1‑149 (MQ=255)
                                            acCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTAAACAACTTGGGTTTTCGCAAGAGCAGGCtttttt                                                                                                             >  1:319989/1‑145 (MQ=255)
                                              ccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTTAAAAACCTTGGTTCTCGCCAGAGCAGGGTTTTTTTACttt                                                                                                       >  1:754271/1‑151 (MQ=255)
                                              ccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTGTTTAAAAGCTTTGGTTTTCGCAAGAGCAGGGGTTTTTCACttt                                                                                                       >  1:386616/1‑151 (MQ=255)
                                                aTGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTAAAAAGCCTGGAATCTCGCAAGAGAAGGCttttttttccttt                                                                                                      >  1:986983/1‑141 (MQ=255)
                                                aTGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGATTTTTTTGTTAAAAACGCTGCGTTTTCGCAAGAGCAGGCTTTTTTTGCttttt                                                                                                     >  1:231754/1‑150 (MQ=255)
                                                    ggATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGATTTTTTCGTTTAAAAGACTGGATTTTTGCCAAAACAGGGTTTTGTGACTTTTTGt                                                                                                   >  2:1086563/1‑149 (MQ=255)
                                                               caACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTTTTTAACAGCTTGGGTTTTCGCAAAAGCAGGCGTTTTTTACTTTTTGGGCtttt                                                                                             >  2:148649/1‑144 (MQ=255)
                                                                  ccTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAACAACTTTGATTCTCGCCAGAGCAGGGTTTTTTCACTTTTTGTGCTTTttct                                                                                          >  1:439435/1‑144 (MQ=255)
                                                                    tCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAAAAGCTTTGGTTTTCGCAAGAGCAGGCGTTTTTTTCTTTTTGTGCTTTTTTTGCGCTgg                                                                                   >  2:437613/1‑149 (MQ=255)
                                                                     cGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAACAGCTTTGGATTTCCCAAAAGGAGGCTTTTTTTCCTTTTTGGGCtttt                                                                                             >  2:378091/1‑138 (MQ=255)
                                                                     cGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAACAGCTTGGGTTTTCGCAAGAGCCGGCTTTTGTTACTTTTTGTGCtttt                                                                                             >  2:87149/1‑138 (MQ=255)
                                                                     cGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGt                                                                                  <  2:211922/149‑1 (MQ=255)
                                                                     cGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTc                                                                                 <  2:916654/150‑1 (MQ=255)
                                                                     cGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCt                                                                                <  2:959783/151‑1 (MQ=255)
                                                                     cGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCt                                                                                <  2:218901/151‑1 (MQ=255)
                                                                               tCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCa                                                                      <  1:1084060/150‑1 (MQ=255)
                                                                                     tCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTCGTTTAAAAGCTTGGATTTTCGCCAGAGCAGGCTTTTTTTTCTTTTTGTGCTTTTTCTCCGCTTGGCTATTGCCTCCttt                                                                   >  2:246124/1‑147 (MQ=255)
                                                                                        tttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTcc                                                                 <  2:774111/146‑1 (MQ=255)
                                                                                        tttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCa                                                                <  2:66124/147‑1 (MQ=255)
                                                                                              tCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAAAAGCTTGGGTTCTTGCCAGAGCAGGCTTTTTTTACTTTTTGTGCTTTTTCTTCGCTGGGCTTTTGCCCCCttttcctttgggggtt                                                      >  2:1162029/1‑138 (MQ=255)
                                                                                               cGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTa                                                      <  1:94155/150‑1 (MQ=255)
                                                                                               cGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTGTTTAAAAGCCTGGGTTCTCGCAAGAGCAGGGTTTTTTCACTTTTTGTGCTTTTTTTGCGCCGGGCTATTGCCTCAttttcctt                                                              >  1:59574/1‑137 (MQ=255)
                                                                                                   ttATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAAAAACTTGGATTTTCGCAAGAGGAGGCTTTTTTGCCTTTTTGTGGCTTTTTTTCGCTGGTCTATTGCCCCCTTTTCATTtgcggggttgt                                                   >  1:1067899/1‑143 (MQ=255)
                                                                                                   ttATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTCGGTTAAAAGCTTGGGTTTTCCCCAGAGCAGGGCTTTTTTACTTTTTGGGCTTTTTCTGCGCTGGGCTATTGGCTTATTTTCCATTggggggat                                                      >  1:1147863/1‑144 (MQ=255)
                                                                                                    tATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTTTTTGTTTAAAAACCTTGGTTTTCGCCAAAGGCGGGTTTTTTTACTTTTTGGTTTTTTTTTGTGCTGGGCTTTGGGCTCATTTTCATTtgcgggctt                                                     >  2:155895/1‑142 (MQ=255)
                                                                                                     aTGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCgg                                                <  1:431524/150‑1 (MQ=255)
                                                                                                     aTGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGc                                               <  2:900047/151‑1 (MQ=255)
                                                                                                      tGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGGTTTTTTTGTTTAAAAGCTTTGGTTTTCGCCAGAGCCGGGTTTTTTTACTTTTTGTGCTTTTTTTTTGCTGGGCTATTGGCTCATTTTCATTTGGGGCTATGTTGGCg                                              >  1:612403/1‑151 (MQ=255)
                                                                                                                  ggCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGc                                               <  1:1036987/138‑1 (MQ=255)
                                                                                                                  ggCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTgg                                  <  2:374463/151‑1 (MQ=255)
                                                                                                                      actacGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGa                               <  2:591080/150‑1 (MQ=255)
                                                                                                                           ggCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAAc                         <  2:375125/151‑1 (MQ=255)
                                                                                                                              gTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCa                      <  1:848044/151‑1 (MQ=255)
                                                                                                                                            ggCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTATTgg         <  1:1027225/150‑1 (MQ=255)
                                                                                                                                              cGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGc      <  1:597941/151‑1 (MQ=255)
                                                                                                                                               ggTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCa     <  1:1086896/151‑1 (MQ=255)
                                                                                                                                               agTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGc      >  2:221496/2‑150 (MQ=255)
                                                                                                                                                gTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCAc    >  1:879535/1‑151 (MQ=255)
                                                                                                                                                   ttcttcGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACtt  >  2:885442/1‑150 (MQ=255)
                                                                                                                                                    tcttcGTTAATCAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCTCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACtt  <  2:395086/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
TTTCAGGATGGTATTGGAAGGTTAATAAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGTTTTCTTCGTTAAACAGCTTGGATTCTCGCAAGAGCAGGCTTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGCGACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTT  >  NC_000913/3640832‑3641128

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