Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 909,400 C→A 26.6% M550I (ATG→ATT 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,4000CA26.6% 139.0 / 44.2 74M550I (ATG→ATTpoxBpyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (23/31);  new base A (11/9);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

CGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTT  >  NC_000913/909252‑909543
                                                                                                                                                    |                                                                                                                                               
cgcgGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCgcagc                                                                                                                                                 <  5:242497/149‑1 (MQ=255)
    gCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCatat                                                                                                                                             <  6:395325/149‑1 (MQ=255)
         gTCGTGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAgg                                                                                                                                        <  1:376330/149‑1 (MQ=255)
          tCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGc                                                                                                                                       <  7:409740/149‑1 (MQ=255)
                    gTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCttt                                                                                                                             <  2:274555/149‑1 (MQ=255)
                    gTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCttt                                                                                                                             <  6:73257/149‑1 (MQ=255)
                     tATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTg                                                                                                                            <  5:398173/149‑1 (MQ=255)
                       tCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAata                                                                                                                                              >  2:88512/1‑129 (MQ=255)
                       tCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAata                                                                                                                                              <  1:88512/129‑1 (MQ=255)
                          cTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTg                                                                                                                       >  2:115573/1‑149 (MQ=255)
                             cGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGtt                                                                                                                     <  8:239692/148‑1 (MQ=255)
                              gTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGGCCGCTGATTATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTc                                                                                                                    <  4:349874/148‑1 (MQ=255)
                               tAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTc                                                                                                                    >  2:231971/1‑147 (MQ=255)
                                  atcaatcaTGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGtt                                                                                                               >  2:35276/1‑149 (MQ=255)
                                  atcaatcaTGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGtt                                                                                                               <  1:231971/149‑1 (MQ=255)
                                  atcaatcaTGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGtt                                                                                                               >  3:113368/1‑149 (MQ=255)
                                          tGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCt                                                                                                         >  4:59019/1‑147 (MQ=255)
                                            gCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCGTGTTCGAGTTTGATCTGCg                                                                                                      <  2:14859/148‑1 (MQ=255)
                                             cATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGt                                                                                                    <  8:116580/149‑1 (MQ=255)
                                                gTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTgg                                                                                                  >  3:133051/1‑148 (MQ=255)
                                                 tCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGaa                                                                                                <  8:86644/149‑1 (MQ=255)
                                                   cTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAAt                                                                                               <  4:113368/148‑1 (MQ=255)
                                                   cTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATg                                                                                              >  7:429617/1‑149 (MQ=255)
                                                    ttattaTGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATAACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATgg                                                                                             >  7:227984/1‑149 (MQ=255)
                                                    ttattaTGACGGGAAATGCCACCCTTTTTACCCTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATgg                                                                                             >  3:114559/1‑149 (MQ=255)
                                                      attatGACGGGAAATGCCACCCTTTTTACCTTATCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAACCCTTTGGCCTGTTCGAGTTTGATCTTCGGTTGAATGGc                                                                                            >  2:363008/1‑148 (MQ=255)
                                                         atGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGCCCTGTTCGAGTTTGATCTGCGGTGGAAt                                                                                               >  5:97377/1‑142 (MQ=255)
                                                         atGACGGGAAATGCCAACCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAAt                                                                                               <  6:97377/142‑1 (MQ=255)
                                                          tGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAAc                                                                                        <  8:143087/148‑1 (MQ=255)
                                                          tGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTg                                                                                                   >  7:461100/1‑137 (MQ=255)
                                                          tGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTg                                                                                                   <  8:461100/137‑1 (MQ=255)
                                                          tGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTACGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTCGCCTGTTCGAGTATGATCTGCGGTGGAATGGCAAAct                                                                                       >  6:39896/1‑149 (MQ=255)
                                                               ggAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACtcttct                                                                                  <  8:429617/149‑1 (MQ=255)
                                                                  aaTGCCACCCTTTTTACCTTAGCGAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTTATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTg                                                                               <  2:287916/149‑1 (MQ=255)
                                                                  aaTGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTg                                                                               <  4:133051/149‑1 (MQ=255)
                                                                  aaTGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTg                                                                               <  8:98623/149‑1 (MQ=255)
                                                                        aCCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGc                                                                                            <  4:78129/130‑1 (MQ=255)
                                                                        aCCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGc                                                                                            >  3:78129/1‑130 (MQ=255)
                                                                         cccTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGacca                                                                        >  7:86648/1‑149 (MQ=255)
                                                                            tttttACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGa                                                                                                                  >  5:304946/1‑104 (MQ=255)
                                                                            tttttACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGa                                                                                                                  <  6:304946/104‑1 (MQ=255)
                                                                            tttttACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGAccacca                                                                     >  6:76273/1‑149 (MQ=255)
                                                                               ttACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACcaca                                                                   <  2:308821/148‑1 (MQ=255)
                                                                                     tAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGacc                                                                         <  5:151098/136‑1 (MQ=255)
                                                                                     tAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGacc                                                                         >  6:151098/1‑136 (MQ=255)
                                                                                      aGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATccacca                                                            <  8:86648/148‑1 (MQ=255)
                                                                                      aGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATccacca                                                            <  3:59019/148‑1 (MQ=255)
                                                                                               gtttTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACATTCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGAc                                                   >  4:78895/1‑148 (MQ=255)
                                                                                               gtttTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGAc                                                   <  1:363008/148‑1 (MQ=255)
                                                                                                ttttCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCt                                                                                                                               <  4:104605/71‑1 (MQ=255)
                                                                                                ttttCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCt                                                                                                                               >  3:104605/1‑71 (MQ=255)
                                                                                                    caccAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGa                                             >  7:436256/3‑149 (MQ=255)
                                                                                                    caccAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCCCAGCATATACAGGCTGCAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGCACCCTCGa                                             >  1:365559/3‑149 (MQ=255)
                                                                                                      acatttCGATCACTTCATCTCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGTCCACCACATCCACCAACACCGGACCGTCGATg                                           <  8:227984/145‑1 (MQ=255)
                                                                                                                     catcaCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGGTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTg                             >  5:139968/1‑148 (MQ=255)
                                                                                                                     catcaCCGCGTCCGCTGATCATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGGTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTAGGCGACCACCACATCCACCAACACCGGACCGACGATGGAGAAGGCGCGTGGc                            >  6:451890/1‑149 (MQ=255)
                                                                                                                      atcaCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCa                           <  5:129292/149‑1 (MQ=255)
                                                                                                                      atcaCCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATg                                           <  1:138451/133‑1 (MQ=255)
                                                                                                                      atcaCCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATg                                           >  2:138451/1‑133 (MQ=255)
                                                                                                                      atcaCCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTg                             >  1:322980/1‑147 (MQ=255)
                                                                                                                       tcaCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAg                          <  1:115573/149‑1 (MQ=255)
                                                                                                                         aCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTc                                               >  7:445803/1‑126 (MQ=255)
                                                                                                                         aCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTc                                               <  8:445803/126‑1 (MQ=255)
                                                                                                                             cgTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTtc                    <  5:76273/149‑1 (MQ=255)
                                                                                                                              gTCGGCTGATGATTGCGCGCAGCATATACAGTCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACCCCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTtca                   <  5:39896/149‑1 (MQ=255)
                                                                                                                              gTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTTACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTtca                   <  4:114559/149‑1 (MQ=255)
                                                                                                                              gTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTg                                                                               <  6:130518/89‑1 (MQ=255)
                                                                                                                              gTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTg                                                                               >  5:130518/1‑89 (MQ=255)
                                                                                                                               tCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAg                          <  2:447171/141‑1 (MQ=255)
                                                                                                                               tCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAg                          >  1:447171/1‑141 (MQ=255)
                                                                                                                                    tgatgaTTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACAATTCTTTGGCGACAACCAAATACACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACt             >  5:230348/1‑149 (MQ=255)
                                                                                                                                       tgatTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCa          <  2:322980/149‑1 (MQ=255)
                                                                                                                                        gatTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACtt            >  4:314651/1‑146 (MQ=255)
                                                                                                                                              gcgcAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCtt   >  6:405681/1‑149 (MQ=255)
                                                                                                                                               cgcAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCttt  >  6:434007/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CGCGGCTTGGTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTT  >  NC_000913/909252‑909543

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: