Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,284,989 (C)5→4 coding (1932/2511 nt) yejO ← putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,284,9850C.100.0% 478.4 / NA 105coding (1936/2511 nt)yejOputative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (51/54);  total (51/54)

Only 100 of 117 total aligned reads displayed.

ACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGAATT  >  NC_000913/2284893‑2285073
                                                                                            |                                                                                        
aCTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGc                                                                                          >  1:4162787/1‑93 (MQ=255)
  tCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGccc                                                                                        >  1:3707434/1‑93 (MQ=255)
  ccATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGccc                                                                                        <  2:2986798/92‑1 (MQ=255)
   cATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTGTGATGCTCCCAACCGGCATAAGcccc                                                                                       >  1:2079972/1‑93 (MQ=255)
    aTCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGcccca                                                                                      >  1:3546991/1‑92 (MQ=255)
     tCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGccccag                                                                                     <  1:887537/93‑3 (MQ=255)
     tCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGccccag                                                                                     <  1:2547336/93‑3 (MQ=255)
       ttGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGccccaggg                                                                                   <  2:2100455/93‑5 (MQ=255)
       ttGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGccccaggg                                                                                   >  1:198978/1‑89 (MQ=255)
        tGCCATGGATGGTTTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGt                                                                                 <  2:1346777/93‑1 (MQ=255)
        tGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGt                                                                                 >  1:4617016/1‑93 (MQ=255)
         tccATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGta                                                                                <  1:468195/92‑1 (MQ=255)
         gCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGta                                                                                <  2:4288724/93‑1 (MQ=255)
          ccATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGtat                                                                               >  1:2490526/1‑93 (MQ=255)
           cATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGtata                                                                              <  2:215630/93‑1 (MQ=255)
           cATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGtata                                                                              <  1:342311/93‑1 (MQ=255)
            atggatggTGTTGGCAAAACGTTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAg                                                                             >  1:1083941/1‑93 (MQ=255)
            atggatggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAg                                                                             <  2:332788/93‑1 (MQ=255)
            atggatggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAg                                                                             >  2:2063018/1‑93 (MQ=255)
             tggatggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGc                                                                            >  1:4186377/1‑93 (MQ=255)
             tggatggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGc                                                                            <  1:2867029/93‑1 (MQ=255)
             tagatggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGc                                                                            >  2:2703791/3‑93 (MQ=255)
               gatggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTa                                                                          <  1:4112618/93‑1 (MQ=255)
                atggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTAt                                                                         <  2:2431572/93‑1 (MQ=255)
                atggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTAt                                                                         <  1:2213995/93‑1 (MQ=255)
                  ggTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCg                                                                       >  1:1765894/1‑93 (MQ=255)
                     gtTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGata                                                                    <  1:1067961/93‑1 (MQ=255)
                      ttGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGatat                                                                   >  1:3592684/1‑93 (MQ=255)
                       tGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATAtt                                                                  >  2:1503289/1‑93 (MQ=255)
                       tGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATAtt                                                                  <  2:1202424/93‑1 (MQ=255)
                        ggCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTa                                                                 >  1:3450323/1‑93 (MQ=255)
                        ggCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTa                                                                 <  1:2272525/93‑1 (MQ=255)
                           aaaaCGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCt                                                              >  1:2730266/1‑93 (MQ=255)
                            aaaCGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCtt                                                             <  1:826648/93‑1 (MQ=255)
                            aaaCGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCtt                                                             <  2:4112601/93‑1 (MQ=255)
                             aaCGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCttt                                                            <  2:3747089/93‑1 (MQ=255)
                              aCGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTg                                                           >  1:1133761/1‑93 (MQ=255)
                              aCGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTg                                                           >  1:213901/1‑93 (MQ=255)
                              aCGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTg                                                           >  2:1919761/1‑93 (MQ=255)
                              aCGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTg                                                           <  1:3895201/93‑1 (MQ=255)
                                 gTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTgccg                                                        >  1:2692615/1‑93 (MQ=255)
                                 gTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGCATAGCTATCGATATTACCTTTgccg                                                        >  1:1547042/1‑93 (MQ=255)
                                  tCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTgccgc                                                       >  2:4251246/1‑93 (MQ=255)
                                  tCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTgccgc                                                       <  1:3904319/93‑1 (MQ=255)
                                    aaCTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGccgccg                                                     >  2:4047505/1‑93 (MQ=255)
                                      cTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATAATACCTTTGCCGCcgcg                                                   >  2:4173769/1‑93 (MQ=255)
                                       tttCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGa                                                  <  1:3809412/93‑1 (MQ=255)
                                        ttCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGAt                                                 >  1:2640215/1‑93 (MQ=255)
                                             caccCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCaaaa                                            <  2:3861807/93‑1 (MQ=255)
                                             caccCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCaaaa                                            <  1:4468389/93‑1 (MQ=255)
                                              accCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAAc                                           <  2:3520435/93‑1 (MQ=255)
                                              accCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAAc                                           <  1:2641213/93‑1 (MQ=255)
                                               ccccATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAAcc                                          <  2:479494/93‑1 (MQ=255)
                                                cccATGAACATAGGCACCGTTCTGATGCTCCTAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCa                                         >  2:487903/1‑93 (MQ=255)
                                                 ccATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCaa                                        <  2:1930654/93‑1 (MQ=255)
                                                  cATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAat                                       <  1:123239/93‑1 (MQ=255)
                                                  cATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAat                                       <  1:1679278/93‑1 (MQ=255)
                                                   aTCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAata                                      >  2:149573/1‑93 (MQ=255)
                                                    tCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAatat                                     >  1:1790109/1‑93 (MQ=255)
                                                     cAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATc                                    >  1:2029199/1‑93 (MQ=255)
                                                      aaCATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGAATAGCTATCGATATTACCTTTGCGGCCGCGATCAAAACCAATATCa                                   >  2:1047881/1‑93 (MQ=255)
                                                       aCATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAg                                  <  1:1870331/93‑1 (MQ=255)
                                                        cATAGGCACCGTTCTGAGGCTCCCAAAGGGCATAAGGCCCCGGGTAGAGCTATCGTTATTACCGGTGCCGACGCGGCAAAAAGCAATAGCAGa                                 >  2:3904048/1‑93 (MQ=255)
                                                         aTAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGaa                                <  1:149412/93‑1 (MQ=255)
                                                          tAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAAt                               >  2:1813568/1‑93 (MQ=255)
                                                          tAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAAt                               >  2:4057785/1‑93 (MQ=255)
                                                           aGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATg                              >  1:3333662/1‑93 (MQ=255)
                                                             gCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATgag                            <  2:3455256/93‑1 (MQ=255)
                                                              cACCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATgaga                           <  1:3925981/93‑1 (MQ=255)
                                                               aCCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGagag                          <  1:3318269/93‑1 (MQ=255)
                                                               aCCGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGagag                          <  2:1072307/93‑1 (MQ=255)
                                                                ccGTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGagaga                         >  1:1055922/1‑92 (MQ=255)
                                                                  gTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTaa                       <  1:1075235/93‑1 (MQ=255)
                                                                  gTTCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTaa                       <  1:599047/93‑1 (MQ=255)
                                                                   ttCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAAc                      <  1:1578054/93‑1 (MQ=255)
                                                                   ttCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAAc                      <  1:3350538/93‑1 (MQ=255)
                                                                   ttCTGATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAAc                      >  1:2314836/1‑93 (MQ=255)
                                                                       gATGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCaaa                  >  1:739294/1‑93 (MQ=255)
                                                                        aTGCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAg                 >  1:4238732/1‑93 (MQ=255)
                                                                          gCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGAt               >  2:387264/1‑93 (MQ=255)
                                                                          gCTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGAt               <  1:3642504/93‑1 (MQ=255)
                                                                           cTCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATc              <  2:4040976/93‑1 (MQ=255)
                                                                            tCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCa             >  2:3702362/1‑93 (MQ=255)
                                                                            tCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCa             >  2:3034995/1‑93 (MQ=255)
                                                                            tCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCa             >  1:1664802/1‑93 (MQ=255)
                                                                            tCCCAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCa             >  1:1877777/1‑93 (MQ=255)
                                                                              ccAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAg           >  2:4050141/1‑93 (MQ=255)
                                                                              ccAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAg           >  2:3353589/1‑93 (MQ=255)
                                                                              ccAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAg           >  1:1470907/1‑93 (MQ=255)
                                                                               cAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGc          <  2:3391888/93‑1 (MQ=255)
                                                                               cAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGc          >  2:3385558/1‑93 (MQ=255)
                                                                               cAACCGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGc          >  2:3329108/1‑93 (MQ=255)
                                                                                  ccGGCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCcgc       >  2:2053109/1‑93 (MQ=255)
                                                                                   cGGCATAAGCCCCAGGGTACAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCcgcg      >  2:4584657/1‑93 (MQ=255)
                                                                                     gCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGaa    <  2:149088/93‑1 (MQ=255)
                                                                                     gCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGaa    <  1:1957420/93‑1 (MQ=255)
                                                                                     gCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGaa    <  1:2913092/93‑1 (MQ=255)
                                                                                     gCATAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGaa    >  1:3160434/1‑93 (MQ=255)
                                                                                       aTAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGAAtt  >  2:2981828/1‑93 (MQ=255)
                                                                                       aTAAGCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCACGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGAAtt  >  2:2926928/1‑93 (MQ=255)
                                                                                            |                                                                                        
ACTCATCTTGCCATGGATGGTGTTGGCAAAACGGTCAACTTTCACCACCCCATCAACATAGGCACCGTTCTGATGCTCCCAACCGGCATAAGCCCCCAGGGTATAGCTATCGATATTACCTTTGCCGCCGCGATCAAAACCAATATCAGAATGAGAGTAACCAAAGATCAAGCCGCGAATT  >  NC_000913/2284893‑2285073

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