Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913436,5160CA44.7% 102.9 / 35.5 77Q295K (CAA→AAA) thiLthiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (31/11);  new base A (0/34);  total (31/46)
Fisher's exact test for biased strand distribution p-value = 2.36e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTG  >  NC_000913/436367‑436661
                                                                                                                                                        |                                                                                                                                                 
tttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGc                                                                                                                                                   >  1:327525/1‑150 (MQ=255)
 ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGc                                                                                                                                                   <  2:768462/149‑1 (MQ=255)
     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAt                                                                                                                                               <  1:49227/149‑1 (MQ=255)
     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAt                                                                                                                                               <  1:117262/149‑1 (MQ=255)
     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAt                                                                                                                                               <  2:117623/149‑1 (MQ=255)
     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGa                                                                                                                                             <  1:897475/151‑1 (MQ=255)
        cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccg                                                                                                                                          >  1:843418/1‑151 (MQ=255)
           gTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccgccg                                                                                                                                       >  1:441685/1‑151 (MQ=255)
            ttGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccgccg                                                                                                                                       <  2:614483/150‑1 (MQ=255)
              gAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGat                                                                                                                                     >  2:530448/1‑150 (MQ=255)
                aCCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATc                                                                                                                                  >  2:1023874/1‑151 (MQ=255)
                   ggAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGa                                                                                                                                >  2:273138/1‑150 (MQ=255)
                      aCAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGAGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAggg                                                                                                                            >  1:396654/1‑151 (MQ=255)
                       cAGGCGCTGCGCTGGGCTCTCTCGGGGGGTGAAGTTTACGTGTTGTTTTTCACTGTGCCGAAACTGACCCGTGGCGCGCGGTATGGGGCCCCCGGCCCCCTGGGCGTACCGTTTACCTGTATCGGGAAAATGACCGCCGATATCGAAGGGc                                                                                                                           <  2:233934/151‑1 (MQ=255)
                       cAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGc                                                                                                                           >  2:1039797/1‑151 (MQ=255)
                         ggCGCTGCGTGGGGCGCTCTCGGGCGGTGAAGATTACGAGTTGTGTTTCCCTGGGCCGGAACTGACCCGGGGCGCGCTGGATGTGGCTCTCGGACACCGGGGCGTCCCGTTTCCCTGTACCGGGAAAATGACCGCCGATATCGAAGGGc                                                                                                                           <  1:287447/149‑1 (MQ=255)
                          gcgctgcgctGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCttt                                                                                                                        >  1:1046800/1‑151 (MQ=255)
                            gctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTg                                                                                                                       >  2:182443/1‑150 (MQ=255)
                            gctgcgcggGGCGCTCTCTGGCGGTGAAGTTTACGGGTTGTGTTTCCCTGTGCCGGAACTGACCCGTGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCTTACCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCttt                                                                                                                        <  1:1022422/149‑1 (MQ=255)
                               gcgcggGGCGCTCTGTGGGGGTAAAGATTACGTGTTGTTTTTCACTGTGCCGAACCGGAACCGGGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCGTCCCGTTTACTTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGtt                                                                                                                     <  1:126984/149‑1 (MQ=255)
                                    gggCGCTCTCTGGCGGTGAAGATTACGTGTTGTGTTTCGCTGTGCGGAACCTGAACCGGGGCGCGCTGGATGTGGCTCTCGGACACGTGGGCGTCCTTTTTTCTTGTATCGGACAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCg                                                                                                              <  2:344233/151‑1 (MQ=255)
                                    gggCGCTCCCGGGGGGGGAAGATTACGGGTTGTGTTCCACGGTGCCGGAACTGAACCGTGGCTCGGTGGATGGGTCTCCCGGCCACCTGGGCGTCCCTTTTCCCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTc                                                                                                               <  2:829136/150‑1 (MQ=255)
                                    gggCGCGCTCGGGCGGTGAAGTTTACGAGTTGTGTTTCAGCTTCCCGGAACCGACCTGTGGGGCGCGGGATGTGGCGCTGGAACACGTGGCCGGACTGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCtt                                                                                                                         <  2:21332/140‑1 (MQ=255)
                                     ggCTCTTTCTGGGGGAGAGGATGACGGGTTGTGTTTCTCTGTGCGGGAACTGAACGGGGGCGCGCTGGGTGGGGCTCACGAACCCCGGGGAATACCTTGTCCTTGTTTCGGTAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGt                                                                                                             <  1:835099/151‑1 (MQ=255)
                                          tttgtGGGGGTAAAGATTATGGTTTGTGTTACGCTGGGCCGGAATTAAACCGTGGCGCGCTGGATGTGCCTCTCGGCCACCTGGGCGTCCCGTTTCCTTGTATCGGCAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGAc                                                                                                          <  1:676509/145‑1 (MQ=255)
                                               gggggTGAAGATTAGGTGTTGTGTTTCTCTGTGCCGGAACTGAACCGTGGCGCGTTGGTTGTGGCTCCCGGCCACGTGGGCGTTCCTTTTTCCTGTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                      <  2:726525/145‑1 (MQ=255)
                                               gggggTAAAGATTAGGAGTTGTGTTTCCCTGTGCCGAAACTGACCCGGGGCGCGCTGGATGGGGCTCCCGGCCCCCGGGGCGTGCTTTTTCCTTGTACGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg                                                                                                        <  2:321924/143‑1 (MQ=255)
                                                  ggTGAAGTTTAGGATTTGTGTTTCACGGTGCCGTAGCGAACCCGGGGGGGGTGGGGTTGGCCTCGGGGACACCGGGGCGTCCTTTTTACTTTTTGCGGGAAATTGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt                                                                                                 <  2:727511/150‑1 (MQ=255)
                                                  ggTAAAGATTAGGGGTTGTGTTTCCTTGTGCGGGAACTGAACCGTGGCGGGCTGGTTGTGGCTCTGGGACCCCGGGGCGTCCCTTTTACCTGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt                                                                                                 <  1:656858/150‑1 (MQ=255)
                                                   gtaaagttgcggatttGTGTTTCACTGTGCCGAAAATGACCCGGGGCGTGTTATGTGGGGCCCCGGGAAGCGTGGGTGTTTTTTTTCCCTGTACCGGCCAATAGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                      <  2:576772/130‑1 (MQ=255)
                                                    tGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGtt                                                                                              >  1:1068873/1‑151 (MQ=255)
                                                    tGAAGATTACGAGTTGGGTTTCACTGTGTCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGAGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTAGTGACGGCGAACCTGtt                                                                                              >  2:37178/1‑151 (MQ=255)
                                                      aaGTTTACGTGTTGTTTTTCCCTGTGCCGAAACTGACCGGTGGCGCGCTGGATGTGGCTCTCGGCCACCGGGGCGTCCCTTTTCCCTGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTa                                                                                             <  1:1010935/150‑1 (MQ=255)
                                                       aGATTAGGTGTTGTTTTTCACTGGGCCGGAACTGACCGGTGGCGCGCTGGATGTGGCCCTCGGACACCTGGGCGTCCTTTTTTCCTGTACCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACa                                                                                           <  2:539011/151‑1 (MQ=255)
                                                       aGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                            >  2:938876/1‑150 (MQ=255)
                                                        gATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAt                                                                                          >  2:265747/1‑151 (MQ=255)
                                                          ttACGTGTTGTTTTGCCCTGTGCGGGAACGAACCGGGGGCGCGCTGGGTTTGGCCCGCTGACACCTGGGCGCCCTGTTTCCCTGTTCCGGAAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAt                                                                                          <  2:951963/149‑1 (MQ=255)
                                                           taggggtggtttttCCCTGGGCCGAAACTGACCGGGGGCGCGTGGGTTGTGGCTCTCGAACCCGTGGGCTTACTGTTTCCTTGTAGGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt                                                                                         <  2:330685/138‑1 (MQ=255)
                                                              gttttGTGTTTCCCGGTGCCGAACTTGACCGGTGGCGCGGGGGGTGTGTCTCTGGCACTCGGGGGCGTACTTTTTCCTTGTTGCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAg                                                                                       <  2:638726/145‑1 (MQ=255)
                                                                 ttGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg                                                                                  >  1:823959/1‑150 (MQ=255)
                                                                 ttGTGTTTACTTGTGCCGGAACTAAGCGTGGCCGCGTTGGAGGGGTCTCCGGGACACTTGGAGGATCTGTTTTCTGGACGCGGGAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGa                                                                                      <  2:872312/146‑1 (MQ=255)
                                                                      tttCCCTGTGCCGAAATTGAGGCGGGGGGCGGTGGGTTGGCCTCCGGGACGCTTGGGCGTACTTTTTTCTTGTCGGGGAGAAGGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAg                                                                              <  2:501052/149‑1 (MQ=255)
                                                                      tttCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTCCCTTTTCCCTGTTTCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAg                                                                              <  1:100935/149‑1 (MQ=255)
                                                                         ccctgtccgggacgtGACCCGTGGCGGGGTGGATTTGCCCCCCGGCCACGGGGCCGTGCTTTTTCCCGGTAGCGGGAAAGTGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaa                                                                                <  2:951521/130‑1 (MQ=255)
                                                                         cACTGTGCCGAAACTGACCCGTGGCGCGCTGGATGGGCCTCTCGGACCCCGGGGCGTCCCTTTTCCCTGTTTCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGatat                                                                         <  1:740261/151‑1 (MQ=255)
                                                                          aatgtgCCGGAACTGACGGGTGGCGCGGTGAATGTGTGTCCCGGACGCTTGGGCGTCCTGTTTCCCTGTAGGGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATg                                                                        <  2:34935/149‑1 (MQ=255)
                                                                            tgtgCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc                                                                      >  1:816898/1‑151 (MQ=255)
                                                                                ccGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCAGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtt                                                                  >  2:514323/1‑151 (MQ=255)
                                                                                   gAACTGAGCCGTGGCGCGCTGGATGTGGCTCGCGGACCCCTGGGCGTTCCGTTTCCCTGTAGCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAttt                                                                 <  2:396445/149‑1 (MQ=255)
                                                                                   gAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGCTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTg                                                               >  2:738450/1‑151 (MQ=255)
                                                                                    aaCTGACCGGTGGCGCGCTGGATGTGGCTCTCGGACCCGGGGGCGTCCCTTTTCCTTGTTTCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTg                                                               <  1:682682/150‑1 (MQ=255)
                                                                                    aaCTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                              >  1:711728/1‑151 (MQ=255)
                                                                                     aCTGACCGGGGGCGCGCTGGATGTGGCTCTCGGCCCCGGGGGCGTCCCGTTTACCTGTGTGGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGcc                                                             <  1:315630/151‑1 (MQ=255)
                                                                                          aCCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGcc                                                        >  1:175981/1‑151 (MQ=255)
                                                                                           ccGGGGCGCGCTGGTTGTGGCTCTCGGACCCCTGGGCTTCCCTTTTCCCTGTAGCGGGAAATTGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                                          <  1:642716/148‑1 (MQ=255)
                                                                                               tGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATa                                                     >  2:706513/1‑150 (MQ=255)
                                                                                                ggCGCGTTGGTTGTGGCTCTGGGACCCGGGGGCTTCCTTTTTCCCTGTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa                                                    <  1:506156/150‑1 (MQ=255)
                                                                                                ggCGCGCGGGTTGTGGCTCTCGGACACCGGGGCGTCCCGTTTACCTGTATGGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCAt                                                      <  2:843103/148‑1 (MQ=255)
                                                                                                    gcgcTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa                                                    >  2:78894/1‑147 (MQ=255)
                                                                                                    gcgcTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa                                                    >  1:158489/1‑147 (MQ=255)
                                                                                                             gtgGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaa                                        >  2:301847/1‑150 (MQ=255)
                                                                                                                  cgcgcggcgcttggggCTTACTTTTTCCCTGTATGTAGAAAAGGTCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaa                                        <  2:591104/133‑1 (MQ=255)
                                                                                                                       ggACCCCGGGGCGTTCGTTTTCCCTTTATGGGAAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTc                                   <  1:33828/146‑1 (MQ=255)
                                                                                                                       ggACACTGGGGCGTGCCGTTTCCCTTTATGGGGAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTc                                   <  1:370213/146‑1 (MQ=255)
                                                                                                                           acCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGa                             >  1:538586/1‑148 (MQ=255)
                                                                                                                               gggCGTCCCTTTTCCTTGTATGGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATg                        <  1:734806/149‑1 (MQ=255)
                                                                                                                                ggCGTCCCTTTTCCTTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGt                     <  1:1037857/151‑1 (MQ=255)
                                                                                                                                ggCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGt                     >  1:937922/1‑151 (MQ=255)
                                                                                                                                   gTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGCGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAAt                  >  2:167459/1‑151 (MQ=255)
                                                                                                                                      ctttttACCTGTAGCGGAAAAATTCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATcc                <  2:711441/147‑1 (MQ=255)
                                                                                                                                      cctttatccTGTACGGGGCAATGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTaa                   <  2:243323/140‑1 (MQ=255)
                                                                                                                                        gTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg             >  2:743527/1‑151 (MQ=255)
                                                                                                                                         tttCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATcc                <  1:1024229/147‑1 (MQ=255)
                                                                                                                                            aCCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCa          >  1:621238/1‑150 (MQ=255)
                                                                                                                                            aCCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGCGGCAt         >  2:320482/1‑151 (MQ=255)
                                                                                                                                                gTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTAc     >  1:920607/1‑151 (MQ=255)
                                                                                                                                                 tATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACt    >  2:278118/1‑151 (MQ=255)
                                                                                                                                                     gggAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACt    <  1:530704/147‑1 (MQ=255)
                                                                                                                                                     gggAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTg  <  1:589599/149‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                 
TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTG  >  NC_000913/436367‑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: