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

TTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTG  >  NC_000913/436368‑436661
                                                                                                                                                      |                                                                                                                                                  
ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGc                                                                                                                                                    <  2:853305/149‑1 (MQ=255)
ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGc                                                                                                                                                    <  1:771345/149‑1 (MQ=255)
    cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGa                                                                                                                                             <  1:15198/151‑1 (MQ=255)
       cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGCAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGCGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATg                                                                                                                                              >  2:501521/1‑147 (MQ=255)
       cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccg                                                                                                                                          <  1:313348/151‑1 (MQ=255)
        aTGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccgc                                                                                                                                         <  1:440476/151‑1 (MQ=255)
             gAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGat                                                                                                                                     >  1:183680/1‑150 (MQ=255)
                 cGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCg                                                                                                                                 >  1:467745/1‑150 (MQ=255)
                    aaCAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAgg                                                                                                                             >  2:397539/1‑151 (MQ=255)
                    aaCAGGCGCTGCCCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGAGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAgg                                                                                                                             >  2:756710/1‑151 (MQ=255)
                       aGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCt                                                                                                                          >  1:692400/1‑151 (MQ=255)
                            ctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGtt                                                                                                                     >  2:625405/1‑151 (MQ=255)
                            ctgcgcggGGCGCTCTCGGGGGGTGACGTTTGCGGGTTGTGTTTCACTGTGCCGGAACTGCCCGGTGTCGCGTTGGATTTGGTTCGCGGCCAGGTGGGAGCCCCGTTTACCGGTATCGGCAAACTTCCCGCCGATATCGAAGGGCttt                                                                                                                        <  2:242381/148‑1 (MQ=255)
                             tgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGttt                                                                                                                    >  2:120581/1‑151 (MQ=255)
                                gctgGGCCCTCTCTGGGGGGAAAGATTAGGAGTTGTTTTTCACTGTGCCGAAACTGACCCGTGGCGCGCTGGATGGGGCTCGCGGCCCCCGGGGCGTCCCGTTTACCTGTAGCGGGCAAATGACCGCCGATATCGAAGGGCTTTg                                                                                                                       <  2:348775/145‑1 (MQ=255)
                                   gggcgtgttttGGGGGTGAAGATTAGGGGTTGGGTTTCACTGTGCGGGAACTGACCGGTGGGGCGTGGGGTGTGCCTCCCGGACACTTGGGGGTCCTGTTTACCTGTAGGGGGAAAGTGACCGCCGATATCGAAGGGCTTTg                                                                                                                       <  2:646858/132‑1 (MQ=255)
                                      cgcTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg                                                                                                            >  1:4843/1‑150 (MQ=255)
                                        ctctcgGGCGGTGAAGATTAGGGGTTGTTTTCCCCTGTCCGGGAACGGATGCGGGGCGCGCGGGAGGTGGCTCGCGCACACCTGGGCGGACTTTTTTCTTGTCGTGGAAAAGTGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGa                                                                                                           <  2:327133/149‑1 (MQ=255)
                                             tGGAGGTGAAGCTTACGAGGTGGGTTTAACTGCGCCGGAACTGCACCGCGGCGCGCTGGATGTGGGTTTTGGACACCTGGGCGTACCGTTTAACTGTATCGAGCAAATGACCGCCGATATTGAAGGGCTTTGCTTTATTTGTGACGGCGa                                                                                                     >  2:563543/1‑150 (MQ=255)
                                             ggggggTGAAGTTTAGGGTTTGTGTTTCCTTGGGCCGCAACTGACCGGTGGCGCGCTGGTTGTGGCTCTCGGCCCCGTGGGCGTACTTTTTTCTTGTAGGGAGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTg                                                                                                            <  2:638039/139‑1 (MQ=255)
                                             ggggggTAAAGATTACGATTTGTTTTCCGCTGTGCCGGAACTGACCGGGGGGCCGCTGGATGTGGCTCCCGGGCTTTTGGGCCATCGTTTACCTTGTATGGGGAAAAGGGCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg                                                                                                        <  1:613361/143‑1 (MQ=255)
                                              gggggTGAAGTTTACGGGTTGTGTTTCACTGTGCCGGAACTGAACCGGGGCGCGCTGGTTGTGGCTCTCGGACCCCGGGGCGTCCCTTTTCCCTGTTTGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                      <  1:464764/145‑1 (MQ=255)
                                              gggggTGAAGATTACGATTTGTGTTTCACTGTGCCGGAACGGAACCGTGGCGCGCTGGTTGTGGCTCTCGGACACGTGGGCTTCCCGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGaa                                                                                                    <  1:824260/147‑1 (MQ=255)
                                               ggggtagaGATTACGAGTTGGGTGTCAAGGCGGTGGAACGAACCCGGGGCGCGCTGATTGTGGCTGTGGAGCGCCGGGGGGTCCGGTTTCACTGCACAGGACAAACCGCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                      <  1:86448/140‑1 (MQ=255)
                                                     aaGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                            >  1:876443/1‑151 (MQ=255)
                                                              aGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                                   >  1:546476/1‑151 (MQ=255)
                                                                ttGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg                                                                                  >  2:564178/1‑150 (MQ=255)
                                                                   tgtTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaa                                                                                >  2:535146/1‑149 (MQ=255)
                                                                    tttgctgtgtgCCGAACATGACCCGGGGCGGGTGGGATGTTGCCCTCGGACCTCTGGGCGTCCGTTTTTGTTTATCGGGGAAAAGGTCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAg                                                                              <  2:648111/143‑1 (MQ=255)
                                                                     tttCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                             >  1:47710/1‑150 (MQ=255)
                                                                     tttCACGGTGCCGGACCTGAACGGGGGCGCGCTGGATGTGGCCCTCGGCCACCTGGGCGTTCCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAg                                                                              <  2:328470/149‑1 (MQ=255)
                                                                     ttcccttgtgCCGGAAGAGAGCGGGGGCGGGCGGGGTGGGGCTCGGGGACCCGGGGGCGTACCTGTTCCCTGTTTCGGCCAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                                                 <  1:449816/140‑1 (MQ=255)
                                                                      ttCACTGTGCGGGAACTGACCGGTGGCGCGCTGGTTGGGGCTCTCGGACCCTGGGGCGTCCCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGa                                                                            <  1:501875/150‑1 (MQ=255)
                                                                       tCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGat                                                                           >  1:341967/1‑150 (MQ=255)
                                                                         aCTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATg                                                                        >  1:811843/1‑151 (MQ=255)
                                                                             tgCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc                                                                      >  2:651151/1‑149 (MQ=255)
                                                                              gCCGGAACTGAACCGGGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTACCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAt                                                                   <  1:305429/151‑1 (MQ=255)
                                                                               cgggAACTGAACGGTGGCGCGCTGGTTGTGGCTCTGGGACCCCTGGGCGTCCCTTTTCCTTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCa                                                                    <  1:933070/147‑1 (MQ=255)
                                                                               ccGGTACTGACCCGTGGCGCGGTGGTGGGCGCTCGCGGTCCCGGGGTTGGACGTTTTACTTGTGTCGGACAATGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGata                                                                          <  2:370404/143‑1 (MQ=255)
                                                                                 ggAACTGACCCGGGGCGCGCGGGATGTGGCTCCCGCCCCCTGGGGCGTCCCTTTTCCCTTTATCGGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt                                                                <  2:883536/151‑1 (MQ=255)
                                                                                            ggggcgcgtgggtggtgGCTCGCGGACCCCGGGGCGTCCCGTTTCCCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCTAt                                                      <  1:817613/136‑1 (MQ=255)
                                                                                            gTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCAt                                                      >  2:727507/1‑150 (MQ=255)
                                                                                             gggCGCGTGGGATGTGGCTCTGGGACCCTGGGGCGTCCCTTTTCCCTGTCCCGGGAAAAGGACCGCCGATATCGAAGGGCTTTTTTTTTTTCTTTCCGTCGAACCTGTTACATTATCCTGGGAAGTATATGACCATTTTGCCACGCCCTaa                                                    <  1:791075/150‑1 (MQ=255)
                                                                                              ggCGCGGTGTATGTGGCTCTCGGACCCTGGGGCGTACTTTTTCCCTGTTGGGGAAAAGTGTCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCa                                                       <  2:229943/147‑1 (MQ=255)
                                                                                              ggCGCGGGGGGTGTGTCCCTCGGACCCTGGGGCGTCCTGTTTCCTTTTCTGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa                                                    <  2:366346/150‑1 (MQ=255)
                                                                                                   gttGGTTGTGGCTCTCGGACCCCGGGGCGTCCCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATg                                               <  1:243933/148‑1 (MQ=255)
                                                                                                   gtgggATGTGTCTCTCGGACACTTGTTCGTCCCTTTTCCCTGTTTGGGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCa                                                       <  2:563122/139‑1 (MQ=255)
                                                                                                               tctcGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg                                  >  1:100169/1‑151 (MQ=255)
                                                                                                                cccGGACCCCGGGGCGTCCCTTTTCCCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCg                                  <  1:582517/148‑1 (MQ=255)
                                                                                                                   gccccccggggCGTCCCGTTTCCTTGTTGGGGGAAAATTACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCt                               <  1:860009/142‑1 (MQ=255)
                                                                                                                    gacccggggggGTACCTTTTCCCTGTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTc                                   <  1:765311/138‑1 (MQ=255)
                                                                                                                     acacCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                             >  2:208426/1‑150 (MQ=255)
                                                                                                                      cacCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                             >  1:472972/1‑149 (MQ=255)
                                                                                                                        ccTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                             >  2:714938/1‑147 (MQ=255)
                                                                                                                             ggcgtcctttttccgtGTACCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGt                     <  1:259387/136‑1 (MQ=255)
                                                                                                                             gCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGt                     >  1:43154/1‑150 (MQ=255)
                                                                                                                                     tttCCCTGTACGGGCCAATTGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg             <  2:912888/150‑1 (MQ=255)
                                                                                                                                      tttccttGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg             <  1:59459/144‑1 (MQ=255)
                                                                                                                                      ttACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGc           >  2:535661/1‑151 (MQ=255)
                                                                                                                                        tccttGTATGGGCAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGc           <  1:283786/146‑1 (MQ=255)
                                                                                                                                            tGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCAt         <  1:32152/148‑1 (MQ=255)
                                                                                                                                                   gggAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACtt   <  1:502187/148‑1 (MQ=255)
                                                                                                                                                   gggAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTg  <  2:917733/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTG  >  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: