Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,270,398 A→G N54S (AAC→AGC)  yjbQ → thiamine phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,270,3980AG98.4% 226.4 / NA 62N54S (AAC→AGC) yjbQthiamine phosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (32/29);  minor base T (0/1);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGAACATCGCTTGTATTGCC  >  NC_000913/4270251‑4270545
                                                                                                                                                   |                                                                                                                                                   
aGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTg                                                                                                                                                  <  1:185800/151‑1 (MQ=255)
 gACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGc                                                                                                                                                 <  2:142081/151‑1 (MQ=255)
   cGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGa                                                                                                                                               <  1:97664/151‑1 (MQ=255)
   cGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGa                                                                                                                                               <  2:161629/151‑1 (MQ=255)
      tCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTAAATGAGAGCTGCGATcc                                                                                                                                            >  1:333395/1‑151 (MQ=255)
       cACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATc                                                                                                                                             <  2:826128/149‑1 (MQ=255)
         cGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCAc                                                                                                                                         >  1:119500/1‑151 (MQ=255)
             tAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGt                                                                                                                                      <  2:447928/150‑1 (MQ=255)
             tAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATAGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCg                                                                                                                                       >  2:121174/1‑149 (MQ=255)
                cgcCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGc                                                                                                                                  >  1:464843/1‑151 (MQ=255)
                 gcCAAATCACGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGcc                                                                                                                                 <  1:150172/151‑1 (MQ=255)
                    aaaTCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACg                                                                                                                              <  1:1047523/151‑1 (MQ=255)
                          cGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATg                                                                                                                         <  1:330240/150‑1 (MQ=255)
                                  tCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGtttt                                                                                                                >  2:867055/1‑151 (MQ=255)
                                  tCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGtttt                                                                                                                >  1:102671/1‑151 (MQ=255)
                                    aTCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGtttttt                                                                                                              >  2:769325/1‑151 (MQ=255)
                                     tCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTtc                                                                                                             >  1:513074/1‑151 (MQ=255)
                                        ggTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCCCCGTACGCCACGACATGGAGCGTTTTTtc                                                                                                             <  2:282055/148‑1 (MQ=255)
                                              ggATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCAcc                                                                                                    <  2:1044087/151‑1 (MQ=255)
                                              ggATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCAcc                                                                                                    >  2:589609/1‑151 (MQ=255)
                                              ggATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCAcc                                                                                                    >  1:876088/1‑151 (MQ=255)
                                                aTGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGt                                                                                                  >  1:850264/1‑151 (MQ=255)
                                                  gAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGtt                                                                                                 <  1:345720/150‑1 (MQ=255)
                                                       tCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGAc                                                                                           >  1:1040856/1‑151 (MQ=255)
                                                         tGAATCAGCTGGCTGATATGCCGCGCGTTAACATCTGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGAc                                                                                           >  1:873848/1‑149 (MQ=255)
                                                         tGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACa                                                                                          >  2:792836/1‑150 (MQ=255)
                                                          gAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACaa                                                                                         <  1:580085/150‑1 (MQ=255)
                                                             tCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACgg                                                                                      >  2:360799/1‑150 (MQ=255)
                                                              cAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGa                                                                                     <  1:412480/150‑1 (MQ=255)
                                                                 ctggctgATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAAt                                                                                  <  2:935581/150‑1 (MQ=255)
                                                                          atGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCAt                                                                         <  1:187606/150‑1 (MQ=255)
                                                                          atGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCAt                                                                         <  2:457936/150‑1 (MQ=255)
                                                                            gCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGa                                                                       <  2:651133/150‑1 (MQ=255)
                                                                              cgcgcgTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTa                                                                    <  1:495130/151‑1 (MQ=255)
                                                                                 gcgTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATg                                                                  >  1:662914/1‑150 (MQ=255)
                                                                                  cgTTAACATCGGCTTACTGCATCTATCGTTGCAACAGGCCCCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCTCCGTACGCCACGACAAGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGa                                                                 <  2:15437/150‑1 (MQ=255)
                                                                                   gTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAgg                                                               >  2:102596/1‑151 (MQ=255)
                                                                                      aaCATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGa                                                             <  2:119496/150‑1 (MQ=255)
                                                                                         aTCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGc                                                           >  2:792983/1‑149 (MQ=255)
                                                                                                tACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGatat                                                  >  2:487086/1‑151 (MQ=255)
                                                                                                    gCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGcc                                               <  1:272920/150‑1 (MQ=255)
                                                                                                    gCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCt                                              <  1:647006/151‑1 (MQ=255)
                                                                                                         tGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTGCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCtt                                             <  1:551830/147‑1 (MQ=255)
                                                                                                           ttGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCat                                        >  2:1026467/1‑150 (MQ=255)
                                                                                                              ctgcAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCat                                        <  2:242591/147‑1 (MQ=255)
                                                                                                              ctgcAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATc                                     >  1:846037/1‑150 (MQ=255)
                                                                                                                gcAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCaa                                   >  1:112182/1‑150 (MQ=255)
                                                                                                                gcAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCaa                                   >  2:489304/1‑150 (MQ=255)
                                                                                                                 cAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAAt                                 >  2:153597/1‑151 (MQ=255)
                                                                                                                 cAACATACCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAAt                                 >  1:907556/1‑151 (MQ=255)
                                                                                                                   aCATACCTCCGCCTCTCTGACGCTTAATGAGATCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCAACGTTCCCGACAACGTAAATTATGAGCATGACTATGAGGGAGCACACGATATGCCTTCTCATATCAAAt                                 <  2:439422/149‑1 (MQ=255)
                                                                                                                       aCCTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCa                            <  1:289478/150‑1 (MQ=255)
                                                                                                                         cTCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATg                         >  2:52933/1‑151 (MQ=255)
                                                                                                                          tCCGCCTCTCTGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGc                        >  2:745141/1‑151 (MQ=255)
                                                                                                                                tctctGACACTTAATGAGAGCTGCGATCCCACTGGACGCCACACCATGGAGCGTTATTTCCTCCGCACCGTTCACGAAAACGGAAATTATGAGAATGACTATGAGAGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGa                   >  1:435106/1‑150 (MQ=255)
                                                                                                                                  tctGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGAAc                 >  1:652840/1‑150 (MQ=255)
                                                                                                                                  tctGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGAAc                 >  2:790990/1‑150 (MQ=255)
                                                                                                                                   ctGACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGaa                  <  2:41322/148‑1 (MQ=255)
                                                                                                                                     gACACTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGAACa                <  2:137197/148‑1 (MQ=255)
                                                                                                                                       cacTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAg                                                                            <  2:582385/86‑1 (MQ=255)
                                                                                                                                       cacTTAATGAGAGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAg                                                                            >  1:582401/1‑86 (MQ=255)
                                                                                                                                               gagaGCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGAACATCGCTTGTATTGc   >  2:85850/1‑151 (MQ=255)
                                                                                                                                                agagCTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGAACATCGCTTGTATTGcc  >  2:345635/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGACGCTCACGCTTAGCGCCAAATCTCGTGGGTTTCATCTGGTAACGGATGAAATTCTGAATCAGCTGGCTGATATGCCGCGCGTTAACATCGGCTTACTGCATCTGTTGCTGCAACATACCTCCGCCTCTCTGACACTTAATGAGAACTGCGATCCCACCGTACGCCACGACATGGAGCGTTTTTTCCTCCGCACCGTTCCCGACAACGGAAATTATGAGCATGACTATGAGGGAGCAGACGATATGCCTTCTCATATCAAATCCTCAATGCTGGGAACATCGCTTGTATTGCC  >  NC_000913/4270251‑4270545

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: