Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,227,150 G→C P191R (CCG→CGG)  pepE ← peptidase E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,227,1500GC100.0% 953.2 / NA 292P191R (CCG→CGG) pepEpeptidase E
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (152/140);  total (152/140)

Only 100 of 294 total aligned reads displayed.

ACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGC  >  NC_000913/4227060‑4227239
                                                                                          |                                                                                         
aCCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGt                                                                                         <  1:292278/93‑1 (MQ=255)
aCCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGt                                                                                         <  2:1619381/93‑1 (MQ=255)
aCCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGt                                                                                         >  1:1656336/1‑93 (MQ=255)
 ccGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTa                                                                                        <  1:3857746/93‑1 (MQ=255)
  cGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAg                                                                                       <  2:1201287/93‑1 (MQ=255)
  cGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAg                                                                                       <  2:4815862/93‑1 (MQ=255)
   gCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGa                                                                                      <  2:2233622/93‑1 (MQ=255)
    cttcttCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGAc                                                                                     <  1:5190832/93‑1 (MQ=255)
    cttcttCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGAc                                                                                     <  2:3779160/93‑1 (MQ=255)
     ttcttcACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGAcc                                                                                    <  1:4189735/93‑1 (MQ=255)
     ttcttcACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGAcc                                                                                    >  1:4656039/1‑93 (MQ=255)
       cttcACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCaa                                                                                  <  1:2541709/93‑1 (MQ=255)
       cttcACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCaa                                                                                  <  1:5026616/93‑1 (MQ=255)
        ttcACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCaat                                                                                 >  1:5082042/1‑93 (MQ=255)
             cAGCCTTAAACACATAAGTGGTGTTCGGTCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCg                                                                            <  1:6059697/93‑1 (MQ=255)
             cAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCg                                                                            >  1:1309030/1‑93 (MQ=255)
               gCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTc                                                                          <  1:2681253/93‑1 (MQ=255)
                ccTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCa                                                                         <  2:3444494/93‑1 (MQ=255)
                ccTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCa                                                                         <  2:2645658/93‑1 (MQ=255)
                 cTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGCCCATTAATCGTCAg                                                                        >  2:5245540/1‑93 (MQ=255)
                 cTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAg                                                                        <  1:5750093/93‑1 (MQ=255)
                 cTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAg                                                                        >  1:3959764/1‑93 (MQ=255)
                 cTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAg                                                                        <  1:2070414/93‑1 (MQ=255)
                 cTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAg                                                                        <  1:700020/93‑1 (MQ=255)
                     aaCACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCt                                                                    <  1:2182104/93‑1 (MQ=255)
                     aaCACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCt                                                                    <  1:197847/93‑1 (MQ=255)
                        acaTAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGgc                                                                 >  1:1018981/1‑93 (MQ=255)
                        acaTAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGgc                                                                 >  1:401749/1‑93 (MQ=255)
                         caTAAGTGATGTTCGGGCCTCCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGCTCTGgcg                                                                >  2:6120126/1‑93 (MQ=255)
                          aTAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTCCCTTCCCGTAGACCAATAATCGTCAGTTCTGgcgc                                                               <  1:1039609/93‑1 (MQ=255)
                           tAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGcgcg                                                              >  1:2337530/1‑93 (MQ=255)
                           tAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGcgcg                                                              >  1:1802499/1‑93 (MQ=255)
                            aaGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGcga                                                             <  1:3602894/93‑1 (MQ=255)
                             aGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGcgac                                                            <  2:5462540/93‑1 (MQ=255)
                               tggtgTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGCTCCAGTTCCCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGcgacga                                                          <  1:6309898/93‑1 (MQ=255)
                                  tgtTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACca                                                       <  2:2945998/93‑1 (MQ=255)
                                   gtTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACcag                                                      <  2:1080616/93‑1 (MQ=255)
                                   gtTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACcag                                                      <  2:4145004/93‑1 (MQ=255)
                                    ttCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACcagc                                                     >  2:1349951/1‑93 (MQ=255)
                                    ttCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACcagc                                                     <  2:1450168/93‑1 (MQ=255)
                                    ttCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACcagc                                                     >  1:5876603/1‑93 (MQ=255)
                                     tCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACcagca                                                    >  1:883770/1‑93 (MQ=255)
                                     tCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACcagca                                                    <  1:1185079/93‑1 (MQ=255)
                                       gggCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGt                                                  >  1:323670/1‑93 (MQ=255)
                                       gggCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGt                                                  <  2:1510141/93‑1 (MQ=255)
                                         gCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTc                                                <  2:3075769/93‑1 (MQ=255)
                                         gCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTc                                                <  1:2840901/93‑1 (MQ=255)
                                            aCCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTcgcg                                             >  2:3694670/1‑93 (MQ=255)
                                             cccAGCGCAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGa                                            >  2:1694381/1‑93 (MQ=255)
                                              ccAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGaa                                           <  2:5748779/93‑1 (MQ=255)
                                                agcgtagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATa                                         <  2:1320747/93‑1 (MQ=255)
                                                agcgtagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATa                                         >  1:3708725/1‑93 (MQ=255)
                                                 gcgtagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATAc                                        >  1:4819944/1‑93 (MQ=255)
                                                   gtagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACgc                                      <  1:3576401/93‑1 (MQ=255)
                                                   gtagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACgc                                      >  2:1283364/1‑93 (MQ=255)
                                                    tagcgtGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACgct                                     >  1:411121/1‑93 (MQ=255)
                                                         tGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCa                                >  2:3830400/1‑93 (MQ=255)
                                                         tGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCa                                >  2:981280/1‑93 (MQ=255)
                                                         tGACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCa                                >  1:6183811/1‑93 (MQ=255)
                                                          gACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCCGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCAc                               <  1:4441337/93‑1 (MQ=255)
                                                          gACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCAc                               >  2:4895685/1‑93 (MQ=255)
                                                          gACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCAc                               >  2:4321066/1‑93 (MQ=255)
                                                          gACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCAc                               <  2:3967033/93‑1 (MQ=255)
                                                          gACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCAc                               <  2:3227711/93‑1 (MQ=255)
                                                          gACCTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCAc                               <  1:4543192/93‑1 (MQ=255)
                                                            ccTTTACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACgg                             <  2:1794389/93‑1 (MQ=255)
                                                                tACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGttt                         >  2:1041966/1‑93 (MQ=255)
                                                                tACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGttt                         <  1:6130668/93‑1 (MQ=255)
                                                                tACTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGttt                         >  2:4364052/1‑93 (MQ=255)
                                                                 aCTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTc                        >  1:1553311/1‑93 (MQ=255)
                                                                 aCTCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTc                        <  2:4500695/93‑1 (MQ=255)
                                                                   tCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCAc                      >  1:3630511/1‑93 (MQ=255)
                                                                   tCACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCAc                      <  2:768316/93‑1 (MQ=255)
                                                                    cACTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCAcc                     <  1:1742161/93‑1 (MQ=255)
                                                                     aCTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCt                    >  2:140485/1‑93 (MQ=255)
                                                                     aCTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCt                    <  1:3144012/93‑1 (MQ=255)
                                                                     aCTGTGATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCt                    >  2:3478616/1‑93 (MQ=255)
                                                                       tgtgATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCttt                  >  2:1716222/1‑93 (MQ=255)
                                                                        ttgATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTa                 <  2:5584664/92‑1 (MQ=255)
                                                                         tgATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTAt                >  1:3083405/1‑93 (MQ=255)
                                                                         tgATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTAt                <  1:1507307/93‑1 (MQ=255)
                                                                          gATCCAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATg               <  1:1352084/93‑1 (MQ=255)
                                                                             ccAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGcc            >  2:2769290/1‑93 (MQ=255)
                                                                              cAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCt           >  2:6017824/1‑93 (MQ=255)
                                                                              cAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCt           >  2:1590453/1‑93 (MQ=255)
                                                                              cAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCt           >  1:1825646/1‑93 (MQ=255)
                                                                              cAGTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCt           >  1:1145332/1‑93 (MQ=255)
                                                                                gTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGTCTTc         >  1:3484736/1‑93 (MQ=255)
                                                                                gTTACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTc         >  2:648984/1‑93 (MQ=255)
                                                                                 ttACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTcc        >  2:2315599/1‑93 (MQ=255)
                                                                                 ttACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTcc        <  1:5345123/93‑1 (MQ=255)
                                                                                  tACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCg       <  1:758679/93‑1 (MQ=255)
                                                                                  tACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCg       >  2:1970306/1‑93 (MQ=255)
                                                                                  tACCTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCg       <  1:5571409/93‑1 (MQ=255)
                                                                                    ccTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGc     >  2:4192619/1‑93 (MQ=255)
                                                                                    ccTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGc     <  1:3413389/93‑1 (MQ=255)
                                                                                    ccTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGc     >  1:6082378/1‑93 (MQ=255)
                                                                                     cTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCa    <  2:5025054/93‑1 (MQ=255)
                                                                                     cTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCa    <  1:228637/93‑1 (MQ=255)
                                                                                     cTTCCCGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCa    <  1:4262309/93‑1 (MQ=255)
                                                                                          |                                                                                         
ACCGCTTCTTCACCAGCCTTAAACACATAAGTGGTGTTCGGGCCACCCAGCGTAGCGTGACCTTTACTCACTGTGATCCAGTTACCTTCCGGTAGACCAATAATCGTCAGTTCTGGCGCGACGACCAGCAGTTCGCGAATACGCTGCTCACGGGTTTCACCTTTATGGCCTTCCGGCAGC  >  NC_000913/4227060‑4227239

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