Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,631,540 C→T K141K (AAG→AAA)  SAUSA300_RS08065 ← integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,631,5400CT99.2% 344.2 / ‑5.8 122K141K (AAG→AAA) SAUSA300_RS08065integrase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (3/118);  total (3/119)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 122 total aligned reads displayed.

TCTTCACGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTCTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATACCCTGTAT  >  NC_007793/1631447‑1631637
                                                                                             |                                                                                                 
tcttcACGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCCCGCTTTTCAGAAGACTTATTATATTTTTTGCCt                                                                                             <  1:5345881/100‑1 (MQ=255)
tcttcACGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCt                                                                                             <  1:901210/100‑1 (MQ=255)
tcttcACGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCt                                                                                             <  1:187410/100‑1 (MQ=255)
tcttcACGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCt                                                                                             <  1:2552592/100‑1 (MQ=255)
 cttcACGTTTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTg                                                                                            <  1:3555825/100‑1 (MQ=255)
 cttcACGCTTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTg                                                                                            <  1:412746/100‑1 (MQ=255)
 cttcACGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTg                                                                                            <  1:2203758/100‑1 (MQ=255)
  ttcACGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTgg                                                                                           <  1:430871/100‑2 (MQ=255)
  ttcACGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGAACCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTga                                                                                           <  1:3916688/100‑2 (MQ=255)
   tccgatctattaATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCt                                                                                          <  1:108222/93‑1 (MQ=255)
   tcACGGTTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGcg                                                                                          <  1:2033966/100‑2 (MQ=255)
   tcACGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTgag                                                                                          <  1:4653133/100‑3 (MQ=255)
    cACGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCtg                                                                                         <  1:3355411/100‑2 (MQ=255)
     aCGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTg                                                                                        <  1:2886230/100‑1 (MQ=255)
     aCGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTg                                                                                        <  1:3322738/100‑1 (MQ=255)
     aCGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTg                                                                                        <  1:5348603/100‑1 (MQ=255)
     aCGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTg                                                                                        <  1:5778070/100‑1 (MQ=255)
     aCGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTg                                                                                        <  1:3807177/100‑1 (MQ=255)
     aCGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGCCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTg                                                                                        <  1:4250051/100‑1 (MQ=255)
      cgatctgttaTACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTgg                                                                                       <  1:2119346/91‑2 (MQ=255)
      cGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTgg                                                                                       <  1:5051083/100‑2 (MQ=255)
      cGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTgc                                                                                       <  1:1214445/100‑2 (MQ=255)
      cGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGt                                                                                       <  1:4734095/100‑1 (MQ=255)
       gATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGtc                                                                                      <  1:3247213/100‑1 (MQ=255)
       gATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAAACTTATTATATTTTTTGCCTGCTTGtt                                                                                      <  1:4371602/100‑2 (MQ=255)
        gttattaATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGtct                                                                                     <  1:5755553/99‑1 (MQ=255)
        attattAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGtct                                                                                     <  1:4220167/100‑1 (MQ=255)
         ttattaATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGtctc                                                                                    <  1:4467393/100‑1 (MQ=255)
         ttattaATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGtctc                                                                                    <  1:1464833/100‑1 (MQ=255)
          tattaATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGtcatt                                                                                   <  1:2230909/100‑4 (MQ=255)
           attaATTCTTTTTGAACCGATTCTTTATAACACGAGACCCTACGTCTAATTTGCCCCGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCat                                                                                  <  1:1439488/100‑1 (MQ=255)
           attaATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCat                                                                                  <  1:4414101/100‑1 (MQ=255)
           atataTACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCat                                                                                  <  1:4372765/96‑1 (MQ=255)
             taATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCatat                                                                                <  1:3238405/100‑1 (MQ=255)
              aaTACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATc                                                                               <  1:5154396/100‑1 (MQ=255)
              aaTACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATACTTTTTGCCTGCTTGTCTCATATc                                                                               <  1:3207730/100‑1 (MQ=255)
                 aCTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCAca                                                                            <  1:1361269/100‑2 (MQ=255)
                 aCTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACg                                                                            <  1:1092194/100‑1 (MQ=255)
                   tttttGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCt                                                                          <  1:1047915/100‑1 (MQ=255)
                    ttttGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCAcacca                                                                         <  1:4961604/100‑5 (MQ=255)
                    ttttGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCtt                                                                         <  1:4652065/100‑1 (MQ=255)
                     tttGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTa                                                                        <  1:747685/100‑1 (MQ=255)
                      ttGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTat                                                                       <  1:626385/100‑2 (MQ=255)
                      ttGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTaa                                                                       <  1:4417136/100‑1 (MQ=255)
                       tGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAAc                                                                      <  1:5337213/100‑1 (MQ=255)
                       tGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAAc                                                                      <  1:2294202/100‑1 (MQ=255)
                       tGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAAc                                                                      <  1:3891771/100‑1 (MQ=255)
                          gCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATc                                                                   <  1:664096/100‑1 (MQ=255)
                           ccGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACAtgc                                                                  <  1:730973/100‑3 (MQ=255)
                           ccGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCa                                                                  <  1:1784994/100‑1 (MQ=255)
                            cGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCag                                                                 <  1:5676599/100‑2 (MQ=255)
                             gATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCaag                                                                <  1:835346/100‑2 (MQ=255)
                                aCTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACtgg                                                             <  1:2169554/100‑3 (MQ=255)
                                    tataACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTAcc                                                         <  1:3742714/100‑1 (MQ=255)
                                    tataACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTAcc                                                         <  1:92678/100‑1 (MQ=255)
                                    tataACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTAcc                                                         <  1:4792242/100‑1 (MQ=255)
                                     ataACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCa                                                        <  1:2742824/100‑1 (MQ=255)
                                     ataACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCa                                                        <  1:5790644/100‑1 (MQ=255)
                                     ataACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCa                                                        <  1:703185/100‑1 (MQ=255)
                                     ataACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCa                                                        <  1:4970022/100‑1 (MQ=255)
                                      taACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCac                                                       <  1:3170023/100‑2 (MQ=255)
                                       aaCACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCatt                                                      <  1:5223735/100‑3 (MQ=255)
                                       aaCACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCaag                                                      <  1:3100708/100‑2 (MQ=255)
                                            gagaCCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCagt                                                 <  1:2714629/100‑3 (MQ=255)
                                            gagaCCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCAtt                                                 <  1:2093250/100‑1 (MQ=255)
                                            gagaCCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCAtg                                                 <  1:4451653/100‑2 (MQ=255)
                                            gagaCCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCAtc                                                 <  1:3329838/100‑2 (MQ=255)
                                             agaCCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCAtgc                                                <  1:3992061/100‑3 (MQ=255)
                                             agaCCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTg                                                <  1:3345565/100‑1 (MQ=255)
                                              gaCCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGc                                               <  1:3150256/100‑1 (MQ=255)
                                              gaCCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGc                                               <  1:2441647/100‑1 (MQ=255)
                                               aCCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCg                                              <  1:871028/100‑1 (MQ=255)
                                               aCCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCg                                              <  1:3782973/100‑1 (MQ=255)
                                               aCCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCg                                              <  1:490418/100‑1 (MQ=255)
                                               aCCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCg                                              <  1:1175322/100‑1 (MQ=255)
                                                cccTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGcca                                             <  1:1271903/100‑3 (MQ=255)
                                                cccTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCgc                                             <  1:2387499/100‑2 (MQ=255)
                                                cccTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGa                                             <  1:3822178/100‑1 (MQ=255)
                                                cccTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGa                                             <  1:1008524/100‑1 (MQ=255)
                                                cccTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGa                                             <  1:3370902/100‑1 (MQ=255)
                                                cccTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTCAACATCAACTTTTACCAATCCATTGCGa                                             <  1:233236/100‑1 (MQ=255)
                                                atcTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGa                                             <  1:1967816/98‑1 (MQ=255)
                                                 ccTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGaa                                            <  1:4652388/100‑1 (MQ=255)
                                                 ccTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGaa                                            >  1:2581116/1‑100 (MQ=255)
                                                  ctgcGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAat                                           <  1:1553802/97‑1 (MQ=255)
                                                  cTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGaaa                                           <  1:3258170/100‑2 (MQ=255)
                                                  cTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAat                                           <  1:4356527/100‑1 (MQ=255)
                                                     cGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAatata                                        <  1:4650465/100‑1 (MQ=255)
                                                                ccACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATAtag                             >  1:2138704/1‑100 (MQ=255)
                                                                ccACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATAtag                             >  1:4425352/1‑100 (MQ=255)
                                                                 cACGCTTTTCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATAtagt                            <  1:5127942/100‑1 (MQ=255)
                                                                      ttttCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGAt                       <  1:5493480/100‑1 (MQ=255)
                                                                       tttCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAATTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATa                      <  1:2755050/100‑1 (MQ=255)
                                                                         tCAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATatc                    <  1:4037184/100‑3 (MQ=255)
                                                                          cAGAAGACTTATTATATTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACt                   <  1:1274442/100‑1 (MQ=255)
                                                                                       atatTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGccac      <  1:2439795/100‑3 (MQ=255)
                                                                                       atatTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCt      <  1:2072346/100‑1 (MQ=255)
                                                                                         atTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGt    <  1:4072528/100‑1 (MQ=255)
                                                                                         atTTTTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGt    <  1:3750953/100‑1 (MQ=255)
                                                                                           atCTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATGCCCTGTag  <  1:5206481/99‑2 (MQ=255)
                                                                                             |                                                                                                 
TCTTCACGATTATTAATACTTTTTGAGCCGATACTTTATAACACGAGACCCTACGTCTAATTTGCCACGCTTTTCAGAAGACTTATTATATTTCTTGCCTGCTTGTCTCATATCACGCTTAACATCAACTTTTACCAATCCATTGCGAATATATCTATATATAGTAGGATAACTTATATAATACCCTGTAT  >  NC_007793/1631447‑1631637

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