Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,580,153 A→G V213A (GTG→GCG)  ydeM ← putative anaerobic sulfatase maturation enzyme YdeM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,580,1530AG100.0% 1509.9 / NA 359V213A (GTG→GCG) ydeMputative anaerobic sulfatase maturation enzyme YdeM
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (176/183);  total (176/183)

Only 100 of 364 total aligned reads displayed.

CAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTT  >  NC_000913/1580055‑1580250
                                                                                                  |                                                                                                 
cAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGcgca                                                                                                 >  1:2575243/1‑101 (MQ=255)
 aaaGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAg                                                                                                >  2:2575244/1‑101 (MQ=255)
 aaaGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAg                                                                                                <  1:2574623/101‑1 (MQ=255)
 aaaGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAg                                                                                                <  2:2574716/101‑1 (MQ=255)
     cTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGaaaa                                                                                            <  2:2574382/101‑1 (MQ=255)
     cTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTCCCATAAGCCGTGGGAGGCGCAGaaaa                                                                                            <  2:2574607/101‑1 (MQ=255)
      tttCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGaaaaa                                                                                           >  2:2575249/1‑101 (MQ=255)
         cAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCa                                                                                        >  2:2575253/1‑101 (MQ=255)
          aaaCTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCaa                                                                                       >  1:2575255/1‑101 (MQ=255)
           aaCTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAAt                                                                                      <  2:2574825/101‑1 (MQ=255)
            aCTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATg                                                                                     >  2:2575258/1‑101 (MQ=255)
             cTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGa                                                                                    >  2:2575262/1‑101 (MQ=255)
             cTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGa                                                                                    <  2:2574699/101‑1 (MQ=255)
              tGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGAt                                                                                   >  2:2575266/1‑101 (MQ=255)
                aCGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATAc                                                                                 <  2:2574690/101‑1 (MQ=255)
                aCGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATAc                                                                                 <  1:2574629/101‑1 (MQ=255)
                    aTGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGaa                                                                             <  2:2574585/101‑1 (MQ=255)
                      gAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAAtg                                                                           >  2:2575272/1‑101 (MQ=255)
                       aaaaTTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAAtgt                                                                          >  1:2575274/1‑101 (MQ=255)
                         aaTTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATgtgt                                                                        <  2:2575028/101‑1 (MQ=255)
                         aaTTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATgtgt                                                                        >  1:2575279/1‑101 (MQ=255)
                         aaTTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATgtgt                                                                        <  1:2574768/101‑1 (MQ=255)
                          aTTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAAtg                                                                           <  1:2574545/97‑1 (MQ=255)
                          aTTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGtt                                                                       <  2:2574800/101‑1 (MQ=255)
                            ttCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTtct                                                                     >  2:2575283/1‑101 (MQ=255)
                             tCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTtctc                                                                    >  1:2575284/1‑101 (MQ=255)
                             tCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTtctc                                                                    <  2:2574734/101‑1 (MQ=255)
                              cACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCa                                                                   <  2:2574465/101‑1 (MQ=255)
                              cACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCa                                                                   >  2:2575287/1‑101 (MQ=255)
                                 ccACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTa                                                                >  2:2575292/1‑101 (MQ=255)
                                 ccACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTa                                                                <  1:2574755/101‑1 (MQ=255)
                                 ccACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTa                                                                >  1:2575294/1‑101 (MQ=255)
                                   acaTCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATg                                                              <  1:2574852/101‑1 (MQ=255)
                                     aTCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGAc                                                            <  2:2574867/101‑1 (MQ=255)
                                     aTCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGAc                                                            <  1:2574816/101‑1 (MQ=255)
                                      tCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGAcc                                                           >  1:2575298/1‑101 (MQ=255)
                                      tCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGAcc                                                           <  2:2574625/101‑1 (MQ=255)
                                      tCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGAcc                                                           <  2:2574894/101‑1 (MQ=255)
                                       cGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCa                                                          >  2:2575301/1‑101 (MQ=255)
                                        gTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCAc                                                         <  1:2574731/101‑1 (MQ=255)
                                        gTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCAc                                                         <  1:2574829/101‑1 (MQ=255)
                                        gTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCAc                                                         >  2:2575305/1‑101 (MQ=255)
                                        gTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCAc                                                         >  1:2575306/1‑101 (MQ=255)
                                        gTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCAc                                                         <  1:2574964/101‑1 (MQ=255)
                                         tttttAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACt                                                        <  1:2574744/101‑1 (MQ=255)
                                            ttAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGaa                                                     >  1:2575313/1‑101 (MQ=255)
                                            ttAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGaa                                                     >  2:2575312/1‑101 (MQ=255)
                                               aCCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATc                                                  >  2:2575319/1‑101 (MQ=255)
                                                  cATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAat                                               >  1:2575322/1‑101 (MQ=255)
                                                    ttGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAatat                                             >  1:2575327/1‑101 (MQ=255)
                                                     tGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATAtt                                            <  1:2575107/101‑1 (MQ=255)
                                                       cATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCg                                          >  1:2575331/1‑101 (MQ=255)
                                                       cATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCg                                          >  2:2575330/1‑101 (MQ=255)
                                                         tAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGc                                        <  2:2574673/101‑1 (MQ=255)
                                                          aaaaaTGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCg                                       >  1:2575339/1‑101 (MQ=255)
                                                           aaaaTGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGt                                      <  2:2574900/101‑1 (MQ=255)
                                                            aaaTGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTc                                     >  2:2575340/1‑101 (MQ=255)
                                                             aaTGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTcc                                    >  2:2575344/1‑101 (MQ=255)
                                                                ggTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCg                                  >  1:2575348/1‑100 (MQ=255)
                                                                 gTTGACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGt                                <  2:2574335/101‑1 (MQ=255)
                                                                    gACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTc                             >  1:2575353/1‑101 (MQ=255)
                                                                    gACATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTc                             >  2:2575354/1‑101 (MQ=255)
                                                                      cATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTa                           <  1:2574899/101‑1 (MQ=255)
                                                                      cATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTa                           <  2:2574775/101‑1 (MQ=255)
                                                                      cATAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTa                           <  1:2574610/101‑1 (MQ=255)
                                                                       aTAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAg                          <  1:2574664/101‑1 (MQ=255)
                                                                       aTAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAg                          <  2:2574843/101‑1 (MQ=255)
                                                                       aTAAACTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAg                          >  2:2575361/1‑101 (MQ=255)
                                                                         aaaCTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCa                        <  1:2574443/101‑1 (MQ=255)
                                                                         aaaCTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCa                        >  1:2575364/1‑101 (MQ=255)
                                                                          aaCTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCaa                       >  1:2575367/1‑101 (MQ=255)
                                                                           aCTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAAt                      >  2:2575369/1‑101 (MQ=255)
                                                                           aCTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAAt                      >  1:2575372/1‑101 (MQ=255)
                                                                            cTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAAtt                     <  1:2575010/101‑1 (MQ=255)
                                                                            cTTGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAAtt                     <  1:2575180/101‑1 (MQ=255)
                                                                             ttGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTc                    <  2:2574519/101‑1 (MQ=255)
                                                                             ttGCCATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATCCGGCGTCCCGGTTTCTAGCAATTc                    >  2:2575377/1‑101 (MQ=255)
                                                                                ccATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGAt                 <  1:2574873/101‑1 (MQ=255)
                                                                                ccATAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGAt                 >  2:2575383/1‑101 (MQ=255)
                                                                                  aTAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATaa               <  1:2574915/101‑1 (MQ=255)
                                                                                  aTAAGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATaa               >  1:2575385/1‑101 (MQ=255)
                                                                                    aaGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAAt             >  2:2575390/1‑101 (MQ=255)
                                                                                    aaGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAAt             >  1:2575388/1‑101 (MQ=255)
                                                                                     aGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAAtt            <  2:2574789/101‑1 (MQ=255)
                                                                                     aGCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAAtt            >  1:2575394/1‑101 (MQ=255)
                                                                                      gCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAATTTg           <  1:2574491/101‑1 (MQ=255)
                                                                                      gCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTg           <  2:2574636/101‑1 (MQ=255)
                                                                                      gCCGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTg           <  2:2574880/101‑1 (MQ=255)
                                                                                       ccGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGc          <  1:2573728/101‑1 (MQ=255)
                                                                                        cGTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCa         >  1:2575398/1‑101 (MQ=255)
                                                                                         gTGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCat        >  1:2575403/1‑101 (MQ=255)
                                                                                          tGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCata       >  2:2575406/1‑101 (MQ=255)
                                                                                          tGGGAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCata       >  1:2575409/1‑101 (MQ=255)
                                                                                           gggAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCatat      <  2:2574655/101‑1 (MQ=255)
                                                                                            ggAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATg     <  1:2574797/101‑1 (MQ=255)
                                                                                            ggAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATg     <  2:2574750/101‑1 (MQ=255)
                                                                                             gAGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGt    <  2:2575124/101‑1 (MQ=255)
                                                                                              aGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGtt   >  1:2575415/1‑101 (MQ=255)
                                                                                              aGGCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGtt   <  2:2574559/101‑1 (MQ=255)
                                                                                               ggCGCAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGttt  >  1:2575419/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
CAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTATGACCACTGAAATCAATATTCGGCGTCCCGGTTTCTAGCAATTCGATAAATTGCATATGTTT  >  NC_000913/1580055‑1580250

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