Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913436,5040AC41.4% 97.6 / 17.2 60T291P (ACC→CCC) thiLthiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base A (25/9);  new base C (0/24);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 3.10e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTG  >  NC_000913/436359‑436650
                                                                                                                                                  |                                                                                                                                                    
tGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGt                                                                                                                                                >  2:119682/1‑151 (MQ=255)
 gATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGt                                                                                                                                                <  1:937498/150‑1 (MQ=255)
    gcgcTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATc                                                                                                                                             <  2:905586/150‑1 (MQ=255)
     cgcTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑gg                                                                                                                                          <  1:904026/151‑1 (MQ=255)
         ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGc                                                                                                                                        <  2:736376/149‑1 (MQ=255)
           ctcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCaaa                                                                                                                                     >  1:38546/1‑150 (MQ=255)
             cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATg                                                                                                                                   <  2:360976/150‑1 (MQ=255)
             cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATg                                                                                                                                   <  2:389676/150‑1 (MQ=255)
             cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGa                                                                                                                                  <  1:993324/151‑1 (MQ=255)
                cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGAccg                                                                                                                               >  1:585847/1‑151 (MQ=255)
                             aaCAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAgg                                                                                                                  >  2:385822/1‑151 (MQ=255)
                                aGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGc                                                                                                                >  1:522194/1‑150 (MQ=255)
                                  gcgctgcgctGGGGTCTCTCGGGCGGTGAATTTTACGGTTTGTGTTTCCCTGGGCCGAAACTAACCCGTGGCGCGATGGGTGTGGCTCCGGCACACGGGGG‑GGTCCTTTTTCCC‑TGTAGC‑GGAAAAATGACCGCCGATATCGAAGGGCttt                                                                                                             <  1:931814/151‑1 (MQ=255)
                                      tgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGttt                                                                                                         >  2:723746/1‑151 (MQ=255)
                                      tgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGttt                                                                                                         >  2:557610/1‑151 (MQ=255)
                                       gcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGtttt                                                                                                        >  1:993561/1‑151 (MQ=255)
                                                 ctctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                              >  2:1125419/1‑151 (MQ=255)
                                                     ctgggggaaaaGATTAGGGGTTTTTTTTCGCTGCCCCGGAACTGACCGGGGGCGCCGGGGGTGGGGC‑CCCCGACCGCTGGGCCGTCCTTTTTCCC‑TGTTTG‑GGGAAAATGTCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                              <  1:704262/138‑1 (MQ=255)
                                                       gggggTGAAGTTTCCGAGTTGTGTTTCACTGTGCCGGAACGGACCGGTGGCGCGCTGGATGTGGCTCTCGGACCCTGGGG‑CGTCCCTTTTCCC‑TGTTCG‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGaa                                                                                         <  2:88668/147‑1 (MQ=255)
                                                       gggggTGAAGATTAGGGGTTGTTTTTCCCTGTGCCGGAACTAACCCGTGGCGCGCGGGAGGTGCCTCTCGGACACTGGGG‑CGTACCTTTTCCC‑TTTACC‑GGGAAAAGGTCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGaa                                                                                         <  1:418993/147‑1 (MQ=255)
                                                       gggggTAAAGATTGCGTTTTTTTTTGCACTGCGCGGGAACTGACCCGGGGCGCGTTGGGGGGGCCTCCCGGACCCGGGGG‑CTTACCTTTTCCC‑TTTACC‑GGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                           <  1:1055157/145‑1 (MQ=255)
                                                       agCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAc                                                                                        >  1:136883/2‑151 (MQ=255)
                                                        ggggTGAAGTTTACGAGTTGTGTTTCCCTGGGCCGAAACTGACCCGGGGCGCGCTGGATGTGGCTCCCGGACCCCGGGG‑CTTCCCTTTTCCT‑TGTACG‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAAcc                                                                                       <  1:777153/149‑1 (MQ=255)
                                                        ggggTGAAGATTACGGTTTGTGTTTCCCTGTGCCGGAACTGAACCGTGGCGCGGTGGGTGGGGCTCTCGGACACGGGGG‑GGTCCCTTTTCCC‑TGTATC‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGAc                                                                                               <  2:151058/141‑1 (MQ=255)
                                                          gggaaagttaaCGATTTGTTTTTCACTGCGCGGAACCGGACCCGGGGCGCGGTGGATGTGGCTCCCGAACCCGGGGG‑CGTCCTTTTTCCT‑TGTATC‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTg                                                                                     <  2:465208/141‑1 (MQ=255)
                                                            tGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGtt                                                                                   >  2:677171/1‑151 (MQ=255)
                                                              aagttgaggggttgttTTTCCCTGGGCCGGACCTGACCGGGGGCGCGTTGGATGTGGCTCCCGGACCCCTGGG‑CGTCCCGTTTCCT‑TGTTTG‑GGAAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTa                                                                                  <  2:963728/140‑1 (MQ=255)
                                                                   tccgtgttttttttCACTGTGCCGAACCTGACCCGTGGCGCGTTGGTTGTGCCCCCCGGACCCGGGGG‑CTTCCTTTTTCCC‑TGTATG‑GGCAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTa                                                                             <  1:637546/139‑1 (MQ=255)
                                                                     cGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGAc                                                                          >  1:147737/1‑151 (MQ=255)
                                                                          tgtgTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                                      >  1:72690/1‑150 (MQ=255)
                                                                           gtgtTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaa                                                                     >  2:988751/1‑150 (MQ=255)
                                                                              tttCCCTGAGCCGGAATGGATCCGCGGCGCGTTGGTTGTGCCTCTGGGACCCTGGGG‑CGTCCTTTTTCCT‑TGTACC‑GGGCAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                  <  1:974211/150‑1 (MQ=255)
                                                                              tttCACTGTGCCGGAACGGACCGGGGGCGCGTGGGATGGGCCTCCGGGCCCCGTGGG‑CGTCCCTTTTCCC‑TGTAGC‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaa                                                                     <  1:119682/147‑1 (MQ=255)
                                                                                  aCTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATg                                                             >  2:172411/1‑151 (MQ=255)
                                                                                  aCTGGGCCGGAACGGACCCGTGGCGCGCTGGTTGTGGCTCTCGGACACCGGGG‑GGTCCCGTTTCCT‑TTTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATg                                                             <  2:55678/151‑1 (MQ=255)
                                                                                     gggCCGGAACTGACCGGGGGCGCGTTGGATGTGTCTCTCGGCCCCTTGGG‑CGTCCTGTTTCCT‑TGTATC‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc                                                           <  2:470269/148‑1 (MQ=255)
                                                                                        ccGAAACTGACCGGGGGCGCGTTGGTTGTGCCTCTCGGACCCGTGGG‑CGTGCCTTTTCCT‑TGTTTG‑GGGAAAAAGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGa                                                            <  1:115565/146‑1 (MQ=255)
                                                                                        ccGAAACTGACCCGGGGCGCGCTGGATGTGGCCCCCGGACCCTGGGG‑CGTACCTTTTACC‑TGTAGC‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATg                                                             <  1:1128558/145‑1 (MQ=255)
                                                                                          ggAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGCAAGGATATGACCAtttt                                                     >  2:421005/1‑151 (MQ=255)
                                                                                          gaaacgaaaCGGTGGCGGGCTGGATGGGGCCCCCGCCCACGGGGG‑CGTCCTTTTTCCT‑TTTATG‑GGCAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtt                                                       <  1:796478/142‑1 (MQ=255)
                                                                                            aaCTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                   >  1:969536/1‑151 (MQ=255)
                                                                                               tGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                >  1:489578/1‑151 (MQ=255)
                                                                                               ggaccggggccgcgcgGGAGGTGCCTCTCGAAC‑CCGGGG‑GCTTCCGTTTTCCTTGTATG‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                <  2:176743/141‑1 (MQ=255)
                                                                                                 accgggggCGCGTTGGATGTGGCTCTCGGACCCGTGGG‑CGTACCGTTTCCT‑TGTATC‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGc                                              <  2:711905/145‑1 (MQ=255)
                                                                                                   ccGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGcc                                             >  1:590044/1‑150 (MQ=255)
                                                                                                    cGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCAt                                           >  1:386494/1‑151 (MQ=255)
                                                                                                     ggggggcgggggttgtgTCTCTGGGACGCTGGGG‑CGTCCCGTTTACT‑TGGAGG‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCAt                                           <  2:1195049/137‑1 (MQ=255)
                                                                                                              tGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcg                                 >  2:273253/1‑151 (MQ=255)
                                                                                                               ggttgtgtcccggagACATCTGGG‑AGTCCTGTTT‑CC‑TGTTTTGGGGAAAATTACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcg                                 <  1:838338/136‑1 (MQ=255)
                                                                                                               ggttgggcctctCGGCCCCCGGGG‑CGTCCCTTTTCCT‑TGTAGG‑GGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATg                                    <  2:1005465/139‑1 (MQ=255)
                                                                                                                      gCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGt                         >  1:472820/1‑151 (MQ=255)
                                                                                                                          tcGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGcc                     >  2:370127/1‑151 (MQ=255)
                                                                                                                            ggaccgggggg‑CGTCCTTTTTCCC‑TTTACG‑GGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTc                        <  1:1160889/138‑1 (MQ=255)
                                                                                                                            gacccccgggg‑CGTACTTTTTCCT‑TGTATG‑GGAAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTg                   <  2:326044/143‑1 (MQ=255)
                                                                                                                               cactgggg‑CGTCCCTTTTCCC‑TGTATG‑GGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCt                    <  1:651639/142‑1 (MQ=255)
                                                                                                                                   tgggggttcgcGTTTCCT‑TGTATG‑GGGCAAGGACCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGa               <  2:311245/139‑1 (MQ=255)
                                                                                                                                    ggg‑CGTCCCTTTTCCC‑TGTAGG‑GGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg             <  2:217607/149‑1 (MQ=255)
                                                                                                                                     gg‑CGTACTTTTTCCT‑TGTATC‑GGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGt          <  2:38546/151‑1 (MQ=255)
                                                                                                                                         gTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAAt       >  1:836559/1‑151 (MQ=255)
                                                                                                                                               tttACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg  >  2:88279/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGATGCGCTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGG‑CGTACCGTTTACC‑TGTATC‑GGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTG  >  NC_000913/436359‑436650

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