Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,631,762 Δ1 bp coding (201/504 nt) SAUSA300_RS08065 ← integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,631,7620C.98.5% 504.8 / ‑0.4 131coding (201/504 nt)SAUSA300_RS08065integrase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base . (5/124);  total (5/126)
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.18e-01

Only 100 of 131 total aligned reads displayed.

GAGACCATTTTAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACGCGATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACG  >  NC_007793/1631668‑1631854
                                                                                              |                                                                                            
gagaCCATTTTAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTGCGTTTAGCAACTTGTTCATCTGAAATACTACG‑GAtttt                                                                                        <  1:3284067/100‑1 (MQ=255)
  gaCCATTTTAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTtat                                                                                      <  1:5455136/100‑1 (MQ=255)
  gaCCATTTTAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTtat                                                                                      <  1:5684561/100‑1 (MQ=255)
     cATTTTAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGc                                                                                   <  1:3617214/100‑1 (MQ=255)
     cATTTTAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGc                                                                                   <  1:1742832/100‑1 (MQ=255)
     cATTTTAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGc                                                                                   <  1:3242321/100‑1 (MQ=255)
     cATTTTAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGc                                                                                   <  1:5753590/100‑1 (MQ=255)
     cATTTTAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGc                                                                                   <  1:4323914/100‑1 (MQ=255)
      aTTTTAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCa                                                                                  <  1:1248158/100‑1 (MQ=255)
       ttttAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCag                                                                                 <  1:168071/100‑2 (MQ=255)
       ttttAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCAt                                                                                 <  1:1255337/100‑1 (MQ=255)
       ttttAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCAt                                                                                 <  1:2726938/100‑1 (MQ=255)
        tttAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCaga                                                                                <  1:1454952/100‑3 (MQ=255)
        tttAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCAta                                                                                <  1:4730621/100‑2 (MQ=255)
        tttAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATc                                                                                <  1:5673012/100‑1 (MQ=255)
         ttAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCagag                                                                               <  1:207056/100‑4 (MQ=255)
         ttAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCAtgt                                                                               <  1:1142481/100‑3 (MQ=255)
          tAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGcgcact                                                                              <  1:1960570/100‑6 (MQ=255)
              gAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTc                                                                          <  1:4724573/100‑1 (MQ=255)
              gAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTc                                                                          <  1:229980/100‑1 (MQ=255)
              gAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTc                                                                          <  1:2323652/100‑1 (MQ=255)
              gAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTc                                                                          <  1:580462/100‑1 (MQ=255)
              gAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTc                                                                          <  1:1649135/100‑1 (MQ=255)
              gAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTc                                                                          <  1:4796495/100‑1 (MQ=255)
               aaaGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCa                                                                         <  1:1989470/100‑1 (MQ=255)
               aaaGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCa                                                                         <  1:2458199/100‑1 (MQ=255)
                aaGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAg                                                                        <  1:2902025/100‑1 (MQ=255)
                 aGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCacc                                                                       <  1:764149/100‑3 (MQ=255)
                 aGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAgg                                                                       <  1:2373911/100‑2 (MQ=255)
                 aGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGc                                                                       <  1:1152725/100‑1 (MQ=255)
                  gTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGcg                                                                      <  1:2749710/100‑2 (MQ=255)
                  gTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCt                                                                      <  1:2817535/100‑1 (MQ=255)
                   tttttCTTCAATATATGCTTTTGTTTCATCATTCCAACGTCTACGTTTAGCAACTTGTTCATCTGAAATCCTACG‑GATTTTATAGCATCTTCCTCAGCta                                                                     <  1:2141205/100‑1 (MQ=255)
                   tttttCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTTCGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCta                                                                     <  1:4667443/100‑1 (MQ=255)
                   tttttCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGcct                                                                     <  1:2159769/100‑3 (MQ=255)
                   tttttCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCta                                                                     <  1:917984/100‑1 (MQ=255)
                   tttttCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCta                                                                     <  1:4098326/100‑1 (MQ=255)
                   tttttCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCta                                                                     <  1:501509/100‑1 (MQ=255)
                   tttttCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCta                                                                     <  1:1347552/100‑1 (MQ=255)
                   tttttCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCta                                                                     <  1:4279440/100‑1 (MQ=255)
                   tttttCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCta                                                                     <  1:522558/100‑1 (MQ=255)
                   tttttCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCta                                                                     <  1:1586849/100‑1 (MQ=255)
                   tttttCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCta                                                                     <  1:5297623/100‑1 (MQ=255)
                   tttttCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCta                                                                     <  1:3124980/100‑1 (MQ=255)
                   tttttCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCta                                                                     <  1:5431532/100‑1 (MQ=255)
                    ttttCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCtag                                                                    <  1:1095323/100‑2 (MQ=255)
                    ttttCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCtaa                                                                    <  1:169780/100‑1 (MQ=255)
                    ttttCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCGCAGccta                                                                    <  1:1710351/100‑4 (MQ=255)
                     tttCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGccaga                                                                   <  1:4854687/100‑5 (MQ=255)
                     tttCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCtaag                                                                   <  1:2613655/100‑2 (MQ=255)
                     tttCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCtaag                                                                   <  1:1984691/100‑2 (MQ=255)
                     tttCTTCAATATATGCCTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCtaag                                                                   <  1:1669684/100‑2 (MQ=255)
                      ttcttcAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCtatta                                                                  <  1:4893502/100‑1 (MQ=255)
                      ttcttcAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCtaata                                                                  <  1:1657549/100‑1 (MQ=255)
                      ttcttcAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCtaata                                                                  <  1:4502763/100‑1 (MQ=255)
                      ttcttcAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCtaata                                                                  <  1:4593978/100‑1 (MQ=255)
                      ttcttcAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCtaata                                                                  <  1:4338607/100‑1 (MQ=255)
                      ttcttcAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCtaaag                                                                  <  1:1310858/100‑3 (MQ=255)
                       tcttcAATATATGCTTTTGTTTCATCATTCAACCGTCTCCGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCtaataa                                                                 <  1:416259/100‑1 (MQ=255)
                       tcttcAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCtaatag                                                                 <  1:3273150/100‑2 (MQ=255)
                       tcttcAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCtaataa                                                                 <  1:2115665/100‑1 (MQ=255)
                       tcttcAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCtaataa                                                                 <  1:244736/100‑1 (MQ=255)
                             atatatGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGAtt                                                           <  1:5459094/100‑1 (MQ=255)
                              tatatGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATta                                                          <  1:2577609/100‑1 (MQ=255)
                               atatGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATtat                                                         <  1:240780/100‑1 (MQ=255)
                               atatGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATtat                                                         >  1:3884592/1‑100 (MQ=255)
                               atatGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATtat                                                         >  1:5330900/1‑100 (MQ=255)
                                         gTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCagt                                               <  1:2930781/100‑3 (MQ=255)
                                         gTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCAtt                                               <  1:5774272/100‑2 (MQ=255)
                                         gTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCAtg                                               <  1:1648927/100‑2 (MQ=255)
                                         gTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATc                                               <  1:3525827/100‑1 (MQ=255)
                                          tttCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCa                                              <  1:4798996/100‑1 (MQ=255)
                                            tcatcaTTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAAc                                            <  1:5016752/100‑1 (MQ=255)
                                            tcatcaTTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAAc                                            <  1:901528/100‑1 (MQ=255)
                                            tcatcaTTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAAc                                            <  1:1120068/100‑1 (MQ=255)
                                            tcatcaTTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAAc                                            <  1:4170153/100‑1 (MQ=255)
                                            tcatcaTTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAAc                                            <  1:2217696/100‑1 (MQ=255)
                                            tcatcaTTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAAc                                            <  1:3878458/100‑1 (MQ=255)
                                            tcatcaTTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAAc                                            <  1:3665083/100‑1 (MQ=255)
                                            tcatcaTTCAAACGTCTACGTTTAGCAACTTGTTCATCTCAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACTCCATCAAc                                            <  1:4512263/100‑1 (MQ=255)
                                              atcatTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAACAg                                          <  1:2524801/100‑1 (MQ=255)
                                                catTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAACAGaa                                        <  1:3064388/100‑1 (MQ=255)
                                                       aCGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAACAGAACGCTTag                                 <  1:141356/100‑2 (MQ=255)
                                                            tACGCTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAACAGAACGCTTATTTACt                            <  1:1725317/100‑1 (MQ=255)
                                                             aCGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAACAGAACGCTTATTTACtt                           <  1:3738987/100‑1 (MQ=255)
                                                             aCGTTTAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAACAGAACGCTTATTTACtt                           <  1:1562560/100‑1 (MQ=255)
                                                                 ttAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAACAGAACGCTTATTTACTTctct                       <  1:3540190/100‑1 (MQ=255)
                                                                  tAGCAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAACAGAACGCTTATTTACTTCTCta                      >  1:3232005/1‑100 (MQ=255)
                                                                     cAACTTGTTCATCTGAAATACTACG‑GCTTTTATAGCATCTTCCTCAGCTAATCAAGGATTATAAACACCATCAACAGAACGCTTATTTACTTCTCTATaa                   <  1:2678471/100‑1 (MQ=255)
                                                                     cAACTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAACAGAACGCTTATTTACTTCTCTATaa                   <  1:4044095/100‑1 (MQ=255)
                                                                       aCTTGTTCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAACAGAACGCTTATTTACTTCTCTATAAAg                 <  1:4062827/100‑1 (MQ=255)
                                                                            ttCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAAccg            <  1:5213151/100‑3 (MQ=255)
                                                                             tCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACCGc           <  1:4610856/100‑1 (MQ=255)
                                                                             tCATCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACCGc           <  1:3221178/100‑1 (MQ=255)
                                                                               aTCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACCGcac         <  1:2714126/100‑2 (MQ=255)
                                                                                tCTGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACCGcatc        <  1:5651660/100‑1 (MQ=255)
                                                                                  tGAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACTCCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACCGcatcat      <  1:1611655/100‑1 (MQ=255)
                                                                                   gAAATACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACCGcatcaga     >  1:4077919/1‑98 (MQ=255)
                                                                                      aTACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAATACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACCGCATCATCACg  <  1:2397381/100‑1 (MQ=255)
                                                                                      aTACTACG‑GATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACCGCATCATCACg  <  1:1755993/100‑1 (MQ=255)
                                                                                              |                                                                                            
GAGACCATTTTAATGAAAGTTTTTCTTCAATATATGCTTTTGTTTCATCATTCAAACGTCTACGTTTAGCAACTTGTTCATCTGAAATACTACGCGATTTTATAGCATCTTCCTCAGCTAATAAAGGATTATAAACACCATCAACAGAACGCTTATTTACTTCTCTATAAAGAACTGCATCATCACG  >  NC_007793/1631668‑1631854

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