Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913436,5160CA44.3% 64.3 / 31.8 61Q295K (CAA→AAA) thiLthiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (25/9);  new base A (0/27);  total (25/36)
Fisher's exact test for biased strand distribution p-value = 7.09e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

TTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTG  >  NC_000913/436368‑436661
                                                                                                                                                      |                                                                                                                                                  
ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGc                                                                                                                                                    <  2:853305/149‑1 (MQ=255)
ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGc                                                                                                                                                    <  1:771345/149‑1 (MQ=255)
    cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGa                                                                                                                                             <  1:15198/151‑1 (MQ=255)
       cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGCAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGCGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑Tg                                                                                                                                              >  2:501521/1‑147 (MQ=255)
       cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGAccg                                                                                                                                          <  1:313348/151‑1 (MQ=255)
        aTGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGAccgc                                                                                                                                         <  1:440476/151‑1 (MQ=255)
             gAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGat                                                                                                                                     >  1:183680/1‑150 (MQ=255)
                 cGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCg                                                                                                                                 >  1:467745/1‑150 (MQ=255)
                    aaCAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAgg                                                                                                                             >  2:397539/1‑151 (MQ=255)
                    aaCAGGCGCTGCCCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGAGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAgg                                                                                                                             >  2:756710/1‑151 (MQ=255)
                       aGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCt                                                                                                                          >  1:692400/1‑151 (MQ=255)
                            ctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGtt                                                                                                                     >  2:625405/1‑151 (MQ=255)
                            ctgcgcggGGCGCTCTCGGGGGGTGACGTTTGCGGGTTGTGTTTCACTGTGCCGGAACTGCCCGGTGTCGCGTTGGATTTGGTTCGCGGCCAGGTGGGAGCCCCGTTTACC‑GGTATC‑GGCAAAC‑TTCCCGCCGATATCGAAGGGCttt                                                                                                                        <  2:242381/148‑1 (MQ=255)
                             tgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGttt                                                                                                                    >  2:120581/1‑151 (MQ=255)
                                gctgGGCCCTCTCTGGGGGGAAAGATTAGGAGTTGTTTTTCACTGTGCCGAAACTGACCCGTGGCGCGCTGGATGGGGCTCGCGGCCCCCGGGGCGTCCCGTTTACC‑TGTAGC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTg                                                                                                                       <  2:348775/145‑1 (MQ=255)
                                   gggcgtgttttGGGGGTGAAGATTAGGGGTTGGGTTTCACTGTGCGGGAACTGACCGGTGGGGCGTGGGGTGTGCCTCCCGGACACTTGGGGGTCCTGTTTACC‑TGTAGG‑GGGAAAG‑TGACCGCCGATATCGAAGGGCTTTg                                                                                                                       <  2:646858/132‑1 (MQ=255)
                                      cgcTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg                                                                                                            >  1:4843/1‑150 (MQ=255)
                                        ctctcgGGCGGTGAAGATTAGGGGTTGTTTTCCCCTGTCCGGGAACGGATGCGGGGCGCGCGGGAGGTGGCTCGCGCACACCTGGGCGGACTTTTTTCT‑TGTCGT‑GGAAAAG‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGa                                                                                                           <  2:327133/149‑1 (MQ=255)
                                             tGGAGGTGAAGCTTACGAGGTGGGTTTAACTGCGCCGGAACTGCACCGCGGCGCGCTGGATGTGGGTTTTGGACACCTGGGCGTACCGTTTAAC‑TGTATC‑GAGCAAA‑TGACCGCCGATATTGAAGGGCTTTGCTTTATTTGTGACGGCGa                                                                                                     >  2:563543/1‑150 (MQ=255)
                                             ggggggTGAAGTTTAGGGTTTGTGTTTCCTTGGGCCGCAACTGACCGGTGGCGCGCTGGTTGTGGCTCTCGGCCCCGTGGGCGTACTTTTTTCT‑TGTAGG‑GAGAAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg                                                                                                            <  2:638039/139‑1 (MQ=255)
                                             ggggggTAAAGATTACGATTTGTTTTCCGCTGTGCCGGAACTGACCGGGGGGCCGCTGGATGTGGCTCCCGGGCTTTTGGGC‑CATCGTTTACCTTGTATG‑GGGAAAA‑GGGCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg                                                                                                        <  1:613361/143‑1 (MQ=255)
                                              gggggTGAAGTTTACGGGTTGTGTTTCACTGTGCCGGAACTGAACCGGGGCGCGCTGGTTGTGGCTCTCGGACCCCGGGGCGTCCCTTTTCCC‑TGTTTG‑GGGAAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                      <  1:464764/145‑1 (MQ=255)
                                              gggggTGAAGATTACGATTTGTGTTTCACTGTGCCGGAACGGAACCGTGGCGCGCTGGTTGTGGCTCTCGGACACGTGGGCTTCCCGTTTCCC‑TGTATC‑GGGAAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGaa                                                                                                    <  1:824260/147‑1 (MQ=255)
                                               ggggtagaGATTACGAGTTGGGTGTCAAGGCGGTGGAACGAACCCGGGGCGCGCTGATTGTGGCTGTGGAGCGCCGGGGGGTCCGGTTTCAC‑TGCACA‑GGACAAA‑CCGCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                      <  1:86448/140‑1 (MQ=255)
                                                     aaGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                            >  1:876443/1‑151 (MQ=255)
                                                              aGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                                   >  1:546476/1‑151 (MQ=255)
                                                                ttGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg                                                                                  >  2:564178/1‑150 (MQ=255)
                                                                   tgtTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaa                                                                                >  2:535146/1‑149 (MQ=255)
                                                                    tttgctgtgtgCCGAACATGACCCGGGGCGGGTGGGATGTTGCCCTCGGACCTCTGGGCGT‑CCGTTTTTG‑TTTATCGGGGAAAA‑GGTCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAg                                                                              <  2:648111/143‑1 (MQ=255)
                                                                     tttCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                             >  1:47710/1‑150 (MQ=255)
                                                                     tttCACGGTGCCGGACCTGAACGGGGGCGCGCTGGATGTGGCCCTCGGCCACCTGGGCGTTCCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAg                                                                              <  2:328470/149‑1 (MQ=255)
                                                                     ttcccttgtgCCGGAAGAGAGCGGGGGCGGGCGGGGTGGGGCTCGGGGACCCGGGGGCGTACCTGTTCCC‑TGTTTC‑GGCCAAA‑GGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                                                 <  1:449816/140‑1 (MQ=255)
                                                                      ttCACTGTGCGGGAACTGACCGGTGGCGCGCTGGTTGGGGCTCTCGGACCCTGGGGCGTCCCTTTTCCC‑TGTATC‑GGGAAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGa                                                                            <  1:501875/150‑1 (MQ=255)
                                                                       tCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGat                                                                           >  1:341967/1‑150 (MQ=255)
                                                                         aCTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATg                                                                        >  1:811843/1‑151 (MQ=255)
                                                                             tgCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc                                                                      >  2:651151/1‑149 (MQ=255)
                                                                              gCCGGAACTGAACCGGGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTACCTTTTCCC‑TGTATC‑GGGAAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAt                                                                   <  1:305429/151‑1 (MQ=255)
                                                                               cgggAACTGAACGGTGGCGCGCTGGTTGTGGCTCTGGGACCCCTGGGCGTCCCTTTTCCT‑TGTATG‑GGGAAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCa                                                                    <  1:933070/147‑1 (MQ=255)
                                                                               ccGGTACTGACCCGTGGCGCGGTGGTGGGCGCTCGCGGTCCCGGGGTTGGACGTTTTACT‑TGTGTC‑GGACAAT‑GGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGata                                                                          <  2:370404/143‑1 (MQ=255)
                                                                                 ggAACTGACCCGGGGCGCGCGGGATGTGGCTCCCGCCCCCTGGGGCGTCCCTTTTCCC‑TTTATC‑GGGAAAA‑TGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt                                                                <  2:883536/151‑1 (MQ=255)
                                                                                            ggggcgcgtgggtggtgGCTCGCGGACCCCGGGGCGTCCCGTTTCCC‑TGTATG‑GGGAAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCTAt                                                      <  1:817613/136‑1 (MQ=255)
                                                                                            gTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCAt                                                      >  2:727507/1‑150 (MQ=255)
                                                                                             gggCGCGTGGGATGTGGCTCTGGGACCCTGGGGCGTCCCTTTTCCC‑TGTCCC‑GGGAAAA‑GGACCGCCGATATCGAAGGGCTTTTTTTTTTTCTTTCCGTCGAACCTGTTACATTATCCTGGGAAGTATATGACCATTTTGCCACGCCCTaa                                                    <  1:791075/150‑1 (MQ=255)
                                                                                              ggCGCGGTGTATGTGGCTCTCGGACCCTGGGGCGTACTTTTTCCC‑TGT‑TG‑GGGAAAAGTGTCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCa                                                       <  2:229943/147‑1 (MQ=255)
                                                                                              ggCGCGGGGGGTGTGTCCCTCGGACCCTGGGGCGTCCTGTTTCCT‑TTTCTG‑GGGAAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa                                                    <  2:366346/150‑1 (MQ=255)
                                                                                                   gttGGTTGTGGCTCTCGGACCCCGGGGCGTCCCTTTTCCC‑TGTATC‑GGGAAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATg                                               <  1:243933/148‑1 (MQ=255)
                                                                                                   gtgggATGTGTCTCTCGGACACTTGTTCGTCCCTTTTCCC‑TGTTTG‑GGGAAAA‑TGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCa                                                       <  2:563122/139‑1 (MQ=255)
                                                                                                               tctcGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg                                  >  1:100169/1‑151 (MQ=255)
                                                                                                                cccGGACCCCGGGGCGTCCCTTTTCCC‑TGTATG‑GGGAAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg                                  <  1:582517/148‑1 (MQ=255)
                                                                                                                   gccccccggggCGTCCCGTTTCCT‑TGTTGG‑GGGAAAA‑TTACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCt                               <  1:860009/142‑1 (MQ=255)
                                                                                                                    gacccggggggGTACCTTTTCCC‑TGTAGG‑GGGAAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTc                                   <  1:765311/138‑1 (MQ=255)
                                                                                                                     acacCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                             >  2:208426/1‑150 (MQ=255)
                                                                                                                      cacCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                             >  1:472972/1‑149 (MQ=255)
                                                                                                                        ccTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                             >  2:714938/1‑147 (MQ=255)
                                                                                                                             ggcgtcctttttccgtGTACC‑GGGAAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGt                     <  1:259387/136‑1 (MQ=255)
                                                                                                                             gCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGt                     >  1:43154/1‑150 (MQ=255)
                                                                                                                                     tttCCC‑TGTACG‑GGCCAAT‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg             <  2:912888/150‑1 (MQ=255)
                                                                                                                                      tttccttGTATG‑GGGAAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg             <  1:59459/144‑1 (MQ=255)
                                                                                                                                      ttACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGc           >  2:535661/1‑151 (MQ=255)
                                                                                                                                        tccttGTATG‑GGCAAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGc           <  1:283786/146‑1 (MQ=255)
                                                                                                                                            tGTATG‑GGGAAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCAt         <  1:32152/148‑1 (MQ=255)
                                                                                                                                                   gggAAAA‑TGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACtt   <  1:502187/148‑1 (MQ=255)
                                                                                                                                                   gggAAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTg  <  2:917733/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACC‑TGTATC‑GGGCAAA‑TGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTG  >  NC_000913/436368‑436661

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