Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913535,9010TG40.7% 104.6 / 11.9 62G64G (GGT→GGG) hyihydroxypyruvate isomerase
Reads supporting (aligned to +/- strand):  ref base T (7/28);  new base G (24/0);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 1.23e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCG‑ATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTC  >  NC_000913/535755‑536048
                                                                                                                                                   |                                                                                                                                                      
gATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACt                                                                                                                                                    <  1:725457/151‑1 (MQ=255)
 aTTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTg                                                                                                                                                   <  1:514794/151‑1 (MQ=255)
   tttCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGACTggg                                                                                                                                                 >  1:383132/1‑151 (MQ=255)
            cGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GCGGATGGGAgggggggg                                                                                                                                        >  2:54467/1‑151 (MQ=255)
             gTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGa                                                                                                                                       <  1:545201/151‑1 (MQ=255)
             gTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGa                                                                                                                                       <  1:479080/151‑1 (MQ=255)
              ttttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGGCTGGGCGGCggggggg                                                                                                                                      >  1:304161/1‑149 (MQ=255)
              ttttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACAGGCGG‑CTGaacaggacgg                                                                                                                                             >  2:498778/1‑135 (MQ=255)
                  gagaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGcgc‑g                                                                                                                                 <  2:539127/151‑1 (MQ=255)
                              cAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGt                                                                                                                      <  1:657451/150‑1 (MQ=255)
                              cAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGGCTGGGCGGCgggggggggggggatt                                                                                                                             >  1:143426/1‑133 (MQ=255)
                                  gtgGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTc                                                                                                                 <  2:714953/150‑1 (MQ=255)
                                        ttCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGGATGGGGGGCGGGGGAGGGC‑GGCGTTGCCTTTA‑TTTCtcggc                                                                                                           >  2:447993/1‑146 (MQ=255)
                                            cggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGGCTGGGGGGGGGGGGAGCGC‑GGGATTGCCTGTATTTCCTGGc                                                                                                            >  2:563535/1‑146 (MQ=255)
                                            cggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GCGTATGGGCCGCGGGGGCGGGC‑GGGATTGCCTCTG‑TTCTTGGCCGTgc                                                                                                       >  1:178207/1‑149 (MQ=255)
                                               cgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGCCTGGGCGGCGGGGGGACGC‑GGGAGTGTCTGTA‑tttcctgccggggaaa                                                                                                    >  1:776026/1‑136 (MQ=255)
                                               cgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGACTGGGGGGCGGGGGGGCGC‑GGTATTTCCTGTA‑TTTCTGGcccggagga                                                                                                    >  2:107647/1‑143 (MQ=255)
                                                gTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGGGTGGGCGGCGGGGGGGCGC‑GGGATGGGCTGGA‑TTTCTGGCCggggag                                                                                                     >  1:30009/1‑143 (MQ=255)
                                                 ttGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑TGGAGTGGGGGGCGGGGGGGGGC‑GGGATTGCCTGGT‑TTTCTgggggggga                                                                                                      >  2:478777/1‑139 (MQ=255)
                                                  tGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAagaaga                                                                                                 <  1:792991/150‑1 (MQ=255)
                                                  tGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGGCAGGGGGGGGGGGGGGGGC‑GGGATTTCCTGTA‑TTTCttg                                                                                                             >  2:237220/1‑136 (MQ=255)
                                                       ttATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTc                                                                                            <  1:134664/150‑1 (MQ=255)
                                                              tCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGTCAGGGGGGCGGGGGGGGGG‑GGGGTTGCTTGCA‑TTTCTAGACGTGGTGAAAAGGTTAggg                                                                                         >  1:70345/1‑146 (MQ=255)
                                                                  tATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGt                                                                                  <  2:624576/149‑1 (MQ=255)
                                                                    tGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGACTGGGCGGGGGGGGGGCGC‑GGGATTGCCTGTA‑TTTCTGGCCGTGGAGAAGAGGTTTGGGGTGGCGGag                                                                                >  1:377033/1‑149 (MQ=255)
                                                                     gACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGCGGGAC‑GGGGGGAGGGGGCGCGCGGGGATCGTCCG‑A‑TTCCTGGCAGCGAAGGAAAGTTTGGGGag                                                                                       >  2:605853/1‑140 (MQ=255)
                                                                         aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTagcagcagc                                                                         <  1:609456/151‑1 (MQ=255)
                                                                         aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTagcagcag                                                                          <  1:582870/150‑1 (MQ=255)
                                                                            aCATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCG‑At                                                                     <  1:523714/151‑1 (MQ=255)
                                                                            aCATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGACTGGGGGGCGGGGGGGCGC‑GGGATTGCCTGTA‑TTTCTGGCCGGGGAGAAGAGGTTTGGGGTGGCGTAGCAGCAGCG‑At                                                                     >  1:329671/1‑151 (MQ=255)
                                                                              aTTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGGATGGGGGGCGGGGGGGGGC‑GGGATTTCCTTTT‑TTTCTGGGCGGGGAGAAGAGGTTTGGGGGGGGGGAGGAGGAGCG‑aat                                                                    >  1:90186/1‑148 (MQ=255)
                                                                               ttGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGAATGGGCGGGGGGGGGGGGC‑GGGAATGCCCGGA‑TTTCTGGGCGGGGAGGAGAGGTTCGGGGTGGGGGTGCAGCAGCG‑Att                                                                    >  1:347441/1‑149 (MQ=255)
                                                                                tGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGGATGGGGGGGGGGGGGGAGC‑GGGATTGCCTGTA‑TTTCTGGCCGTGGAGGAGAGTTTCGGGGTGGGGGagcaggcggggatt                                                                    >  2:683800/1‑140 (MQ=255)
                                                                                  aagaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGAGGTGGCGGCGGGGGAGGGG‑GGGATTGCCTGTA‑TTTCCGGACGTGAAGGAAAGTTTATGGGGAGCGTagcagcagcattttgtt                                                                 >  1:224509/1‑142 (MQ=255)
                                                                                     aaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGGGGGC‑GGGATTGCCTGTATTTCTTGGCCG‑GGAGAAAAGGTTTGGGGTGGGGGAGCAGCAGCG‑ATTTGTTAt                                                              >  1:652751/1‑149 (MQ=255)
                                                                                         aaaaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGAATGGGGGGCGGGGGAGGGC‑GGGGTTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGGTTTGGGATGGGGGAGCAGCAGGG‑GTTCGTTATGCCCg                                                         >  1:644763/1‑150 (MQ=255)
                                                                                          aaaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCG‑ATTCGTTATGCCCGTg                                                       <  2:202452/151‑1 (MQ=255)
                                                                                            aCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGACTGGGGGGGGGGGGAGCGC‑GGGATTGCCTGTT‑TTCCCCGGCCTGAAGAAGGGGTTTGGGGGGGGGTAGGAACAGCGAATTCG‑TATGCCCgag                                                       >  2:566266/1‑147 (MQ=255)
                                                                                             cATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGACTGGGGGGCGGGGGGGGGG‑GGGGATTCCTGGA‑TTTCTGGCCCGGGAAAAGGGTTTTGGGGTGGGGTAGCAGCAGCG‑Atttggtttggcc                                                          >  2:152758/1‑135 (MQ=255)
                                                                                                 tgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGGGCGC‑GGGATTGCCTGTA‑TTTCTGGCCGTGAAGAAGAGATTTGGGGTGGGGGAGCAGCAGCG‑ATTCGTTTTGCCCGGGCGCTggg                                                >  1:260714/1‑151 (MQ=255)
                                                                                                     ggCGAGTAATAAACTGGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGTCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCG‑ATTCGTTATTCCCGTGCGCTGGGtaa                                             <  2:759483/150‑1 (MQ=255)
                                                                                                      gCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCG‑ATTCGTTATGCCCGTGCGCTGGGtaata                                           <  2:39643/151‑1 (MQ=255)
                                                                                                        gAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGGGCGC‑GGGATTGCCTGTA‑TTCCTGGCCGTGGAGAAGAGTTTCGGGGTGGCGTAGCAGCAGCG‑ATTCGTTTTGCCCCTGCGCTGGGgaataa                                          >  1:32757/1‑150 (MQ=255)
                                                                                                        gAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGACTGGGCGGCGGGGGGGCGC‑GGGATTGCCTGTA‑TTCCTGGGCGGGGAGAAAAGGTTCGGGGTGGGGGAGCAGCAGGG‑ATTCGTTATTCCCGTGCGGCGGGtaa                                             >  2:734398/1‑147 (MQ=255)
                                                                                                         aGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GGGAGTGGGAGGGGGGGGAGAAG‑GGGATTGAGTGTA‑TTTATGGGCGAGAGAAAGAGGTTTGGGGGGGCGTTGGAGGAGTG‑ATTCGATATGTCCGTggggaggg                                                >  2:259585/1‑136 (MQ=255)
                                                                                                          gTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCG‑ATTCGTTATGCCCGTGCGCTGGGTAATaaaa                                        <  2:495736/150‑1 (MQ=255)
                                                                                                            aataaACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGTAGCAGCG‑ATTCGTTATGCCCGTGCGCTGGGTAATaaaaaaa                                     <  2:792880/151‑1 (MQ=255)
                                                                                                             ataaACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGGGC‑GGGATTTCCTTTA‑TTCCTGGCCGGGAAGAAGAGTTTCGGGGTGGGGTAGGAGCAGGG‑GTTCGTTATGCCCGTGCGGTGGGTAATaaaaaaa                                     >  2:99261/1‑150 (MQ=255)
                                                                                                                aaCTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCG‑ATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAAtt                                   <  1:126503/149‑1 (MQ=255)
                                                                                                                    cGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCG‑ATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACt                               <  1:683300/149‑1 (MQ=255)
                                                                                                                    cGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCG‑ATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGt                             <  2:267838/151‑1 (MQ=255)
                                                                                                                      aaCATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCG‑ATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCt                           <  1:787431/151‑1 (MQ=255)
                                                                                                                         aTACGCTGCACAATTTACCGGCGG‑GGGACTGGGGGGCGGGGGAGGGC‑GGGATTGCCTGTA‑TTTCTGGGCGTGGAGAAGAGGTTCGGGGGGGGGGAGCAGCAGCG‑ATTCGTTTTGCCCCGGCGCTGGGGAATAAAAAAAATAAATGTCTgg                         >  2:249172/1‑150 (MQ=255)
                                                                                                                                gCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCG‑ATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTa                   <  2:304121/149‑1 (MQ=255)
                                                                                                                                  acaATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCG‑ATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAAc               <  1:319029/151‑1 (MQ=255)
                                                                                                                                          ccGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCG‑ATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTg       <  1:514454/151‑1 (MQ=255)
                                                                                                                                          ccGGCGG‑GCGACTGTGCGGCGGGCGCACGC‑GGCATTGCCTGTA‑TTCCTGGCCGTGAAGCAGAGTTTCCGGACGGCGGAGCAGCAGCG‑CTTCGTTATGTCCGCGCACTCGGTATTAAAAAAATTTACTGTCTGGCCCGGAAAACGCCGGGTg       >  2:456881/1‑151 (MQ=255)
                                                                                                                                           cggcgg‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCG‑ATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTgg      <  2:96320/151‑1 (MQ=255)
                                                                                                                                               gg‑gTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCG‑ATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTc  <  2:350094/151‑1 (MQ=255)
                                                                                                                                               gg‑gTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCG‑ATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTc  >  2:720088/1‑151 (MQ=255)
                                                                                                                                               gg‑gTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCG‑ATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTc  <  2:122536/151‑1 (MQ=255)
                                                                                                                                                  gTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCG‑ATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTc  >  1:67005/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGG‑GTGACTGGGCGGCGGGGGAGCGC‑GGTATTGCCTGTA‑TTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCG‑ATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTC  >  NC_000913/535755‑536048

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