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

TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTG  >  NC_000913/436367‑436661
                                                                                                                                                        |                                                                                                                                                 
tttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGc                                                                                                                                                   >  1:327525/1‑150 (MQ=255)
 ttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGc                                                                                                                                                   <  2:768462/149‑1 (MQ=255)
     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAt                                                                                                                                               <  1:49227/149‑1 (MQ=255)
     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAt                                                                                                                                               <  1:117262/149‑1 (MQ=255)
     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAt                                                                                                                                               <  2:117623/149‑1 (MQ=255)
     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGa                                                                                                                                             <  1:897475/151‑1 (MQ=255)
        cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccg                                                                                                                                          >  1:843418/1‑151 (MQ=255)
           gTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccgccg                                                                                                                                       >  1:441685/1‑151 (MQ=255)
            ttGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccgccg                                                                                                                                       <  2:614483/150‑1 (MQ=255)
              gAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGat                                                                                                                                     >  2:530448/1‑150 (MQ=255)
                aCCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATc                                                                                                                                  >  2:1023874/1‑151 (MQ=255)
                   ggAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGa                                                                                                                                >  2:273138/1‑150 (MQ=255)
                      aCAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGAGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAggg                                                                                                                            >  1:396654/1‑151 (MQ=255)
                       cAGGCGCTGCGCTGGGCTCTCTCGGGGGGTGAAGTTTACGTGTTGTTTTTCACTGTGCCGAAACTGACCC‑GTGG‑CGCGCGGTATGGGG‑CCCCCGGCCCCCTGGGCGTACCGTTTACCTGTATCGGGAAAATGACCGCCGATATCGAAGGGc                                                                                                                           <  2:233934/151‑1 (MQ=255)
                       cAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGc                                                                                                                           >  2:1039797/1‑151 (MQ=255)
                         ggCGCTGCGTGGGGCGCTCTCGGGCGGTGAAGATTACGAGTTGTGTTTCCCTGGGCCGGAACTGACCC‑GGGG‑CGCGCTGGATGTGG‑CTCTCGGACACCGGGGCGTCCCGTTTCCCTGTACCGGGAAAATGACCGCCGATATCGAAGGGc                                                                                                                           <  1:287447/149‑1 (MQ=255)
                          gcgctgcgctGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCttt                                                                                                                        >  1:1046800/1‑151 (MQ=255)
                            gctgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTg                                                                                                                       >  2:182443/1‑150 (MQ=255)
                            gctgcgcggGGCGCTCTCTGGCGGTGAAGTTTACGGGTTGTGTTTCCCTGTGCCGGAACTGACCC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACCCCGGGGCTTACCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCttt                                                                                                                        <  1:1022422/149‑1 (MQ=255)
                               gcgcggGGCGCTCTGTGGGGGTAAAGATTACGTGTTGTTTTTCACTGTGCCGAACCGGAACC‑GGGG‑CGCGCTGGATGTGG‑CTCTCGGACCCCGGGGCGTCCCGTTTACTTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGtt                                                                                                                     <  1:126984/149‑1 (MQ=255)
                                    gggCGCTCTCTGGCGGTGAAGATTACGTGTTGTGTTTCGCTGTGCGGAACCTGAACC‑GGGG‑CGCGCTGGATGTGG‑CTCTCGGACACGTGGGCGTCCTTTTTTCTTGTATCGGACAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCg                                                                                                              <  2:344233/151‑1 (MQ=255)
                                    gggCGCTCCCGGGGGGGGAAGATTACGGGTTGTGTTCCACGGTGCCGGAACTGAACC‑GTGG‑CTCGGTGGATGGGT‑CTCCCGGCCACCTGGGCGTCCCTTTTCCCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTc                                                                                                               <  2:829136/150‑1 (MQ=255)
                                    gggCGCGCTCGGGCGGTGAAGTTTACGAGTTGTGTTTCAGCTTCCCGGAACCGACCT‑GTGG‑GGCGCGGGATGTGG‑CGCTGGAACACGTGGCCGGACTGTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCtt                                                                                                                         <  2:21332/140‑1 (MQ=255)
                                     ggCTCTTTCTGGGGGAGAGGATGACGGGTTGTGTTTCTCTGTGCGGGAACTGAACG‑GGGG‑CGCGCTGGGTGGGG‑CTCACGAACCCCGGGGAATACCTTGTCCTTGTTTCGGTAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGt                                                                                                             <  1:835099/151‑1 (MQ=255)
                                          tttgtGGGGGTAAAGATTATGGTTTGTGTTACGCTGGGCCGGAATTAAACC‑GTGG‑CGCGCTGGATGTGC‑CTCTCGGCCACCTGGGCGTCCCGTTTCCTTGTATCGGCAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGAc                                                                                                          <  1:676509/145‑1 (MQ=255)
                                               gggggTGAAGATTAGGTGTTGTGTTTCTCTGTGCCGGAACTGAACC‑GTGG‑CGCGTTGGTTGTGG‑CTCCCGGCCACGTGGGCGTTCCTTTTTCCTGTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                      <  2:726525/145‑1 (MQ=255)
                                               gggggTAAAGATTAGGAGTTGTGTTTCCCTGTGCCGAAACTGACCC‑GGGG‑CGCGCTGGATGGGG‑CTCCCGGCCCCCGGGGCGTGCTTTTTCCTTGTACGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg                                                                                                        <  2:321924/143‑1 (MQ=255)
                                                  ggTGAAGTTTAGGATTTGTGTTTCACGGTGCCGTAGC‑GAACC‑CGGG‑GGGGGTGGGGTTGGCCTCGGGGACACCGGGGCGTCCTTTTTACTTTTTGCGGGAAATTGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt                                                                                                 <  2:727511/150‑1 (MQ=255)
                                                  ggTAAAGATTAGGGGTTGTGTTTCCTTGTGCGGGAACTGAACC‑GTGG‑CGGGCTGGTTGTGG‑CTCTGGGACCCCGGGGCGTCCCTTTTACCTGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt                                                                                                 <  1:656858/150‑1 (MQ=255)
                                                   gtaaagttgcggatttGTGTTTCACTGTGCCGAAAATGACCC‑GGGG‑CGTGTTATGTGGGG‑CCCCGGGAAGCGTGGGTGTTTTTTTTCCCTGTACCGGCCAATAGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                      <  2:576772/130‑1 (MQ=255)
                                                    tGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGtt                                                                                              >  1:1068873/1‑151 (MQ=255)
                                                    tGAAGATTACGAGTTGGGTTTCACTGTGTCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGAGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTAGTGACGGCGAACCTGtt                                                                                              >  2:37178/1‑151 (MQ=255)
                                                      aaGTTTACGTGTTGTTTTTCCCTGTGCCGAAACTGACCG‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGCCACCGGGGCGTCCCTTTTCCCTGTATCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTa                                                                                             <  1:1010935/150‑1 (MQ=255)
                                                       aGATTAGGTGTTGTTTTTCACTGGGCCGGAACTGACCG‑GTGG‑CGCGCTGGATGTGG‑CCCTCGGACACCTGGGCGTCCTTTTTTCCTGTACCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACa                                                                                           <  2:539011/151‑1 (MQ=255)
                                                       aGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                            >  2:938876/1‑150 (MQ=255)
                                                        gATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAt                                                                                          >  2:265747/1‑151 (MQ=255)
                                                          ttACGTGTTGTTTTGCCCTGTGCGGGAAC‑GAACCGGGGG‑CGCGCTGGGTTTGG‑CCCGCTGACACCTGGGCGCCCTGTTTCCCTGTTCCGGAAAAATGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAt                                                                                          <  2:951963/149‑1 (MQ=255)
                                                           taggggtggtttttCCCTGGGCCGAAACTGACCG‑GGGG‑CGCGTGGGTTGTGG‑CTCTCGAACCCGTGGGCTTACTGTTTCCTTGTAGGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt                                                                                         <  2:330685/138‑1 (MQ=255)
                                                              gttttGTGTTTCCCGGTGCCGAACTTGACCG‑GTGG‑CGCGGGGGGTGTGT‑CTCTGGCACTCGGGGGCGTACTTTTTCCTTGTTGCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAg                                                                                       <  2:638726/145‑1 (MQ=255)
                                                                 ttGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg                                                                                  >  1:823959/1‑150 (MQ=255)
                                                                 ttGTGTTTACTTGTGCCGGAACT‑AAGC‑GTGGCCGCGTTGGAGGGGT‑CTCCGGGACACTTGGAGGATCTGTTTTCTGGACGCGGGAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGa                                                                                      <  2:872312/146‑1 (MQ=255)
                                                                      tttCCCTGTGCCGAAATTGAGGC‑GGGG‑GGCGGTGGGTTGGC‑CTCCGGGACGCTTGGGCGTACTTTTTTCTTGTCGGGGAGAAGGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAg                                                                              <  2:501052/149‑1 (MQ=255)
                                                                      tttCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTCCCTTTTCCCTGTTTCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAg                                                                              <  1:100935/149‑1 (MQ=255)
                                                                         ccctgtccgggacgtGACCC‑GTGG‑CGGGGTGGATTTGC‑CCCCCGGCCACGGGGCCGTGCTTTTTCCCGGTAGCGGGAAAGTGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaa                                                                                <  2:951521/130‑1 (MQ=255)
                                                                         cACTGTGCCGAAACTGACCC‑GTGG‑CGCGCTGGATGGGC‑CTCTCGGACCCCGGGGCGTCCCTTTTCCCTGTTTCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGatat                                                                         <  1:740261/151‑1 (MQ=255)
                                                                          aatgtgCCGGAACTGACGG‑GTGG‑CGCGGTGAATGTGT‑GTCCCGGACGCTTGGGCGTCCTGTTTCCCTGTAGGGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATg                                                                        <  2:34935/149‑1 (MQ=255)
                                                                            tgtgCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc                                                                      >  1:816898/1‑151 (MQ=255)
                                                                                ccGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCAGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtt                                                                  >  2:514323/1‑151 (MQ=255)
                                                                                   gAACTGAGCC‑GTGG‑CGCGCTGGATGTGG‑CTCGCGGACCCCTGGGCGTTCCGTTTCCCTGTAGCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAttt                                                                 <  2:396445/149‑1 (MQ=255)
                                                                                   gAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGCTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTg                                                               >  2:738450/1‑151 (MQ=255)
                                                                                    aaCTGACCG‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACCCGGGGGCGTCCCTTTTCCTTGTTTCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTg                                                               <  1:682682/150‑1 (MQ=255)
                                                                                    aaCTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                              >  1:711728/1‑151 (MQ=255)
                                                                                     aCTGACCG‑GGGG‑CGCGCTGGATGTGG‑CTCTCGGCCCCGGGGGCGTCCCGTTTACCTGTGTGGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGcc                                                             <  1:315630/151‑1 (MQ=255)
                                                                                          aCC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGcc                                                        >  1:175981/1‑151 (MQ=255)
                                                                                           cc‑GGGG‑CGCGCTGGTTGTGG‑CTCTCGGACCCCTGGGCTTCCCTTTTCCCTGTAGCGGGAAATTGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                                          <  1:642716/148‑1 (MQ=255)
                                                                                               tGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATa                                                     >  2:706513/1‑150 (MQ=255)
                                                                                                gg‑CGCGTTGGTTGTGG‑CTCTGGGACCCGGGGGCTTCCTTTTTCCCTGTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa                                                    <  1:506156/150‑1 (MQ=255)
                                                                                                gg‑CGCGCGGGTTGTGG‑CTCTCGGACACCGGGGCGTCCCGTTTACCTGTATGGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCAt                                                      <  2:843103/148‑1 (MQ=255)
                                                                                                    gcgcTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa                                                    >  2:78894/1‑147 (MQ=255)
                                                                                                    gcgcTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATaa                                                    >  1:158489/1‑147 (MQ=255)
                                                                                                             gtgG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGaa                                        >  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)
                                                                                                                                                        |                                                                                                                                                 
TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGG‑CGCGCTGGATGTGG‑CTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTG  >  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: ‑