Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913436,5160CA40.7% 88.3 / 29.8 60Q295K (CAA→AAA) thiLthiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (29/6);  new base A (0/24);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 3.29e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCT  >  NC_000913/436367‑436663
                                                                                                                                                        |                                                                                                                                                     
tttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGc                                                                                                                                                       >  2:973921/1‑150 (MQ=255)
   ctcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAt                                                                                                                                                  >  2:699799/1‑151 (MQ=255)
     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑g                                                                                                                                                <  1:140318/150‑1 (MQ=255)
        cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GAccg                                                                                                                                            <  1:794250/151‑1 (MQ=255)
        cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GAcc                                                                                                                                             >  2:1189790/1‑150 (MQ=255)
         aTGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GAccgc                                                                                                                                           >  1:631731/1‑151 (MQ=255)
                aCCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATc                                                                                                                                    >  1:747060/1‑151 (MQ=255)
                 ccGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCg                                                                                                                                   >  2:264055/1‑151 (MQ=255)
                              tgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGttt                                                                                                                      >  1:387741/1‑151 (MQ=255)
                               gggttgGGCGCTCTCTGGGGGTGAGGTTTAGGGGTTGTGTTTCACGGTGGGGGAACTGAACCGGGGCGCGC‑TGGTTGTGGCTCGC‑GGACTCCGGGG‑TGTCCCGTTTTCCTTTTTGGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGtt                                                                                                                       <  1:607336/145‑1 (MQ=255)
                                 gctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGtttt                                                                                                                     >  2:535596/1‑149 (MQ=255)
                                 gctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                                                   >  1:341851/1‑151 (MQ=255)
                                  ctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTAtt                                                                                                                  >  2:388675/1‑151 (MQ=255)
                                   tgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTc                                                                                                                 >  2:979187/1‑151 (MQ=255)
                                           ctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg                                                                                                          >  1:147922/1‑150 (MQ=255)
                                              tGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGa                                                                                                       >  1:187133/1‑150 (MQ=255)
                                              tGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGa                                                                                                       >  1:1134244/1‑150 (MQ=255)
                                               gggggTGAAGATTGGGTGTTGTTTTTCCCTGGGCCGGACCTGACCGGGGGCGCGT‑TGGTGGTGGCTCTC‑GGCCACGGGGG‑CGTACTTTTTCCCTGTAGCGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGa                                                                                                       <  2:620911/146‑1 (MQ=255)
                                               gggggTGAAGATTACGAGTTGTTTTTCACTGTGCCGAAACTGAACGGTGGCGCGC‑TGGATGTGGCTCCC‑GGACCCCTGGG‑CGTCCCTTTTCCCTGTATGGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg                                                                                                          <  1:81881/143‑1 (MQ=255)
                                                       aGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTgat                                                                                                >  2:625194/1‑146 (MQ=255)
                                                       aGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                              >  2:103492/1‑150 (MQ=255)
                                                         aTTAGGGGTTTTGTTTCCCGGGGCCGAATCTGACCGGTGGCCCGT‑GGGATGTGGCTCTC‑GGACCCTGGGG‑CTTCCTTTTTACCTGTATGGGGA‑AAGG‑GTCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTa                                                                                               <  1:853603/147‑1 (MQ=255)
                                                           tACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt                                                                                           >  1:1097952/1‑149 (MQ=255)
                                                              gttttGTTTTTCCCTGGGCGCAAACTGAGCCGGGGCGCGC‑TGGTTGTGCCTCCC‑GGACGCCGGGG‑CGTACCTTTTCCCTGTATCGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACt                                                                                      <  2:346362/148‑1 (MQ=255)
                                                              gAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACt                                                                                      >  2:1004727/1‑151 (MQ=255)
                                                               aGTTGTGTTTCACTGTGCCGGGACTGAAACGTGGCGCGC‑TGGATGGGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGCGACGGCGAACCTGTTACATTAGACTg                                                                                     >  2:709301/1‑151 (MQ=255)
                                                                 ttttttttCACTGTGCGGGAACTAACCGGGGGCGCGT‑TGGTTGGGTCTCTC‑GGACACGTGGG‑CTTCCCTTTTACCTGTAGGGGAA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg                                                                                    <  2:776940/145‑1 (MQ=255)
                                                                  tgttTTTCCCTGTGCGGGAACTAACCGGTGGGGCGG‑TGGATGTGGCTCTG‑GGACACTGGGG‑CGTACGTTTTACCTGTATCGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                                     <  2:1122215/148‑1 (MQ=255)
                                                                    tttttcctggggCCGGAAGTGATCGGGGGCGCTG‑GGGATGGGCCCCCC‑GGACCCCTGGG‑CGTACTTTTTACCTTTAGCGGAC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaa                                                                                  <  1:787460/138‑1 (MQ=255)
                                                                      tttCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                               >  1:727343/1‑150 (MQ=255)
                                                                        tcctggtgCCGGAAC‑GGACCGTGGGCCGCGTGGATTTTGC‑CTCGGGACACCTGCG‑TGTCCTTTTTCCCTGTAGGGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAg                                                                                <  2:794778/142‑1 (MQ=255)
                                                                         cACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGata                                                                            >  2:1105878/1‑150 (MQ=255)
                                                                              tgCCGAAACAGACCGGTGGCGCGC‑TGGTTGTGCCTCTC‑GGCCACGTGGG‑TGTCCTTTTTTCCTGTTGCGGAA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATg                                                                          <  1:832892/147‑1 (MQ=255)
                                                                               gcgggAATTGACCCGTGGCGCGG‑TGGTTGTGGCTCCC‑GGCCCCGTGGG‑CGTACTGTTT‑CCTTTTACGGGCAAAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAt                                                                     <  2:1134266/148‑1 (MQ=255)
                                                                                ccGGAACTGACCCGGGGCGCGC‑TGGATGTGTCTCTC‑GGACACCGGGG‑CGTCCCTTTTCCCTGTATCGGGC‑AAAG‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAt                                                                     <  1:944979/150‑1 (MQ=255)
                                                                                  ggACCTGACCGGGGGCGCGT‑TGGATGTGGCTTTG‑GGACCCCGGGG‑TGTCTTTTTTCCCTGTGT‑GGGA‑AAATGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGatat                                                                           <  2:746078/142‑1 (MQ=255)
                                                                                  ggAACTGAACGGGGGCGCGC‑TGGTTGTGGCTCTC‑GGACCCCGGGG‑GGTCCTGTTTACCTGTATGGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt                                                                  <  2:984609/151‑1 (MQ=255)
                                                                                  ggAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGAATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt                                                                  >  1:137559/1‑151 (MQ=255)
                                                                                  ggAACTGAACCGGGGCGCGC‑TGGGTGTGGCTCTG‑GGCCCCTTGGG‑CGTCCTTTTTTCCTTTATGGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc                                                                        <  2:878632/145‑1 (MQ=255)
                                                                                    acctgaccgggggCGCGC‑GGGTTGGGGCTCCC‑GGACCCCGGGG‑CGTCCGTTTTCCCTTTATGGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGata                                                                            <  2:372462/128‑1 (MQ=255)
                                                                                    aaCTGACCGGTGGCGCGC‑TGGGTGTGGCTCTC‑GGACACCGGGG‑CGTCCCTTTTCCCTGTATGGGGA‑AAAG‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                                <  1:494968/151‑1 (MQ=255)
                                                                                         agccGTGGCGCGT‑GGGATGTGGCTCTC‑GGACACGTGGG‑GGTGCCTTTTCCCTGTATCGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                                            <  1:379097/148‑1 (MQ=255)
                                                                                         acccGGGGCGCGT‑TGGATGTGGCTCTC‑GGACCCGGGGG‑CGTCCCTTTTCCCTGTATCGGGA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCa                                                              <  1:293880/146‑1 (MQ=255)
                                                                                         aaCCGAGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGc                                                           >  1:375848/1‑151 (MQ=255)
                                                                                                        gggtgtgGCTCCC‑GGCCCCCGGGG‑CTTCCCTTTTCCTTGTATGGGAA‑AAAT‑AACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTc                                               <  1:914363/146‑1 (MQ=255)
                                                                                                        ggATGTGGCTCTC‑GGCCCCCGGGG‑CGTCCCTTTTCCCTTTATGGGAA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcgc                                             <  1:851153/151‑1 (MQ=255)
                                                                                                          aTGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcgcg                                            >  2:521033/1‑150 (MQ=255)
                                                                                                                  ctc‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGc                                   >  2:607737/1‑151 (MQ=255)
                                                                                                                   cacggaccctgggg‑CGTACTTTTTTCCTGTGCGGGGA‑AAAG‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTc                                     <  1:1069164/138‑1 (MQ=255)
                                                                                                                       ggacccTGGGCCGTACCTTTTCCCTTTATCTGCA‑AAAT‑AACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTc                                               <  2:1156665/132‑1 (MQ=255)
                                                                                                                         ccctgggg‑CGTCCCTTTTCCTTTTATGGGAA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTc                                     <  1:453976/138‑1 (MQ=255)
                                                                                                                         ccccgggg‑CGTCCTTTTTCCCTGTACCGGGA‑AAAG‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGaa                              <  2:387726/145‑1 (MQ=255)
                                                                                                                                    tACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATc                   >  1:340307/1‑151 (MQ=255)
                                                                                                                                         tttccctgttgggggA‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg               <  1:363776/138‑1 (MQ=255)
                                                                                                                                          tccccTGTACGGGGA‑AAGT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg               <  2:236530/146‑1 (MQ=255)
                                                                                                                                            aCCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCAt           >  1:360344/1‑151 (MQ=255)
                                                                                                                                              cTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGCAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGGGGCATc          >  2:1112425/1‑150 (MQ=255)
                                                                                                                                               tGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCt         >  1:803282/1‑150 (MQ=255)
                                                                                                                                                 tATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACt      >  1:1243447/1‑151 (MQ=255)
                                                                                                                                                  aTCTGGG‑AAAT‑GACGGCCGATATCGAAGGGCTTTGATTTATTCTTGATGGCGAACATGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATTTCGCGAAAAGTCGTGTGAAGATGAGTAATCGGTGGCATCTACtt     >  2:174623/1‑151 (MQ=255)
                                                                                                                                                     gggA‑AAAT‑GGCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCt  <  1:1105859/151‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGC‑TGGATGTGGCTCTC‑GGACACCTGGG‑CGTACCGTTTACCTGTATCGGGC‑AAAT‑GACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCT  >  NC_000913/436367‑436663

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