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

TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCT  >  NC_000913/436367‑436663
                                                                                                                                                        |                                                                                                                                                     
tttCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGc                                                                                                                                                       >  2:973921/1‑150 (MQ=255)
   ctcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAAt                                                                                                                                                  >  2:699799/1‑151 (MQ=255)
     cGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATg                                                                                                                                                <  1:140318/150‑1 (MQ=255)
        cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccg                                                                                                                                            <  1:794250/151‑1 (MQ=255)
        cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAcc                                                                                                                                             >  2:1189790/1‑150 (MQ=255)
         aTGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGAccgc                                                                                                                                           >  1:631731/1‑151 (MQ=255)
                aCCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATc                                                                                                                                    >  1:747060/1‑151 (MQ=255)
                 ccGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCg                                                                                                                                   >  2:264055/1‑151 (MQ=255)
                              tgcgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGttt                                                                                                                      >  1:387741/1‑151 (MQ=255)
                               gggttgGGCGCTCTCTGGGGGTGAGGTTTAGGGGTTGTGTTTCACGGTGGGGGAACTGAACCGGGGCGCGCTGGTTGTGGCTCGCGGACTCCGGGGTGTCCCGTTTTCCTTTTTGGGGAAAATGACCGCCGATATCGAAGGGCTTTGtt                                                                                                                       <  1:607336/145‑1 (MQ=255)
                                 gctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGtttt                                                                                                                     >  2:535596/1‑149 (MQ=255)
                                 gctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTAt                                                                                                                   >  1:341851/1‑151 (MQ=255)
                                  ctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTAtt                                                                                                                  >  2:388675/1‑151 (MQ=255)
                                   tgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTc                                                                                                                 >  2:979187/1‑151 (MQ=255)
                                           ctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg                                                                                                          >  1:147922/1‑150 (MQ=255)
                                              tGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGa                                                                                                       >  1:187133/1‑150 (MQ=255)
                                              tGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGa                                                                                                       >  1:1134244/1‑150 (MQ=255)
                                               gggggTGAAGATTGGGTGTTGTTTTTCCCTGGGCCGGACCTGACCGGGGGCGCGTTGGTGGTGGCTCTCGGCCACGGGGGCGTACTTTTTCCCTGTAGCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGa                                                                                                       <  2:620911/146‑1 (MQ=255)
                                               gggggTGAAGATTACGAGTTGTTTTTCACTGTGCCGAAACTGAACGGTGGCGCGCTGGATGTGGCTCCCGGACCCCTGGGCGTCCCTTTTCCCTGTATGGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg                                                                                                          <  1:81881/143‑1 (MQ=255)
                                                       aGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTgat                                                                                                >  2:625194/1‑146 (MQ=255)
                                                       aGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                              >  2:103492/1‑150 (MQ=255)
                                                         aTTAGGGGTTTTGTTTCCCGGGGCCGAATCTGACCGGTGGCCCGTGGGATGTGGCTCTCGGACCCTGGGGCTTCCTTTTTACCTGTATGGGGAAAGGGTCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTa                                                                                               <  1:853603/147‑1 (MQ=255)
                                                           tACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt                                                                                           >  1:1097952/1‑149 (MQ=255)
                                                              gttttGTTTTTCCCTGGGCGCAAACTGAGCCGGGGCGCGCTGGTTGTGCCTCCCGGACGCCGGGGCGTACCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACt                                                                                      <  2:346362/148‑1 (MQ=255)
                                                              gAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACt                                                                                      >  2:1004727/1‑151 (MQ=255)
                                                               aGTTGTGTTTCACTGTGCCGGGACTGAAACGTGGCGCGCTGGATGGGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGCGACGGCGAACCTGTTACATTAGACTg                                                                                     >  2:709301/1‑151 (MQ=255)
                                                                 ttttttttCACTGTGCGGGAACTAACCGGGGGCGCGTTGGTTGGGTCTCTCGGACACGTGGGCTTCCCTTTTACCTGTAGGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg                                                                                    <  2:776940/145‑1 (MQ=255)
                                                                  tgttTTTCCCTGTGCGGGAACTAACCGGTGGGGCGGTGGATGTGGCTCTGGGACACTGGGGCGTACGTTTTACCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                                     <  2:1122215/148‑1 (MQ=255)
                                                                    tttttcctggggCCGGAAGTGATCGGGGGCGCTGGGGATGGGCCCCCCGGACCCCTGGGCGTACTTTTTACCTTTAGCGGACAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaa                                                                                  <  1:787460/138‑1 (MQ=255)
                                                                      tttCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAgg                                                                               >  1:727343/1‑150 (MQ=255)
                                                                        tcctggtgCCGGAACGGACCGTGGGCCGCGTGGATTTTGCCTCGGGACACCTGCGTGTCCTTTTTCCCTGTAGGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAg                                                                                <  2:794778/142‑1 (MQ=255)
                                                                         cACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGata                                                                            >  2:1105878/1‑150 (MQ=255)
                                                                              tgCCGAAACAGACCGGTGGCGCGCTGGTTGTGCCTCTCGGCCACGTGGGTGTCCTTTTTTCCTGTTGCGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATg                                                                          <  1:832892/147‑1 (MQ=255)
                                                                               gcgggAATTGACCCGTGGCGCGGTGGTTGTGGCTCCCGGCCCCGTGGGCGTACTGTTTCCTTTTACGGGCAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAt                                                                     <  2:1134266/148‑1 (MQ=255)
                                                                                ccGGAACTGACCCGGGGCGCGCTGGATGTGTCTCTCGGACACCGGGGCGTCCCTTTTCCCTGTATCGGGCAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAt                                                                     <  1:944979/150‑1 (MQ=255)
                                                                                  ggACCTGACCGGGGGCGCGTTGGATGTGGCTTTGGGACCCCGGGGTGTCTTTTTTCCCTGTGTGGGAAAATGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGatat                                                                           <  2:746078/142‑1 (MQ=255)
                                                                                  ggAACTGAACGGGGGCGCGCTGGTTGTGGCTCTCGGACCCCGGGGGGTCCTGTTTACCTGTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt                                                                  <  2:984609/151‑1 (MQ=255)
                                                                                  ggAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGAATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt                                                                  >  1:137559/1‑151 (MQ=255)
                                                                                  ggAACTGAACCGGGGCGCGCTGGGTGTGGCTCTGGGCCCCTTGGGCGTCCTTTTTTCCTTTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAc                                                                        <  2:878632/145‑1 (MQ=255)
                                                                                    acctgaccgggggCGCGCGGGTTGGGGCTCCCGGACCCCGGGGCGTCCGTTTTCCCTTTATGGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGata                                                                            <  2:372462/128‑1 (MQ=255)
                                                                                    aaCTGACCGGTGGCGCGCTGGGTGTGGCTCTCGGACACCGGGGCGTCCCTTTTCCCTGTATGGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                                <  1:494968/151‑1 (MQ=255)
                                                                                         agccGTGGCGCGTGGGATGTGGCTCTCGGACACGTGGGGGTGCCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACg                                                            <  1:379097/148‑1 (MQ=255)
                                                                                         acccGGGGCGCGTTGGATGTGGCTCTCGGACCCGGGGGCGTCCCTTTTCCCTGTATCGGGAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCa                                                              <  1:293880/146‑1 (MQ=255)
                                                                                         aaCCGAGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGc                                                           >  1:375848/1‑151 (MQ=255)
                                                                                                        gggtgtgGCTCCCGGCCCCCGGGGCTTCCCTTTTCCTTGTATGGGAAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTc                                               <  1:914363/146‑1 (MQ=255)
                                                                                                        ggATGTGGCTCTCGGCCCCCGGGGCGTCCCTTTTCCCTTTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcgc                                             <  1:851153/151‑1 (MQ=255)
                                                                                                          aTGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcgcg                                            >  2:521033/1‑150 (MQ=255)
                                                                                                                  ctcGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGc                                   >  2:607737/1‑151 (MQ=255)
                                                                                                                   cacggaccctggggCGTACTTTTTTCCTGTGCGGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTc                                     <  1:1069164/138‑1 (MQ=255)
                                                                                                                       ggacccTGGGCCGTACCTTTTCCCTTTATCTGCAAAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTc                                               <  2:1156665/132‑1 (MQ=255)
                                                                                                                         ccctggggCGTCCCTTTTCCTTTTATGGGAAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTc                                     <  1:453976/138‑1 (MQ=255)
                                                                                                                         ccccggggCGTCCTTTTTCCCTGTACCGGGAAAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGaa                              <  2:387726/145‑1 (MQ=255)
                                                                                                                                    tACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATc                   >  1:340307/1‑151 (MQ=255)
                                                                                                                                         tttccctgttgggggAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg               <  1:363776/138‑1 (MQ=255)
                                                                                                                                          tccccTGTACGGGGAAAGTGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTg               <  2:236530/146‑1 (MQ=255)
                                                                                                                                            aCCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCAt           >  1:360344/1‑151 (MQ=255)
                                                                                                                                              cTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGCAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGGGGCATc          >  2:1112425/1‑150 (MQ=255)
                                                                                                                                               tGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCt         >  1:803282/1‑150 (MQ=255)
                                                                                                                                                 tATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACt      >  1:1243447/1‑151 (MQ=255)
                                                                                                                                                  aTCTGGGAAATGACGGCCGATATCGAAGGGCTTTGATTTATTCTTGATGGCGAACATGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATTTCGCGAAAAGTCGTGTGAAGATGAGTAATCGGTGGCATCTACtt     >  2:174623/1‑151 (MQ=255)
                                                                                                                                                     gggAAAATGGCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCt  <  1:1105859/151‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
TTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACCGTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGCAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCT  >  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: