Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,613,731 T→C L611S (TTA→TCA)  hyfR → DNA‑binding transcriptional activator HyfR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,613,7310TC97.8% 1096.1 / 8.7 273L611S (TTA→TCA) hyfRDNA‑binding transcriptional activator HyfR
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base C (136/131);  total (139/134)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.26e-01

Only 100 of 274 total aligned reads displayed.

TCGTGGTAACAGTCTGAATTTACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTTACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCAATGGCATTGTTGCCGGGCCCCGTGGCGCA  >  NC_000913/2613633‑2613831
                                                                                                  |                                                                                                    
   tGGTAACAGTCTGAATTTACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTTACTGc                                                                                                 >  2:3963908/1‑101 (MQ=255)
   tGGTAACAGTCTGAATTTACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGc                                                                                                 <  2:3963461/101‑1 (MQ=255)
    ggTAACAGTCTGAATTTACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCa                                                                                                >  1:3963909/1‑101 (MQ=255)
    ggTAACAGTCTGAATTTACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCa                                                                                                <  1:3963633/101‑1 (MQ=255)
     gTAACAGTCTGAATTTACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCAc                                                                                               <  1:3963463/101‑1 (MQ=255)
       aaCAGTCTGAATTTACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCAccc                                                                                             >  1:3963913/1‑101 (MQ=255)
           gTCTGAATTTACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCcgacg                                                                                         >  2:3963917/1‑101 (MQ=255)
            tCTGAATTTACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCcgacga                                                                                        <  1:3963606/101‑1 (MQ=255)
             cTGAATTTACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCgacgac                                                                                       >  1:3963920/1‑101 (MQ=255)
             cTGAATTTACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCgacgac                                                                                       <  1:3963493/101‑1 (MQ=255)
               gAATTTACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGc                                                                                     <  1:3963514/101‑1 (MQ=255)
               gAATTTACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGc                                                                                     <  2:3963496/101‑1 (MQ=255)
                aaTTTACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCcgacga                                                                                        <  1:3963395/97‑1 (MQ=255)
                aaTTTACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGcc                                                                                    >  2:3963926/1‑101 (MQ=255)
                aaTTTACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGcc                                                                                    >  1:3963924/1‑101 (MQ=255)
                 aTTTACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCa                                                                                   >  2:3963928/1‑101 (MQ=255)
                  tttACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCag                                                                                  >  1:3963931/1‑101 (MQ=255)
                    tACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCagag                                                                                >  1:3963935/1‑101 (MQ=255)
                     aCATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAgaga                                                                               <  2:3963538/101‑1 (MQ=255)
                     aCATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAgaga                                                                               >  2:3963937/1‑101 (MQ=255)
                       aTCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAAt                                                                             <  2:3963568/101‑1 (MQ=255)
                       aTCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAAt                                                                             >  2:3963940/1‑101 (MQ=255)
                        tCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATg                                                                            <  2:3963450/101‑1 (MQ=255)
                        tCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATg                                                                            <  1:3963503/101‑1 (MQ=255)
                         cTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAAt                                                                             <  2:3963502/99‑1 (MQ=255)
                          tAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGAc                                                                          <  2:3963657/101‑1 (MQ=255)
                          tAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGAc                                                                          >  1:3963944/1‑101 (MQ=255)
                            aaTGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACga                                                                        >  2:3963947/1‑101 (MQ=255)
                             aTGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACgaa                                                                       <  1:3963726/101‑1 (MQ=255)
                             aTGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACgaa                                                                       >  1:3963948/1‑101 (MQ=255)
                               gTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACgaaga                                                                     >  1:3963951/1‑101 (MQ=255)
                                tCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACgaagaa                                                                    <  2:3963520/101‑1 (MQ=255)
                                tCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACgaagaa                                                                    <  2:3963555/101‑1 (MQ=255)
                                tCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACgaagaa                                                                    <  2:3963601/101‑1 (MQ=255)
                                tCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACgaagaa                                                                    >  2:3963956/1‑101 (MQ=255)
                                 ccGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACgaagaag                                                                   <  2:3963323/101‑1 (MQ=255)
                                 ccGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACgaagaag                                                                   >  2:3963959/1‑101 (MQ=255)
                                  cGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACgaagaaga                                                                  <  1:3963592/101‑1 (MQ=255)
                                  cGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACgaagaaga                                                                  >  1:3963960/1‑101 (MQ=255)
                                   gACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACgaagaagaa                                                                 <  1:3963685/101‑1 (MQ=255)
                                     cAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACg                                                               <  1:3963831/101‑1 (MQ=255)
                                     cAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACg                                                               >  2:3963965/1‑101 (MQ=255)
                                     cAAAGCCGGTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACg                                                               >  1:3963966/1‑101 (MQ=255)
                                      aaaGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGt                                                              <  1:3963650/101‑1 (MQ=255)
                                         gCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAAc                                                                >  1:3963968/1‑96 (MQ=255)
                                         gCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAg                                                           <  1:3963583/101‑1 (MQ=255)
                                          ccGTTTACTGCCGACGCTAAATGAAGATTAAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAgc                                                          <  2:3963695/101‑1 (MQ=255)
                                           cGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAgcg                                                         >  1:3963972/1‑101 (MQ=255)
                                            gTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAgcgc                                                        >  2:3963976/1‑101 (MQ=255)
                                            gTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAgcgc                                                        <  2:3963709/101‑1 (MQ=255)
                                            gTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAgcgc                                                        <  2:3963559/101‑1 (MQ=255)
                                            gTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAgcgc                                                        <  1:3963622/101‑1 (MQ=255)
                                             tttACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACgaagaaga                                                                  >  2:3963978/1‑90 (MQ=255)
                                              ttACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCAt                                                      <  1:3963678/101‑1 (MQ=255)
                                                 cTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGt                                                   >  2:3963981/1‑101 (MQ=255)
                                                  tGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGtt                                                  >  1:3963983/1‑101 (MQ=255)
                                                  tGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGtt                                                  <  2:3963749/101‑1 (MQ=255)
                                                  tGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGtt                                                  >  2:3963985/1‑101 (MQ=255)
                                                    ccGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCa                                                <  1:3963546/101‑1 (MQ=255)
                                                     cGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAg                                               <  1:3963440/101‑1 (MQ=255)
                                                       aCGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGt                                             >  1:3963991/1‑101 (MQ=255)
                                                       aCGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGt                                             >  2:3963990/1‑101 (MQ=255)
                                                        cGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTa                                            >  2:3963993/1‑101 (MQ=255)
                                                        cGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTa                                            <  2:3963693/101‑1 (MQ=255)
                                                          cTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTAtt                                          >  2:3963999/1‑101 (MQ=255)
                                                          cTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTAtt                                          <  2:3963575/101‑1 (MQ=255)
                                                          cTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTAtt                                          >  1:3963997/1‑101 (MQ=255)
                                                           tAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTTACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTg                                         >  2:3964005/1‑101 (MQ=255)
                                                           tAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTa                                            <  1:3963528/98‑1 (MQ=255)
                                                           tAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTg                                         <  2:3963394/101‑1 (MQ=255)
                                                           tAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTg                                         >  1:3964002/1‑101 (MQ=255)
                                                            aaaTGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGc                                        <  2:3963672/101‑1 (MQ=255)
                                                             aaTGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCg                                       >  1:3964010/1‑101 (MQ=255)
                                                             aaTGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCg                                       <  2:3963584/101‑1 (MQ=255)
                                                             aaTGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCg                                       >  2:3964009/1‑101 (MQ=255)
                                                              aTGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCga                                      >  2:3964012/1‑101 (MQ=255)
                                                                 aaGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGaa                                   >  1:3964016/1‑101 (MQ=255)
                                                                 aaGATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGaa                                   >  2:3964018/1‑101 (MQ=255)
                                                                   gATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAAc                                 >  1:3964021/1‑101 (MQ=255)
                                                                   gATTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAAc                                 <  1:3963094/101‑1 (MQ=255)
                                                                    aTTCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAAcc                                >  2:3964024/1‑101 (MQ=255)
                                                                     ttCGGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCa                               <  2:3963615/101‑1 (MQ=255)
                                                                     ttCAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCa                               <  1:3963561/101‑1 (MQ=255)
                                                                       cAGCGCTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCAAt                             >  2:3964027/1‑101 (MQ=255)
                                                                         gcgcTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCa                                                       <  2:3963473/73‑1 (MQ=255)
                                                                          cgcTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCAATGGc                          <  1:3963459/101‑1 (MQ=255)
                                                                            cTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCAATGGCAt                        <  1:3963664/101‑1 (MQ=255)
                                                                            cTTCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCAATGGCAt                        <  2:3963578/101‑1 (MQ=255)
                                                                             ttCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCAATGGCAtt                       >  1:3964029/1‑101 (MQ=255)
                                                                              tCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCAATGGCAttg                      >  2:3964031/1‑101 (MQ=255)
                                                                              tCGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCAATGGCAttg                      <  2:3963594/101‑1 (MQ=255)
                                                                               cGCAGTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCAATGGCAttgt                     <  1:3963614/101‑1 (MQ=255)
                                                                                   gTTCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCAATGGCATTGTTGcc                 >  1:3964035/1‑101 (MQ=255)
                                                                                    ttCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCAATGGCATTGTTGCCg                >  2:3964037/1‑101 (MQ=255)
                                                                                    ttCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCAATGGCATTGTTGCCg                <  2:3963839/101‑1 (MQ=255)
                                                                                     tCAATGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCAATGGCATTGTTGCCgg               <  2:3963635/101‑1 (MQ=255)
                                                                                       aaTGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCAATGGCATTGTTGCCGGGc             <  1:3963765/101‑1 (MQ=255)
                                                                                        aTGGCGCAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCAATGGCATTGTTGCCGGGcc            >  1:3964041/1‑101 (MQ=255)
                                                                                           gcgcAGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCAATGGCATTGTTGCCGGGCCCCg         >  1:3964044/1‑101 (MQ=255)
                                                                                               aGTCACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCAATGGCATTGTTGCCGGGCCCCGTGgc     >  2:3964046/1‑101 (MQ=255)
                                                                                                  |                                                                                                    
TCGTGGTAACAGTCTGAATTTACATCTAAATGTCCGACAAAGCCGTTTACTGCCGACGCTAAATGAAGATTCAGCGCTTCGCAGTTCAATGGCGCAGTTACTGCACCCGACGACGCCAGAGAATGACGAAGAAGAACGTCAGCGCATTGTTCAGGTATTGCGAGAAACCAATGGCATTGTTGCCGGGCCCCGTGGCGCA  >  NC_000913/2613633‑2613831

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