Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,102,221 T→C E107E (GAA→GAG)  csgE ← curlin secretion specificity factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,102,2210TC100.0% 255.9 / NA 72E107E (GAA→GAG) csgEcurlin secretion specificity factor
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (37/35);  total (37/35)

CGAAATAAGCATGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCC  >  NC_000913/1102086‑1102351
                                                                                                                                       |                                                                                                                                  
cGAAATAAGCATGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGcctct                                                                                                                                 >  1:1639326/1‑139 (MQ=255)
  aaaTAAGCATGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCctcttc                                                                                                                               >  2:525529/1‑139 (MQ=255)
  aaaTAAGCATGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTAATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCctcttc                                                                                                                               >  1:522895/1‑139 (MQ=255)
           tGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTGAGTTTGAAt                                                                                                                      <  2:1262417/139‑1 (MQ=255)
           tGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                                      <  2:1639326/139‑1 (MQ=255)
           tGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                                      <  2:1004748/139‑1 (MQ=255)
           tGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                                      <  2:1473464/139‑1 (MQ=255)
           tGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                                      <  1:1307695/139‑1 (MQ=255)
           tGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAAt                                                                                                                      <  2:522895/139‑1 (MQ=255)
                   actactGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCaa                                                                                                              >  2:649833/1‑139 (MQ=255)
                       ctGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGcctctt                                                                                                                                >  1:297941/1‑117 (MQ=255)
                       ctGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGcctctt                                                                                                                                <  2:297941/117‑1 (MQ=255)
                       ctGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAgac                                                                                                          >  2:86290/1‑139 (MQ=255)
                         gCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAgacga                                                                                                        >  2:61561/1‑139 (MQ=255)
                                 gACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTtc                                                                                                >  1:1048810/1‑139 (MQ=255)
                                        aTTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGt                                                                                         >  1:558418/1‑139 (MQ=255)
                                                  ccGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCaa                                                                               >  1:918599/1‑139 (MQ=255)
                                                   cGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGa                                                                                            >  2:286454/1‑125 (MQ=255)
                                                   cGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGa                                                                                            <  1:286454/125‑1 (MQ=255)
                                                          ttttttATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATa                                                                       >  2:1410464/1‑139 (MQ=255)
                                                                   taGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCt                                                              <  1:1691750/139‑1 (MQ=255)
                                                                          catcatGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGcctctt                                                                                                                                <  1:842749/66‑1 (MQ=255)
                                                                          catcatGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGcctctt                                                                                                                                >  2:842749/1‑66 (MQ=255)
                                                                           atcatGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCGAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATaaaaaa                                                                  >  2:1248079/1‑127 (MQ=255)
                                                                           atcatGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCGAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATaaaaaa                                                                  <  1:1248079/127‑1 (MQ=255)
                                                                           atcatGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCg                                                                                             >  2:639563/1‑100 (MQ=255)
                                                                           atcatGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCg                                                                                             <  1:639563/100‑1 (MQ=255)
                                                                           atcatGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATa                                                      >  1:1764205/1‑139 (MQ=255)
                                                                              atGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACg                                                   <  1:940082/139‑1 (MQ=255)
                                                                                gcgcCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTc                                                               >  2:1246320/1‑125 (MQ=255)
                                                                                gcgcCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTc                                                               <  1:1246320/125‑1 (MQ=255)
                                                                                gcgcCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTc                                                 >  2:369072/1‑139 (MQ=255)
                                                                                 cgcCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCaaa                                                                                                             <  2:559200/78‑1 (MQ=255)
                                                                                 cgcCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCaaa                                                                                                             >  1:559200/1‑78 (MQ=255)
                                                                                   ccAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAAc                                                                              <  1:9384/107‑1 (MQ=255)
                                                                                   ccAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAAc                                                                              >  2:9384/1‑107 (MQ=255)
                                                                                   ccAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTg                                              >  1:1685880/1‑139 (MQ=255)
                                                                                     aaaTCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGcctctt                                                                                                                                <  2:384718/55‑1 (MQ=255)
                                                                                     aaaTCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGcctctt                                                                                                                                >  1:384718/1‑55 (MQ=255)
                                                                                      aaTCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                                           <  1:86290/139‑1 (MQ=255)
                                                                                      aaTCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                                           <  2:1685880/139‑1 (MQ=255)
                                                                                      aaTCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                                           >  1:708836/1‑139 (MQ=255)
                                                                                       aTCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTg                                          >  1:944094/1‑139 (MQ=255)
                                                                                       aTCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTg                                          <  2:918599/139‑1 (MQ=255)
                                                                                         cGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGAc                                        <  1:61561/139‑1 (MQ=255)
                                                                                           cccGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGa                                                                                            >  2:210826/1‑85 (MQ=255)
                                                                                           cccGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGa                                                                                            <  1:210826/85‑1 (MQ=255)
                                                                                           cccGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCg                                      <  1:1410464/139‑1 (MQ=255)
                                                                                              gTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAGAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTta                                   >  2:592273/1‑139 (MQ=255)
                                                                                                aCTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAg                                                                                                     >  1:1650281/1‑71 (MQ=255)
                                                                                                aCTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAg                                                                                                     <  2:1650281/71‑1 (MQ=255)
                                                                                                aCTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAAt                                                       >  1:1691557/1‑117 (MQ=255)
                                                                                                aCTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAAt                                                       <  2:1691557/117‑1 (MQ=255)
                                                                                                     ataaCGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGAt                            >  2:1577795/1‑139 (MQ=255)
                                                                                                         cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAg                                                                                                     >  2:1262200/1‑62 (MQ=255)
                                                                                                         cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAg                                                                                                     <  1:1262200/62‑1 (MQ=255)
                                                                                                         cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACgg                                                                            >  1:1541920/1‑87 (MQ=255)
                                                                                                         cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACgg                                                                            <  2:1541920/87‑1 (MQ=255)
                                                                                                         cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                                           >  1:1645370/1‑120 (MQ=255)
                                                                                                         cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGAtt                                           <  2:1645370/120‑1 (MQ=255)
                                                                                                         cGCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAg                        <  2:708836/139‑1 (MQ=255)
                                                                                                          gCCTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGc                       >  2:565424/1‑139 (MQ=255)
                                                                                                           ccTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTg                                          >  2:1828043/1‑119 (MQ=255)
                                                                                                           ccTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTg                                          <  1:1828043/119‑1 (MQ=255)
                                                                                                           ccTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCt                      >  2:1224581/1‑139 (MQ=255)
                                                                                                            cTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGt                                                                                         >  1:738378/1‑71 (MQ=255)
                                                                                                            cTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGt                                                                                         <  2:738378/71‑1 (MQ=255)
                                                                                                            cTGATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCtt                     >  1:494890/1‑139 (MQ=255)
                                                                                                              gATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAGAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTcc                   <  1:592273/139‑1 (MQ=255)
                                                                                                              gATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGa                                                           <  2:1698592/99‑1 (MQ=255)
                                                                                                              gATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGa                                                           >  1:1698592/1‑99 (MQ=255)
                                                                                                              gATTTATCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTcc                   <  1:1577795/139‑1 (MQ=255)
                                                                                                                   aTCTGGCGACGATTTAGTGCCTCTCCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATc              <  2:1048810/139‑1 (MQ=255)
                                                                                                                    tCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAgacga                                                                                                        >  2:1441532/1‑48 (MQ=255)
                                                                                                                    tCTGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAgacga                                                                                                        <  1:1441532/48‑1 (MQ=255)
                                                                                                                      tGGCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATCGTg           >  1:734241/1‑139 (MQ=255)
                                                                                                                        gCGACGATTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCa         <  1:649833/139‑1 (MQ=255)
                                                                                                                        gCGACGATTTAGTGCCTCTTCAGTTTGAATCAATGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGTAAAATAACGACCTGATTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCa         <  1:369072/139‑1 (MQ=255)
                                                                                                                         cgacgaTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAAc                                                                              >  1:125357/1‑69 (MQ=255)
                                                                                                                         cgacgaTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAAc                                                                              <  2:125357/69‑1 (MQ=255)
                                                                                                                         cgacgaTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCAc        >  1:656042/1‑139 (MQ=255)
                                                                                                                              aTTTAGTGCCTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGG‑AAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGcc  >  1:972888/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                  
CGAAATAAGCATGAGTAGAACTACTGCATGTTTGACACGCATTGCAGCCTCCGAACAATTTTTTATTTAGAATTCATCATGCGCCAAATCGCCCGTACTTAATAACGCCTGATTTATCTGGCGACGATTTAGTGCTTCTTCAGTTTGAATCAGTGCAAAGACGACAGTTTTCTCGAAGTCTCTTTTCAACGGAAATAAAAAAGTCTGGAAAATAACGTCCTGATTGACCGTTATAGTGATCCAGCTTCCCCATCGTGCACTGGGCC  >  NC_000913/1102086‑1102351

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