Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 2,150,261 +C coding (766/960 nt) SAUSA300_RS10960 ← fructokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077932,150,2571.C96.0% 207.3 / 3.3 74coding (770/960 nt)SAUSA300_RS10960fructokinase
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base C (3/68);  total (3/71)
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.82e-01

GCTAAAATTTCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑CCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGGCAGCACCATCTGCTCCCTTCGTAT  >  NC_007793/2150158‑2150353
                                                                                                    |                                                                                                 
gCTAAAATTTCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTg                                                                                                    <  1:3018143/100‑1 (MQ=255)
 ctcttccgatcTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑c                                                                                                  <  1:1639192/90‑1 (MQ=255)
 cTAAAATTTCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTgg                                                                                                   <  1:5434511/100‑2 (MQ=255)
  tAAAATTTCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑cg                                                                                                 <  1:2203077/100‑2 (MQ=255)
  tAAAATTTCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑cc                                                                                                 <  1:140026/100‑1 (MQ=255)
  tAAAATTTCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑ca                                                                                                 <  1:3230564/100‑2 (MQ=255)
  tAAAATTTCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑ca                                                                                                 <  1:2186649/100‑2 (MQ=255)
       tttCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑ccc‑ccag                                                                                           <  1:2392129/100‑4 (MQ=255)
       tttCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑ccc‑ccag                                                                                           <  1:2166916/100‑4 (MQ=255)
        ttCCTCGCCTTCATTTTCAAATAACTCAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCATTAAATGCATCCCCTG‑ccc‑ccagt                                                                                          <  1:4350985/100‑5 (MQ=255)
        ttCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑ccc‑ccagt                                                                                          <  1:4572004/100‑5 (MQ=255)
        ttCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑ccc‑ccagt                                                                                          <  1:3034826/100‑5 (MQ=255)
        ttCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑ccc‑ccagt                                                                                          <  1:2091097/100‑5 (MQ=255)
        ttCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑ccc‑ccagg                                                                                          <  1:3783179/100‑5 (MQ=255)
        tctcTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑ccc‑ccagt                                                                                          <  1:5506521/97‑5 (MQ=255)
         tCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCTGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑ccc‑ccagtt                                                                                         <  1:1219221/100‑6 (MQ=255)
         gatctaCCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑ccc‑ccagtg                                                                                         <  1:3985181/95‑6 (MQ=255)
         ctcTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑ccc‑ccagta                                                                                         <  1:122336/98‑6 (MQ=255)
          ccTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTg                                                                                         <  1:545796/100‑1 (MQ=255)
          ccTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTg                                                                                         <  1:2196192/100‑1 (MQ=255)
           ttCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGt                                                                                        <  1:3284022/99‑1 (MQ=255)
           ttCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGt                                                                                        <  1:3213581/99‑1 (MQ=255)
           tctaCCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGt                                                                                        <  1:5406234/97‑1 (MQ=255)
           gatctcctCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGt                                                                                        <  1:1650537/93‑1 (MQ=255)
           cgatctttCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGt                                                                                        <  1:3897448/94‑1 (MQ=255)
           ccgatcTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATTACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGt                                                                                        <  1:922282/95‑1 (MQ=255)
           cTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTga                                                                                        <  1:1745189/100‑2 (MQ=255)
           cTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGt                                                                                        <  1:3561677/100‑1 (MQ=255)
           cTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGt                                                                                        <  1:807526/100‑1 (MQ=255)
           cTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGt                                                                                        <  1:846410/100‑1 (MQ=255)
           cTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGt                                                                                        <  1:973296/100‑1 (MQ=255)
           atctccTTCGTTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGt                                                                                        <  1:1006212/96‑1 (MQ=255)
           atctccTTCATTTTCATATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGt                                                                                        <  1:1699936/96‑1 (MQ=255)
           atctccTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGt                                                                                        <  1:1931783/96‑1 (MQ=255)
             tctcTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAGTCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTAt                                                                                      <  1:2458564/97‑1 (MQ=255)
             cACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTag                                                                                      <  1:3826372/100‑2 (MQ=255)
             cACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTaa                                                                                      <  1:1995210/100‑2 (MQ=255)
              tccTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATc                                                                                     <  1:1193784/99‑1 (MQ=255)
              ctcTTCATTATCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATc                                                                                     <  1:3353042/98‑1 (MQ=255)
              aCCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTAta                                                                                     <  1:4712415/100‑2 (MQ=255)
              aCCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATc                                                                                     <  1:1072875/100‑1 (MQ=255)
              aCCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATc                                                                                     <  1:3526047/100‑1 (MQ=255)
              aCCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATc                                                                                     <  1:746615/100‑1 (MQ=255)
              aCCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATc                                                                                     <  1:3254800/100‑1 (MQ=255)
              aCCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATc                                                                                     <  1:4941702/100‑1 (MQ=255)
               ccTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAGTAAATGCATCCCCTGCCCC‑CAGTTGTATCa                                                                                    <  1:250641/100‑1 (MQ=255)
               ccTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCa                                                                                    <  1:2081437/100‑1 (MQ=255)
               ccTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCa                                                                                    <  1:2529013/100‑1 (MQ=255)
               ccTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCa                                                                                    <  1:2987927/100‑1 (MQ=255)
                cTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCac                                                                                   <  1:1181944/100‑2 (MQ=255)
                cTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCaa                                                                                   <  1:1944197/100‑1 (MQ=255)
                     ttttCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAgga                                                                              <  1:1646395/100‑2 (MQ=255)
                               aCTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCGTCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTtat                                                                    <  1:1414882/100‑1 (MQ=255)
                               aCTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTtat                                                                    <  1:5482870/100‑1 (MQ=255)
                                cTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTCACCTtatg                                                                   <  1:455165/100‑2 (MQ=255)
                                   aGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGcc                                                                <  1:3594105/100‑1 (MQ=255)
                                   aGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGcc                                                                <  1:1750164/100‑1 (MQ=255)
                                   aGTTAAATTTAACACATCGGTCGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTGCCTTATAGcc                                                                <  1:2847468/100‑1 (MQ=255)
                                     ttAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCAc                                                              <  1:4757332/100‑1 (MQ=255)
                                      tAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACt                                                             <  1:2981665/100‑1 (MQ=255)
                                      tAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACt                                                             <  1:3286619/100‑1 (MQ=255)
                                      tAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACt                                                             <  1:4248510/100‑1 (MQ=255)
                                      tAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTCTTACCTTATAGccgca                                                             <  1:2004135/100‑4 (MQ=255)
                                       aaaTTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCagat                                                            <  1:3437886/100‑4 (MQ=255)
                                       aaaTTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCAcga                                                            <  1:3962756/100‑3 (MQ=255)
                                       aaaTTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACtg                                                            <  1:3832339/100‑2 (MQ=255)
                                       aaaTTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTa                                                            <  1:3141153/100‑1 (MQ=255)
                                       aaaTTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTa                                                            <  1:5216329/100‑1 (MQ=255)
                                       aaaTTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑CCCTCAGTTGTATCAACAGGTTTTACCTTATAGCCACTa                                                            <  1:2828179/100‑1 (MQ=255)
                                        aaTTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTAt                                                           <  1:4964075/100‑1 (MQ=255)
                                                  aTCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGAtt                                                 <  1:2147141/100‑1 (MQ=255)
                                                      gTTGCTAAAATCCTGCTAATAACTGCACCAATAGATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATa                                             >  1:2933826/1‑100 (MQ=255)
                                                        tGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATAcc                                           <  1:736631/100‑1 (MQ=255)
                                                        tGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATAcc                                           <  1:334426/100‑1 (MQ=255)
                                                                 ccTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCaa                                  >  1:2037950/1‑100 (MQ=255)
                                                                          atcTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgc                         <  1:2677764/98‑1 (MQ=255)
                                                                          aaCTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgc                         <  1:5040634/100‑1 (MQ=255)
                                                                          aaCTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgc                         <  1:169308/100‑1 (MQ=255)
                                                                          aaCTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgc                         <  1:5318183/100‑1 (MQ=255)
                                                                          aaCTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgc                         <  1:87780/100‑1 (MQ=255)
                                                                          aaCTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgc                         <  1:3200168/100‑1 (MQ=255)
                                                                          aaCTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgc                         <  1:1396661/100‑1 (MQ=255)
                                                                           aCTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACggag                        <  1:3341603/100‑3 (MQ=255)
                                                                            cTGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAAAGATTAATACCATTTTTCAAATATACGgcaa                       <  1:261284/100‑2 (MQ=255)
                                                                             tGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgcagc                      <  1:5381832/100‑1 (MQ=255)
                                                                             tGCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgcagc                      <  1:3468489/100‑1 (MQ=255)
                                                                              gCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTTCCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgcagat                     <  1:647004/100‑3 (MQ=255)
                                                                              gCACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgcagca                     >  1:3655011/1‑100 (MQ=255)
                                                                               cACCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgcagcaa                    <  1:1766250/100‑2 (MQ=255)
                                                                                aCCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGGCAGCAcc                   <  1:4695368/100‑1 (MQ=255)
                                                                                aCCAATAAATGCATCCCCTGCCCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGGCAGCAcc                   <  1:3500209/100‑1 (MQ=255)
                                                                                           tccgatcTG‑CCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGGCAGCACCATCTGCTCCCtt       <  1:1431623/94‑1 (MQ=255)
                                                                                                tcTG‑CCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGGCAGCACCATCTGCTCCCTTCGtag  <  1:3132012/99‑2 (MQ=255)
                                                                                                tcTG‑CCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGGCAGCACCATCTGCTCCCTTCGtag  <  1:2677955/99‑2 (MQ=255)
                                                                                                    |                                                                                                 
GCTAAAATTTCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑CCC‑CAGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGGCAGCACCATCTGCTCCCTTCGTAT  >  NC_007793/2150158‑2150353

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑