Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,102,451 A→G intergenic (‑32/+435) csgD ← / ← b4518 CsgD/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,102,4510AG100.0% 1323.9 / NA 403intergenic (‑32/+435)csgD/b4518CsgD/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (177/226);  total (177/226)

Only 100 of 414 total aligned reads displayed.

GAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATGGGTA  >  NC_000913/1102361‑1102540
                                                                                          |                                                                                         
gaagaTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGcc                                                                                         >  1:5539647/1‑93 (MQ=255)
    aTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGAc                                                                                     >  1:1664496/1‑93 (MQ=255)
    aTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGAc                                                                                     >  2:4185725/1‑93 (MQ=255)
    aTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGAc                                                                                     >  2:9706125/1‑93 (MQ=255)
    aTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGAc                                                                                     <  2:9313578/93‑1 (MQ=255)
    aTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGAc                                                                                     >  1:2406930/1‑93 (MQ=255)
    aTTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGAc                                                                                     >  1:633992/1‑93 (MQ=255)
      ttAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAg                                                                                   >  1:9139169/1‑93 (MQ=255)
       tAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGc                                                                                  <  2:4740621/93‑1 (MQ=255)
       tAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGc                                                                                  <  1:2970339/93‑1 (MQ=255)
          tGATCAACAATAATGTATGACCATGAATTCTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGc                                                                               <  2:4029635/93‑1 (MQ=255)
          tGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGc                                                                               <  1:6731970/93‑1 (MQ=255)
          tGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGc                                                                               >  1:8651420/1‑93 (MQ=255)
            aTCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCct                                                                             >  2:5252267/1‑93 (MQ=255)
            aTCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCct                                                                             >  2:8940530/1‑93 (MQ=255)
             tCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCctc                                                                            <  1:7680673/93‑1 (MQ=255)
             tCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCctc                                                                            <  1:8863468/93‑1 (MQ=255)
              caacaaTAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCctct                                                                           >  2:1084950/1‑93 (MQ=255)
                 caaTAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTaaa                                                                        >  1:7480100/1‑93 (MQ=255)
                 caaTAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTaaa                                                                        <  1:5329150/93‑1 (MQ=255)
                 caaTAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTaaa                                                                        <  2:8117969/93‑1 (MQ=255)
                  aataatGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTaaaa                                                                       <  2:9098351/93‑1 (MQ=255)
                    taatGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATa                                                                     >  1:4239033/1‑93 (MQ=255)
                      atGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGa                                                                   <  1:3390904/93‑1 (MQ=255)
                      atGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGa                                                                   >  2:9900331/1‑93 (MQ=255)
                      aggTATCACCAGGAATACTATGGACTTCATTAAACATAATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGa                                                                   <  2:6295050/91‑1 (MQ=255)
                        gTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAg                                                                 <  1:1842555/93‑1 (MQ=255)
                        gTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAg                                                                 <  2:800334/93‑1 (MQ=255)
                             aCCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCa                                                            >  1:10018273/1‑93 (MQ=255)
                              ccATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAg                                                           >  1:3914478/1‑93 (MQ=255)
                                aTGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGaa                                                         >  2:3283690/1‑93 (MQ=255)
                                aTGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGaa                                                         <  2:1562754/93‑1 (MQ=255)
                                aTGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGaa                                                         <  2:6605299/93‑1 (MQ=255)
                                aTGAATACTATGGACTACATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGaa                                                         >  2:6709048/1‑93 (MQ=255)
                                  gAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGt                                                       <  1:3154741/93‑1 (MQ=255)
                                  gAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGt                                                       <  1:8283746/93‑1 (MQ=255)
                                     tACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACt                                                    >  1:7815844/1‑93 (MQ=255)
                                      aCTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTg                                                   <  1:3543403/93‑1 (MQ=255)
                                        tATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGAc                                                 <  1:6533883/93‑1 (MQ=255)
                                        tATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGAc                                                 <  1:5849889/93‑1 (MQ=255)
                                        tATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGGACTGAc                                                 >  1:5805338/1‑93 (MQ=255)
                                         aTGGACTTCATTAAACATGATGAACCCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACa                                                <  2:9388864/93‑1 (MQ=255)
                                         aTGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACa                                                >  1:4092681/1‑93 (MQ=255)
                                          tGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAg                                               >  2:8463208/1‑93 (MQ=255)
                                           ggACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGa                                              >  1:7197534/1‑93 (MQ=255)
                                           ggACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGa                                              <  2:9473347/93‑1 (MQ=255)
                                            gACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGAt                                             <  2:7295110/93‑1 (MQ=255)
                                             aCTTCATTAAACATGATGAACCCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATg                                            <  2:4850965/93‑1 (MQ=255)
                                              cTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGt                                           <  2:7194760/93‑1 (MQ=255)
                                               ttCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGtt                                          >  2:3660684/1‑93 (MQ=255)
                                               ttCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGtt                                          <  1:4456633/93‑1 (MQ=255)
                                               ttCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGtt                                          >  2:1814905/1‑93 (MQ=255)
                                                   ttAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCAc                                      <  1:8318743/93‑1 (MQ=255)
                                                   ttAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCAc                                      >  1:5475792/1‑93 (MQ=255)
                                                     aaaCATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCActg                                    >  2:2289908/1‑93 (MQ=255)
                                                      aaCATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCActgc                                   <  1:2072023/93‑1 (MQ=255)
                                                        cATGATGAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCActgctg                                 >  2:8724884/1‑93 (MQ=255)
                                                         atgatgAAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCtgt                                <  1:2755440/93‑1 (MQ=255)
                                                          tgatgaAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCtgtg                               >  2:7666120/1‑93 (MQ=255)
                                                           gatgaAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCtgtgt                              >  2:4575301/1‑93 (MQ=255)
                                                             tgaAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTgtgtgt                            >  1:5079319/1‑93 (MQ=255)
                                                              gaAACCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTgta                           <  1:2589552/93‑1 (MQ=255)
                                                               aaaCCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTgtag                          >  1:1375851/1‑93 (MQ=255)
                                                                aaCCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTgtagt                         >  1:8500281/1‑93 (MQ=255)
                                                                 aCCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTgtagta                        >  2:3415472/1‑93 (MQ=255)
                                                                 aCCCCGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTgtagta                        >  2:5396530/1‑93 (MQ=255)
                                                                  ccccGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGtaa                       >  2:2495153/1‑93 (MQ=255)
                                                                   cccGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGtaat                      >  2:10036017/1‑93 (MQ=255)
                                                                   cccGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGtaat                      <  1:969565/93‑1 (MQ=255)
                                                                   cccGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGtaat                      <  2:2708699/93‑1 (MQ=255)
                                                                    ccGCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGtaata                     >  1:1955322/1‑93 (MQ=255)
                                                                      gCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATaaa                   <  2:1688633/93‑1 (MQ=255)
                                                                      gCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATaaa                   >  2:6532141/1‑93 (MQ=255)
                                                                      gCTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATaaa                   >  1:4562107/1‑93 (MQ=255)
                                                                       cTTTTTTTATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAAt                  >  1:57313/1‑93 (MQ=255)
                                                                        tttttttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATc                 <  2:8359818/93‑1 (MQ=255)
                                                                        tttttttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATc                 <  1:4905409/93‑1 (MQ=255)
                                                                        tttttttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATc                 <  1:4752818/93‑1 (MQ=255)
                                                                        tttttttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATc                 <  2:4312444/93‑1 (MQ=255)
                                                                         ttttttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCa                >  1:1110697/1‑93 (MQ=255)
                                                                         ttttttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCa                >  1:6086416/1‑93 (MQ=255)
                                                                         ttttttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCa                >  1:6278184/1‑93 (MQ=255)
                                                                           ttttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGc              >  2:5000559/1‑93 (MQ=255)
                                                                           ttttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGc              <  2:9583325/93‑1 (MQ=255)
                                                                            tttATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGcc             <  2:6069210/93‑1 (MQ=255)
                                                                              tATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCt           >  2:3786387/1‑93 (MQ=255)
                                                                              tATTGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCt           <  1:6042062/93‑1 (MQ=255)
                                                                                ttGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTaa         <  2:1834057/93‑1 (MQ=255)
                                                                                ttGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTaa         >  1:7059555/1‑93 (MQ=255)
                                                                                ttGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTaa         <  1:6924037/93‑1 (MQ=255)
                                                                                 tGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTaaa        <  2:6949550/93‑1 (MQ=255)
                                                                                 tGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTaaa        >  2:1300278/1‑93 (MQ=255)
                                                                                 tGATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTaaa        <  1:6236017/93‑1 (MQ=255)
                                                                                  gATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAAt       >  2:1152415/1‑93 (MQ=255)
                                                                                  gATCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAAt       >  1:932261/1‑93 (MQ=255)
                                                                                    tCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATgg     >  1:8766906/1‑93 (MQ=255)
                                                                                    tCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATgg     <  1:2258542/93‑1 (MQ=255)
                                                                                    tCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATgg     <  1:8117583/93‑1 (MQ=255)
                                                                                    tCGCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATgg     <  2:5651812/93‑1 (MQ=255)
                                                                                      gCACGCCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATGGGt   <  2:2821830/93‑1 (MQ=255)
                                                                                          |                                                                                         
GAAGATTTAGTGATCAACAATAATGTATGACCATGAATACTATGGACTTCATTAAACATGATGAAACCCCGCTTTTTTTATTGATCGCACACCTGACAGCTGCCTCTAAAATAGAAGCACCAGAAGTACTGACAGATGTTGCACTGCTGTGTGTAGTAATAAATCAGCCCTAAATGGGTA  >  NC_000913/1102361‑1102540

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