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→GGGhyihydroxypyruvate 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

GATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTC  >  NC_000913/535755‑536048
                                                                                                                                                   |                                                                                                                                                      
gATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACt                                                                                                                                                    <  1:725457/151‑1 (MQ=255)
 aTTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTg                                                                                                                                                   <  1:514794/151‑1 (MQ=255)
   tttCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTggg                                                                                                                                                 >  1:383132/1‑151 (MQ=255)
            cGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGCGGATGGGAgggggggg                                                                                                                                        >  2:54467/1‑151 (MQ=255)
             gTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGa                                                                                                                                       <  1:545201/151‑1 (MQ=255)
             gTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGa                                                                                                                                       <  1:479080/151‑1 (MQ=255)
              ttttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGCGGCggggggg                                                                                                                                      >  1:304161/1‑149 (MQ=255)
              ttttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACAGGCGGCTGaacaggacgg                                                                                                                                             >  2:498778/1‑135 (MQ=255)
                  gagaAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGcgcg                                                                                                                                 <  2:539127/151‑1 (MQ=255)
                              cAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGt                                                                                                                      <  1:657451/150‑1 (MQ=255)
                              cAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGCGGCgggggggggggggatt                                                                                                                             >  1:143426/1‑133 (MQ=255)
                                  gtgGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTc                                                                                                                 <  2:714953/150‑1 (MQ=255)
                                        ttCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGATGGGGGGCGGGGGAGGGCGGCGTTGCCTTTATTTCtcggc                                                                                                           >  2:447993/1‑146 (MQ=255)
                                            cggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGGGGGGGAGCGCGGGATTGCCTGTATTTCCTGGc                                                                                                            >  2:563535/1‑146 (MQ=255)
                                            cggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGCGTATGGGCCGCGGGGGCGGGCGGGATTGCCTCTGTTCTTGGCCGTgc                                                                                                       >  1:178207/1‑149 (MQ=255)
                                               cgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGCCTGGGCGGCGGGGGGACGCGGGAGTGTCTGTAtttcctgccggggaaa                                                                                                    >  1:776026/1‑136 (MQ=255)
                                               cgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCGGTATTTCCTGTATTTCTGGcccggagga                                                                                                    >  2:107647/1‑143 (MQ=255)
                                                gTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGGTGGGCGGCGGGGGGGCGCGGGATGGGCTGGATTTCTGGCCggggag                                                                                                     >  1:30009/1‑143 (MQ=255)
                                                 ttGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGTGGAGTGGGGGGCGGGGGGGGGCGGGATTGCCTGGTTTTCTgggggggga                                                                                                      >  2:478777/1‑139 (MQ=255)
                                                  tGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAagaaga                                                                                                 <  1:792991/150‑1 (MQ=255)
                                                  tGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCAGGGGGGGGGGGGGGGGCGGGATTTCCTGTATTTCttg                                                                                                             >  2:237220/1‑136 (MQ=255)
                                                       ttATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTc                                                                                            <  1:134664/150‑1 (MQ=255)
                                                              tCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGTCAGGGGGGCGGGGGGGGGGGGGGTTGCTTGCATTTCTAGACGTGGTGAAAAGGTTAggg                                                                                         >  1:70345/1‑146 (MQ=255)
                                                                  tATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGt                                                                                  <  2:624576/149‑1 (MQ=255)
                                                                    tGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGGGGGGGGGCGCGGGATTGCCTGTATTTCTGGCCGTGGAGAAGAGGTTTGGGGTGGCGGag                                                                                >  1:377033/1‑149 (MQ=255)
                                                                     gACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGCGGGACGGGGGGAGGGGGCGCGCGGGGATCGTCCGATTCCTGGCAGCGAAGGAAAGTTTGGGGag                                                                                       >  2:605853/1‑140 (MQ=255)
                                                                         aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTagcagcagc                                                                         <  1:609456/151‑1 (MQ=255)
                                                                         aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTagcagcag                                                                          <  1:582870/150‑1 (MQ=255)
                                                                            aCATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGAt                                                                     <  1:523714/151‑1 (MQ=255)
                                                                            aCATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGCCGGGGAGAAGAGGTTTGGGGTGGCGTAGCAGCAGCGAt                                                                     >  1:329671/1‑151 (MQ=255)
                                                                              aTTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGATGGGGGGCGGGGGGGGGCGGGATTTCCTTTTTTTCTGGGCGGGGAGAAGAGGTTTGGGGGGGGGGAGGAGGAGCGaat                                                                    >  1:90186/1‑148 (MQ=255)
                                                                               ttGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGCGGGGGGGGGGGGCGGGAATGCCCGGATTTCTGGGCGGGGAGGAGAGGTTCGGGGTGGGGGTGCAGCAGCGAtt                                                                    >  1:347441/1‑149 (MQ=255)
                                                                                tGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGATGGGGGGGGGGGGGGAGCGGGATTGCCTGTATTTCTGGCCGTGGAGGAGAGTTTCGGGGTGGGGGagcaggcggggatt                                                                    >  2:683800/1‑140 (MQ=255)
                                                                                  aagaaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAGGTGGCGGCGGGGGAGGGGGGGATTGCCTGTATTTCCGGACGTGAAGGAAAGTTTATGGGGAGCGTagcagcagcattttgtt                                                                 >  1:224509/1‑142 (MQ=255)
                                                                                     aaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGGGCGGGATTGCCTGTATTTCTTGGCCGGGAGAAAAGGTTTGGGGTGGGGGAGCAGCAGCGATTTGTTAt                                                              >  1:652751/1‑149 (MQ=255)
                                                                                         aaaaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAATGGGGGGCGGGGGAGGGCGGGGTTGCCTGTATTCCTGGCCGTGAAGAAGAGGTTTGGGATGGGGGAGCAGCAGGGGTTCGTTATGCCCg                                                         >  1:644763/1‑150 (MQ=255)
                                                                                          aaaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTg                                                       <  2:202452/151‑1 (MQ=255)
                                                                                            aCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGGGGGGGAGCGCGGGATTGCCTGTTTTCCCCGGCCTGAAGAAGGGGTTTGGGGGGGGGTAGGAACAGCGAATTCGTATGCCCgag                                                       >  2:566266/1‑147 (MQ=255)
                                                                                             cATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGGGGGGGGATTCCTGGATTTCTGGCCCGGGAAAAGGGTTTTGGGGTGGGGTAGCAGCAGCGAtttggtttggcc                                                          >  2:152758/1‑135 (MQ=255)
                                                                                                 tgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTTCTGGCCGTGAAGAAGAGATTTGGGGTGGGGGAGCAGCAGCGATTCGTTTTGCCCGGGCGCTggg                                                >  1:260714/1‑151 (MQ=255)
                                                                                                     ggCGAGTAATAAACTGGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGTCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATTCCCGTGCGCTGGGtaa                                             <  2:759483/150‑1 (MQ=255)
                                                                                                      gCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGtaata                                           <  2:39643/151‑1 (MQ=255)
                                                                                                        gAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTCCTGGCCGTGGAGAAGAGTTTCGGGGTGGCGTAGCAGCAGCGATTCGTTTTGCCCCTGCGCTGGGgaataa                                          >  1:32757/1‑150 (MQ=255)
                                                                                                        gAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGCGCGGGATTGCCTGTATTCCTGGGCGGGGAGAAAAGGTTCGGGGTGGGGGAGCAGCAGGGATTCGTTATTCCCGTGCGGCGGGtaa                                             >  2:734398/1‑147 (MQ=255)
                                                                                                         aGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGAGTGGGAGGGGGGGGAGAAGGGGATTGAGTGTATTTATGGGCGAGAGAAAGAGGTTTGGGGGGGCGTTGGAGGAGTGATTCGATATGTCCGTggggaggg                                                >  2:259585/1‑136 (MQ=255)
                                                                                                          gTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATaaaa                                        <  2:495736/150‑1 (MQ=255)
                                                                                                            aataaACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGTAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATaaaaaaa                                     <  2:792880/151‑1 (MQ=255)
                                                                                                             ataaACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGGGCGGGATTTCCTTTATTCCTGGCCGGGAAGAAGAGTTTCGGGGTGGGGTAGGAGCAGGGGTTCGTTATGCCCGTGCGGTGGGTAATaaaaaaa                                     >  2:99261/1‑150 (MQ=255)
                                                                                                                aaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAAtt                                   <  1:126503/149‑1 (MQ=255)
                                                                                                                    cGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACt                               <  1:683300/149‑1 (MQ=255)
                                                                                                                    cGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGt                             <  2:267838/151‑1 (MQ=255)
                                                                                                                      aaCATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCt                           <  1:787431/151‑1 (MQ=255)
                                                                                                                         aTACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGAGGGCGGGATTGCCTGTATTTCTGGGCGTGGAGAAGAGGTTCGGGGGGGGGGAGCAGCAGCGATTCGTTTTGCCCCGGCGCTGGGGAATAAAAAAAATAAATGTCTgg                         >  2:249172/1‑150 (MQ=255)
                                                                                                                                gCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTa                   <  2:304121/149‑1 (MQ=255)
                                                                                                                                  acaATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAAc               <  1:319029/151‑1 (MQ=255)
                                                                                                                                          ccGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTg       <  1:514454/151‑1 (MQ=255)
                                                                                                                                          ccGGCGGGCGACTGTGCGGCGGGCGCACGCGGCATTGCCTGTATTCCTGGCCGTGAAGCAGAGTTTCCGGACGGCGGAGCAGCAGCGCTTCGTTATGTCCGCGCACTCGGTATTAAAAAAATTTACTGTCTGGCCCGGAAAACGCCGGGTg       >  2:456881/1‑151 (MQ=255)
                                                                                                                                           cggcggGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTgg      <  2:96320/151‑1 (MQ=255)
                                                                                                                                               gggTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTc  <  2:350094/151‑1 (MQ=255)
                                                                                                                                               gggTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTc  >  2:720088/1‑151 (MQ=255)
                                                                                                                                               gggTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTc  <  2:122536/151‑1 (MQ=255)
                                                                                                                                                  gTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTc  >  1:67005/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GATTTTCTCGCCCGTTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGCGGTATTGCCTGTATTCCTGGCCGTGAAGAAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGCGCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTC  >  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: