Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 1,846,346 A→T R186R (CGT→CGA)  SAUSA300_RS09160 ← iron‑regulated surface determinant protein H

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077931,846,3460AT98.9% 247.7 / NA 92R186R (CGT→CGA) SAUSA300_RS09160iron‑regulated surface determinant protein H
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (6/85);  total (6/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTAATTTAGTATAATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAAACGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCATAGACTTCAAACTTCTTCCATGT  >  NC_007793/1846250‑1846441
                                                                                                |                                                                                               
tACTAATTCAGTATAATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATGCTTGTGTGCCATCTGAAACTGGGAATCGa                                                                                              <  1:4624465/100‑1 (MQ=255)
    aaTTTAGTATAATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAatat                                                                                          <  1:4825952/100‑1 (MQ=255)
     aTTTAGTATAATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTGTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAatatc                                                                                         <  1:454508/100‑2 (MQ=255)
      tttAGTATAATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAataggt                                                                                        <  1:655100/100‑4 (MQ=255)
        tAGTATAATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGc                                                                                      <  1:971449/100‑1 (MQ=255)
        tAGTATAATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGc                                                                                      <  1:2397986/100‑1 (MQ=255)
        tAGTATAATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGc                                                                                      <  1:4252363/100‑1 (MQ=255)
        tAGTATAATCATAATGTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGc                                                                                      <  1:2569648/100‑1 (MQ=255)
         aGTATAATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGcc                                                                                     <  1:653643/100‑2 (MQ=255)
          gTATAATCATAATTTGTTTCTTCTCCATCATCAATTTGGGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGctc                                                                                    <  1:4481365/100‑3 (MQ=255)
          gTATAATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCAt                                                                                    <  1:6018759/100‑1 (MQ=255)
          gTATAATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGAGAATCGAATATAGGCAt                                                                                    <  1:1740141/100‑1 (MQ=255)
          gTATAATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAAGTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCAt                                                                                    <  1:5827380/100‑1 (MQ=255)
           tataATCATAATTTGTTTTTTCTCCATCATCAATTTGAGTCGACGAAACAATTTTCAATTCTAGTGGGCGATCTGAAACTGTGAATCGAATATAGGCAta                                                                                   <  1:2382689/100‑2 (MQ=255)
           tataATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCAtt                                                                                   <  1:6547680/100‑2 (MQ=255)
           tataATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATg                                                                                   <  1:6351094/100‑1 (MQ=255)
           tataATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATg                                                                                   <  1:5484774/100‑1 (MQ=255)
           tataATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATg                                                                                   <  1:4793642/100‑1 (MQ=255)
           tataATCATAATTTGTTTCTCCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATg                                                                                   <  1:2072522/100‑1 (MQ=255)
            ataATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATgg                                                                                  <  1:1586061/100‑1 (MQ=255)
            ataATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATgg                                                                                  <  1:3377008/100‑1 (MQ=255)
            ataATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATgg                                                                                  <  1:5775132/100‑1 (MQ=255)
            ataATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATgg                                                                                  <  1:3088089/100‑1 (MQ=255)
            ataATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAAACGAATATAAGCATgc                                                                                  <  1:1412335/100‑2 (MQ=255)
            ataATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTCGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATgg                                                                                  <  1:6335154/100‑1 (MQ=255)
            ataATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCCTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATgg                                                                                  <  1:1082087/100‑1 (MQ=255)
            ataATCATAATTTGTTTCTGCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATgg                                                                                  <  1:1620008/100‑1 (MQ=255)
             taATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATggc                                                                                 <  1:2649121/100‑2 (MQ=255)
              aaTCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGtt                                                                                <  1:2383456/100‑2 (MQ=255)
               aTCATAATTTGTTTCTTCTCCATCATCGATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGtct                                                                               <  1:936280/100‑1 (MQ=255)
               aTCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGtct                                                                               <  1:5831562/100‑1 (MQ=255)
               aTCATAATTGGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGtct                                                                               <  1:3545896/100‑1 (MQ=255)
                tCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGtctg                                                                              <  1:2552495/100‑2 (MQ=255)
                tCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGtcta                                                                              <  1:1762384/100‑2 (MQ=255)
                  aTAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATggggcgta                                                                            <  1:4309186/100‑7 (MQ=255)
                    aaTTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTgg                                                                          <  1:6213963/100‑1 (MQ=255)
                     gctcttccgatctcCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGt                                                                         <  1:3951556/90‑1 (MQ=255)
                     aTTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGt                                                                         <  1:3199467/100‑1 (MQ=255)
                     aTTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGt                                                                         <  1:419743/100‑1 (MQ=255)
                     aTTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGt                                                                         <  1:3927779/100‑1 (MQ=255)
                     aTTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGt                                                                         <  1:5992350/100‑1 (MQ=255)
                     aTTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGt                                                                         <  1:3104425/100‑1 (MQ=255)
                      tttGTTTCTTCTCCATCATCAATTTTAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTa                                                                        <  1:6172715/100‑1 (MQ=255)
                       ttGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTAc                                                                       <  1:3591750/100‑1 (MQ=255)
                       ttGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTAc                                                                       <  1:509624/100‑1 (MQ=255)
                       ttGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTAc                                                                       <  1:3065746/100‑1 (MQ=255)
                       ttGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTAc                                                                       <  1:4031690/100‑1 (MQ=255)
                       ttGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTAc                                                                       <  1:3205098/100‑1 (MQ=255)
                       ttGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTAc                                                                       <  1:1443983/100‑1 (MQ=255)
                       ttGTTTCTTCTCCACCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTAc                                                                       <  1:4786249/100‑1 (MQ=255)
                       tcttccgatctcCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTAc                                                                       <  1:6380923/92‑1 (MQ=255)
                        tGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTAcc                                                                      <  1:2781377/100‑2 (MQ=255)
                         gTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACaa                                                                     <  1:3521214/100‑2 (MQ=255)
                         gTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAg                                                                     <  1:1212565/100‑1 (MQ=255)
                         gTTTCTTCTACATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAg                                                                     <  1:1994167/100‑1 (MQ=255)
                          tttCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAgg                                                                    <  1:1220958/100‑1 (MQ=255)
                          tttCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAgg                                                                    <  1:6415256/100‑1 (MQ=255)
                               tctcCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACtat                                                               <  1:1337812/100‑1 (MQ=255)
                                 tcCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACtaggt                                                             <  1:4762398/100‑4 (MQ=255)
                                             ttGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATGTGATACAAGtctc                                                 <  1:5796418/100‑1 (MQ=255)
                                              tGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCa                                                <  1:3299382/100‑1 (MQ=255)
                                                  tCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTgg                                            <  1:3981976/100‑1 (MQ=255)
                                                   cGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTggg                                           <  1:2520791/100‑2 (MQ=255)
                                                   cGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTgga                                           <  1:2123490/100‑2 (MQ=255)
                                                      agaaACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCagg                                        <  1:2967455/100‑3 (MQ=255)
                                                           cAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTg                                   <  1:4055106/100‑1 (MQ=255)
                                                           cAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTg                                   <  1:3139453/100‑1 (MQ=255)
                                                                ttCAATTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGATACAGGACTATATGATACAAGTCTCACTGG‑AATTTTTGATTGtt                             <  1:5805275/100‑1 (MQ=255)
                                                                    aTTCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTc                          <  1:6002268/100‑1 (MQ=255)
                                                                     ttCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGCTTGTTTTCa                         >  1:3427910/1‑100 (MQ=255)
                                                                      tCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCAt                        <  1:4124533/100‑1 (MQ=255)
                                                                      tCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCAt                        <  1:5522542/100‑1 (MQ=255)
                                                                      tCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCAt                        <  1:6473289/100‑1 (MQ=255)
                                                                      tCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCAt                        <  1:6364950/100‑1 (MQ=255)
                                                                      tCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCAt                        <  1:2784413/100‑1 (MQ=255)
                                                                      tCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCAt                        <  1:6049127/100‑1 (MQ=255)
                                                                      tCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCAt                        <  1:6110988/100‑1 (MQ=255)
                                                                      tCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCAt                        <  1:4257797/100‑1 (MQ=255)
                                                                      tCTTGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCAt                        <  1:3073726/100‑1 (MQ=255)
                                                                        ttGTGTGCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCATAg                      <  1:21753/100‑1 (MQ=255)
                                                                              gCCATCTGAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCATAGACTTCa                <  1:3641513/100‑1 (MQ=255)
                                                                                     gAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATCTGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCATAGACTTCAAActtct         >  1:6371815/1‑100 (MQ=255)
                                                                                     gAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCATAGACTTCAAActtct         >  1:3519680/1‑100 (MQ=255)
                                                                                     gAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCATAGACTTCAAActtct         >  1:6585256/1‑100 (MQ=255)
                                                                                     gAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCATAGACTTCAAActtct         <  1:1579368/100‑1 (MQ=255)
                                                                                     gAAACTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCATAGACTTCAAActtct         >  1:3905964/1‑100 (MQ=255)
                                                                                        aCTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCATAGACTTCAAACTTCTTcc      <  1:1710007/100‑1 (MQ=255)
                                                                                        aCTGGGAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCATAGACTTCAAACTTCTTcc      <  1:1695269/100‑1 (MQ=255)
                                                                                           gggAATCGAATATGGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCATAGACTTCAAACTTCTTCCATg   <  1:3972409/100‑1 (MQ=255)
                                                                                           gggAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCATAGACTTCAAACTTCTTCCATg   >  1:6652516/1‑100 (MQ=255)
                                                                                           gggAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCATAGACTTCAAACTTCTTCCATg   <  1:4863769/100‑1 (MQ=255)
                                                                                            ggAATCGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCATAGACTTCAAACTTCTTCCATGt  <  1:4763403/100‑1 (MQ=255)
                                                                                                |                                                                                               
TACTAATTTAGTATAATCATAATTTGTTTCTTCTCCATCATCAATTTGAGTCGAAGAAACAATTTTCAATTCTTGTGTGCCATCTGAAACTGGGAAACGAATATAGGCATGGTCTTCTGGTACAGGACTATATGATACAAGTCTCACTGGCAATTTTTGATTGTTTTCATAGACTTCAAACTTCTTCCATGT  >  NC_007793/1846250‑1846441

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