Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 909,422 G→A 76.6% A543V (GCC→GTC)  poxB ← pyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913909,4220GA76.6% 105.0 / 32.9 62A543V (GCC→GTC) poxBpyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base G (5/9);  new base A (25/22);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01

TCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGC  >  NC_000913/909275‑909566
                                                                                                                                                   |                                                                                                                                                
tCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGc                                                                                                                                                 <  2:388866/149‑1 (MQ=255)
        tAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGa                                                                                                                                         <  4:136371/149‑1 (MQ=255)
        tAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGa                                                                                                                                         <  6:141240/149‑1 (MQ=255)
            tcaatcaTGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGttt                                                                                                                                     >  4:64998/1‑149 (MQ=255)
             caatcaTGGCATGTCCTTATTATGCCGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGttt                                                                                                                                     >  5:47458/1‑148 (MQ=255)
              aatcaTGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGa                                                                                                                                   >  4:75859/1‑149 (MQ=255)
                      cATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGt                                                                                                                           >  7:169523/1‑149 (MQ=255)
                       aTGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTg                                                                                                                          <  1:32376/149‑1 (MQ=255)
                       aTGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTg                                                                                                                          <  5:304794/149‑1 (MQ=255)
                       aTGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTg                                                                                                                          <  3:64998/149‑1 (MQ=255)
                         gTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGa                                                                                                                                         <  3:355912/132‑1 (MQ=255)
                         gTCATTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGa                                                                                                                                         >  4:355912/1‑132 (MQ=255)
                          tCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGaa                                                                                                                       <  3:75859/149‑1 (MQ=255)
                           ccTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAAt                                                                                                                      >  4:248778/1‑149 (MQ=255)
                            cTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCAGCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATg                                                                                                                     <  4:358760/149‑1 (MQ=255)
                            cTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCATCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAAt                                                                                                                      <  8:232757/148‑1 (MQ=255)
                              tattatGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGt                                                                                                                                       >  5:306665/1‑129 (MQ=255)
                              tattatGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGt                                                                                                                                       <  6:306665/129‑1 (MQ=255)
                               attatGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAAGCTTTGACCTGTTCGAGTTTGATc                                                                                                                                 <  6:294221/134‑1 (MQ=255)
                               attatGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGc                                                                                                                   <  5:439299/148‑1 (MQ=255)
                               attatGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAAACTTTGACCTGTTCGAGTTTGATc                                                                                                                                 >  5:294221/1‑134 (MQ=255)
                               attatGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCGGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGc                                                                                                                   <  6:47458/148‑1 (MQ=255)
                                  atGACGGGAAATGCCACCCTTTTTACCTTAGCCTGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTaa                                                                                                                <  6:206871/148‑1 (MQ=255)
                                    gACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAActc                                                                                                             >  1:452209/1‑149 (MQ=255)
                                       gggAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGa                                                                                                                                         >  2:59282/1‑118 (MQ=255)
                                       gggAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGa                                                                                                                                         <  1:59282/118‑1 (MQ=255)
                                           aaTGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGCTCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTg                                                                                                      <  8:169523/149‑1 (MQ=255)
                                               ccACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGa                                                                                                  >  2:457586/1‑149 (MQ=255)
                                                cACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGtt                                                                                                                                      <  1:302622/112‑1 (MQ=255)
                                                cACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGtt                                                                                                                                      >  2:302622/1‑112 (MQ=255)
                                                cACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGa                                                                                                  >  6:417472/1‑148 (MQ=255)
                                                 aCCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACtcttc                                                                                                          <  8:275825/139‑1 (MQ=255)
                                                 aCCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACtcttc                                                                                                          >  7:275825/1‑139 (MQ=255)
                                                  cccTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGacca                                                                                               >  8:90300/1‑149 (MQ=255)
                                                   ccTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTTCGCGCAGCATATACAGGCTGAAACCTTTGTCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGaccac                                                                                              >  1:465775/1‑149 (MQ=255)
                                                   ccTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGac                                                                                                 <  6:144908/146‑1 (MQ=255)
                                                   ccTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGac                                                                                                 >  5:144908/1‑146 (MQ=255)
                                                           ccTTAGCCAGTATGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATcca                                                                                      >  4:286849/1‑149 (MQ=255)
                                                                    gtttgtttTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGTTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCg                                                                             >  1:125620/1‑149 (MQ=255)
                                                                     tttgtttTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCg                                                                             >  6:6063/1‑148 (MQ=255)
                                                                       tgtttTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGa                                                                           <  1:457586/148‑1 (MQ=255)
                                                                       tgtttTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGTAATGTCTAATTCTTCTTTTGCGACCACCACATAGACCAACACCGGa                                                                           <  8:444599/148‑1 (MQ=255)
                                                                                      gATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAACCCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAgg                                                           >  5:95615/1‑149 (MQ=255)
                                                                                            ttCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGtt                                                     >  8:417115/1‑149 (MQ=255)
                                                                                               atcaCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGc                                                   <  5:6063/148‑1 (MQ=255)
                                                                                                   ccGCGTCCGCTGATTATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTATTCTTTGGCGACCACCACATCCACCAACACCGGACAGTCGATGGAGAAGGCGCGTTGCAGTGc                                              >  7:201382/1‑149 (MQ=255)
                                                                                                   ccGCGTCCGCTGATGATTGCGCGCAGCATATACAGGTTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGc                                              <  2:125620/149‑1 (MQ=255)
                                                                                                   ccGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGc                                              <  7:417115/149‑1 (MQ=255)
                                                                                                      attCCGCTGATGATTGCTCGCAGCATATACAGGCTGAATCCTTTTACCTGTTAGAGTTTGATCTGCGGTGGAATTGCTAACTCTTCTTTGGCGTCCACCACATAGACCAACACAGGACCGTCGATTGAGATTGCGCGTTGCAGGGATtc                                           <  8:201382/147‑1 (MQ=255)
                                                                                                       gTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACATTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTtca                                          >  6:35643/1‑149 (MQ=255)
                                                                                                               atgatTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACtt                                   >  7:114415/1‑148 (MQ=255)
                                                                                                                      cgcgcAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCTGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGc                            <  2:452209/148‑1 (MQ=255)
                                                                                                                        cgcAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCttt                         >  6:305764/1‑149 (MQ=255)
                                                                                                                         gcagcaTATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCtttt                        >  2:212526/1‑149 (MQ=255)
                                                                                                                                 taCAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTacac                 <  3:248778/148‑1 (MQ=255)
                                                                                                                                   cAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGa              <  7:90300/149‑1 (MQ=255)
                                                                                                                                      gCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTac                   >  8:147971/1‑141 (MQ=255)
                                                                                                                                      gCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTac                   <  7:147971/141‑1 (MQ=255)
                                                                                                                                        tGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGAt             >  2:453426/1‑145 (MQ=255)
                                                                                                                                        tGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGAt             <  1:453426/145‑1 (MQ=255)
                                                                                                                                        tGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCATCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATAccc         >  6:100730/1‑149 (MQ=255)
                                                                                                                                               tttGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGc  >  6:19789/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGC  >  NC_000913/909275‑909566

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