Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,138,732 T→C H277H (CAT→CAC yegH → inner membrane protein YegH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,138,7320TC90.6% 1112.9 / 85.0 330H277H (CAT→CACyegHinner membrane protein YegH
Reads supporting (aligned to +/- strand):  ref base T (15/16);  new base C (123/176);  total (138/192)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

Only 100 of 330 total aligned reads displayed.

CTGCGCCAGCGGACTACCGAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCATGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCGTCAGCAGCATTATGACGTCGCGCC  >  NC_000913/2138634‑2138829
                                                                                                  |                                                                                                 
cTGCGCCAGCGGACTACCGAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACgg                                                                                                 <  2:3851449/101‑1 (MQ=255)
cTGCGCCAGCGGACTACCGAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACgg                                                                                                 >  1:3851686/1‑101 (MQ=255)
 tGCGCCAGCGGACTACCGAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGt                                                                                                <  2:3850996/101‑1 (MQ=255)
 tGCGCCAGCGGACTACCGAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGt                                                                                                >  2:3851687/1‑101 (MQ=255)
  gcgcCAGCGGACTACCGAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTa                                                                                               <  2:3851206/101‑1 (MQ=255)
      cAGCGGACTACCGAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACca                                                                                           >  2:3851691/1‑101 (MQ=255)
      cAGCGGACTACCGAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACca                                                                                           <  2:3851161/101‑1 (MQ=255)
       aGCGGACTACCGAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACcag                                                                                          <  2:3851307/101‑1 (MQ=255)
        gCGGACTACCGAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCATGGTAACcagc                                                                                         <  2:3851078/101‑1 (MQ=255)
          ggACTACCGAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAACCCGCGTCCATGCTGGTGGATCACgg                                                                                                 <  1:3851160/91‑1 (MQ=255)
           gACTACCGAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGa                                                                                      >  1:3851697/1‑101 (MQ=255)
            aCTACCGAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGAt                                                                                     >  2:3851699/1‑101 (MQ=255)
             cTACCGAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAGGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCATGGTAACCAGCAGATc                                                                                    <  1:3851173/101‑1 (MQ=255)
             cTACCGAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACcagc                                                                                         <  2:3851184/96‑1 (MQ=255)
              tACCGAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCt                                                                                   >  1:3851702/1‑101 (MQ=255)
                 cGAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTa                                                                                <  1:3850755/101‑1 (MQ=255)
                 cGAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTa                                                                                >  2:3851707/1‑101 (MQ=255)
                  gAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTaa                                                                               <  1:3851119/101‑1 (MQ=255)
                  gAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTaa                                                                               >  1:3851708/1‑101 (MQ=255)
                    aGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATc                                                                             >  1:3851712/1‑101 (MQ=255)
                      cGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCg                                                                           <  2:3851038/101‑1 (MQ=255)
                      cGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCg                                                                           <  1:3851346/101‑1 (MQ=255)
                      cGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCg                                                                           >  2:3851715/1‑101 (MQ=255)
                       ggTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCg                                                                           <  1:3851330/100‑1 (MQ=255)
                         tAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGAt                                                                                     >  2:3851720/1‑88 (MQ=255)
                          aaTGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAgg                                                                       >  1:3851721/1‑101 (MQ=255)
                           aTGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGa                                                                      >  2:3851725/1‑101 (MQ=255)
                           aTGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGa                                                                      <  1:3851137/101‑1 (MQ=255)
                              cgcCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAAcg                                                                   <  2:3850864/101‑1 (MQ=255)
                              cgcCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAAcg                                                                   >  2:3851730/1‑101 (MQ=255)
                               gcCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAAcgg                                                                  >  2:3851733/1‑101 (MQ=255)
                                ccTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCATGGTAACCAGCAGATCTTTAATCCGCAGGAAcggc                                                                 <  2:3851312/101‑1 (MQ=255)
                                ccTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAAcggc                                                                 >  1:3851735/1‑101 (MQ=255)
                                  tGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAAcggcgg                                                               >  2:3851739/1‑101 (MQ=255)
                                    cTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGgat                                                             <  2:3851113/101‑1 (MQ=255)
                                    cTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGgat                                                             >  1:3851743/1‑101 (MQ=255)
                                      tAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGgatga                                                           <  2:3851285/101‑1 (MQ=255)
                                        gCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGAtt                                                         <  2:3851129/101‑1 (MQ=255)
                                          gggCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGa                                                       <  2:3851067/101‑1 (MQ=255)
                                          gggCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGa                                                       <  1:3851143/101‑1 (MQ=255)
                                            gCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGc                                                     >  1:3851748/1‑101 (MQ=255)
                                            gCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGc                                                     >  2:3851749/1‑101 (MQ=255)
                                             cAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAAc                                                                    >  1:3851753/1‑85 (MQ=255)
                                             cAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCg                                                    <  2:3851481/101‑1 (MQ=255)
                                             cAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCg                                                    <  2:3851267/101‑1 (MQ=255)
                                              aaaaaGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGAtt                                                         <  1:3851262/95‑1 (MQ=255)
                                              aaaaaGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCgg                                                   >  1:3851756/1‑101 (MQ=255)
                                               aaaaGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCggg                                                  <  1:3851084/101‑1 (MQ=255)
                                               aaaaGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCggg                                                  <  1:3851202/101‑1 (MQ=255)
                                               aaaaGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCggg                                                  <  2:3851099/101‑1 (MQ=255)
                                                aaaGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCggg                                                  >  2:3851759/1‑100 (MQ=255)
                                                aaaGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGt                                                 <  2:3851217/101‑1 (MQ=255)
                                                aaaGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGt                                                 <  1:3851232/101‑1 (MQ=255)
                                                  agaagaTGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTAc                                               <  2:3851194/101‑1 (MQ=255)
                                                   gaagaTGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACt                                              <  2:3851245/101‑1 (MQ=255)
                                                    aagaTGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACtt                                             >  2:3851766/1‑101 (MQ=255)
                                                    aagaTGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACtt                                             >  1:3851765/1‑101 (MQ=255)
                                                     agaTGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTa                                            <  1:3851094/101‑1 (MQ=255)
                                                     agaTGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTa                                            <  2:3851235/101‑1 (MQ=255)
                                                      gaTGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTaa                                           <  1:3851320/101‑1 (MQ=255)
                                                       aTGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAAt                                          >  2:3851772/1‑101 (MQ=255)
                                                       aTGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAAt                                          <  1:3851246/101‑1 (MQ=255)
                                                       aTGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAAt                                          >  1:3851771/1‑101 (MQ=255)
                                                        tGCGGAACTGGACGCCTAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAAcggcg                                                                <  2:3851331/78‑1 (MQ=255)
                                                        tGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATc                                         <  1:3851008/101‑1 (MQ=255)
                                                          cGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCtt                                       <  2:3851300/101‑1 (MQ=255)
                                                           ggAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTa                                      >  1:3851776/1‑101 (MQ=255)
                                                            gAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTaa                                     >  2:3851778/1‑101 (MQ=255)
                                                              aCTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAAcc                                   <  1:3850960/101‑1 (MQ=255)
                                                              aCTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAAcc                                   >  2:3851781/1‑101 (MQ=255)
                                                               cTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCa                                  <  2:3851251/101‑1 (MQ=255)
                                                               cTGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCa                                  <  2:3851417/101‑1 (MQ=255)
                                                                tGGACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAg                                 <  1:3851178/101‑1 (MQ=255)
                                                                 ggACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGc                                <  1:3851389/101‑1 (MQ=255)
                                                                 ggACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGc                                <  1:3851375/101‑1 (MQ=255)
                                                                 ggACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGc                                <  2:3851361/101‑1 (MQ=255)
                                                                 ggACGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGc                                >  2:3851785/1‑101 (MQ=255)
                                                                   aCGCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGt                              <  1:3851216/101‑1 (MQ=255)
                                                                     gCCGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTAc                            >  2:3851789/1‑101 (MQ=255)
                                                                      ccGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTAcc                           >  2:3851793/1‑101 (MQ=255)
                                                                       cGAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCg                          <  1:3851236/101‑1 (MQ=255)
                                                                        gAAACCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCGt                         >  1:3851796/1‑101 (MQ=255)
                                                                         aaaCCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCGTc                        <  1:3851277/101‑1 (MQ=255)
                                                                          aaCCGCGTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCGTca                       >  2:3851800/1‑101 (MQ=255)
                                                                              gcgTCCATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAg                                 >  2:3851803/1‑87 (MQ=255)
                                                                                  ccATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCGTCAGCAGCAtt               <  1:3851306/101‑1 (MQ=255)
                                                                                  ccATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCGTCAGCAGCAtt               <  1:3851515/101‑1 (MQ=255)
                                                                                  ccATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCGTCAGCAGCAtt               >  1:3851806/1‑101 (MQ=255)
                                                                                   cATGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCGTCAGCAGCATTa              <  1:3851457/101‑1 (MQ=255)
                                                                                    aTGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCGTCAGCAGCATTAt             >  1:3851810/1‑101 (MQ=255)
                                                                                     tGCTGGTGGATCATGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCGTCAGCAGCATTATg            <  1:3850914/101‑1 (MQ=255)
                                                                                     tGCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCGTCAGCAGCATTATg            <  1:3851052/101‑1 (MQ=255)
                                                                                      gCTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGGACTTa                                            >  2:3851812/1‑68 (MQ=255)
                                                                                       cTGGTGGATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCGTCAGCAGCATTATGAc          <  2:3851547/101‑1 (MQ=255)
                                                                                        tggtggATCACGGTAACCAGCAGATCTTTAATCCGCACGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCGTCAGCAGCATTATGACg         >  1:3851814/1‑101 (MQ=255)
                                                                                         ggtggATCATGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCGTCAGCAGCATTATGACGt        <  2:3851368/101‑1 (MQ=255)
                                                                                          gtggATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCGTCAGCAGCATTATGACGTc       <  1:3851401/101‑1 (MQ=255)
                                                                                             gATCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCGTCAGCAGCATTATGACGTcgcg    <  2:3851391/101‑1 (MQ=255)
                                                                                              aTCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCGTCAGCAGCATTATGACGTCgcgc   >  1:3851820/1‑101 (MQ=255)
                                                                                               tCACGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCGTCAGCAGCATTATGACGTCGCGcc  >  1:3851824/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
CTGCGCCAGCGGACTACCGAAGCGGTAATGCGCCTGCTTAGCGGGCAAAAAGAAGATGCGGAACTGGACGCCGAAACCGCGTCCATGCTGGTGGATCATGGTAACCAGCAGATCTTTAATCCGCAGGAACGGCGGATGATTGAGCGGGTACTTAATCTTAACCAGCGTACCGTCAGCAGCATTATGACGTCGCGCC  >  NC_000913/2138634‑2138829

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: