Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,356,994 C→T A206T (GCA→ACA)  cadB ← transport of lysine/cadaverine

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,356,9940CT100.0% 1086.8 / NA 291A206T (GCA→ACA) cadBtransport of lysine/cadaverine
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (157/134);  total (157/134)

Only 100 of 296 total aligned reads displayed.

GATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAG  >  NC_000913/4356902‑4357083
                                                                                            |                                                                                         
gATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGCACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGc                                                                                           >  2:3224288/1‑93 (MQ=255)
  tGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTgg                                                                                         <  2:3757619/93‑1 (MQ=255)
   gTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGa                                                                                        >  1:3819265/1‑93 (MQ=255)
   gTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGa                                                                                        <  1:2879215/93‑1 (MQ=255)
    tAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGAt                                                                                       <  2:2095635/93‑1 (MQ=255)
    tAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGAt                                                                                       >  2:6044668/1‑93 (MQ=255)
      aaCAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTc                                                                                     >  1:5614804/1‑93 (MQ=255)
      aaCAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTc                                                                                     >  2:2825784/1‑93 (MQ=255)
      aaCAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTc                                                                                     >  1:3416219/1‑93 (MQ=255)
       aCAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCa                                                                                    >  2:6466089/1‑93 (MQ=255)
          aTACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAaca                                                                                 <  1:3888450/93‑1 (MQ=255)
           tACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAacac                                                                                <  1:1019423/93‑1 (MQ=255)
           tACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAacac                                                                                >  2:5089550/1‑93 (MQ=255)
           tACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAacac                                                                                <  1:4290249/93‑1 (MQ=255)
            aCCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACAcc                                                                               >  1:4806180/1‑93 (MQ=255)
            aCCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACAcc                                                                               <  2:4001214/93‑1 (MQ=255)
             ccTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACAccc                                                                              <  1:1428768/93‑1 (MQ=255)
              cTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCa                                                                             <  2:6304001/93‑1 (MQ=255)
              cTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCa                                                                             <  2:5239326/93‑1 (MQ=255)
               tGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCAc                                                                            >  1:4762713/1‑93 (MQ=255)
               tGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCAc                                                                            >  1:1263125/1‑93 (MQ=255)
               cgCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCAc                                                                            <  2:2239650/92‑1 (MQ=255)
                  tAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGaa                                                                         >  1:5787623/1‑93 (MQ=255)
                   aaaCCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAg                                                                        >  2:5710914/1‑93 (MQ=255)
                   aaaCCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAg                                                                        >  1:3019951/1‑93 (MQ=255)
                   aaaCCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAg                                                                        <  1:6471082/93‑1 (MQ=255)
                     aCCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGc                                                                      >  2:3480362/1‑93 (MQ=255)
                     aCCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGc                                                                      >  1:938212/1‑93 (MQ=255)
                     aCCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGc                                                                      >  2:3825509/1‑93 (MQ=255)
                       cAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGccc                                                                    >  2:3426254/1‑93 (MQ=255)
                       cAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGccc                                                                    >  2:2028834/1‑93 (MQ=255)
                       cAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGccc                                                                    >  2:3163339/1‑93 (MQ=255)
                        aGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCca                                                                   >  1:4570973/1‑93 (MQ=255)
                         gTACCCAGCATGGTTACCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCcac                                                                  >  1:1976944/1‑93 (MQ=255)
                          tACCCGGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTTACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCcaca                                                                 <  2:243490/93‑1 (MQ=255)
                          tACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATTCCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCcaca                                                                 <  2:656376/93‑1 (MQ=255)
                          tACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATATCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCcaca                                                                 <  2:3558503/93‑1 (MQ=255)
                          tACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCcaca                                                                 <  2:3618308/93‑1 (MQ=255)
                           tcccAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAg                                                                <  2:805514/92‑1 (MQ=255)
                           aCCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAg                                                                >  1:6483975/1‑93 (MQ=255)
                           aCCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAg                                                                >  1:75704/1‑93 (MQ=255)
                            cccAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAgg                                                               >  2:6113532/1‑93 (MQ=255)
                             ccAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgc                                                              >  1:1692451/1‑93 (MQ=255)
                                gCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgcaga                                                           <  2:6661600/93‑1 (MQ=255)
                                   tGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgcagagca                                                        >  2:1448948/1‑93 (MQ=255)
                                   tGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgcagagca                                                        >  1:6146371/1‑93 (MQ=255)
                                     gTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGgcagagcaga                                                      <  2:2364514/93‑1 (MQ=255)
                                       tGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAAt                                                    <  1:537153/93‑1 (MQ=255)
                                         ccAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATAc                                                  <  2:5411922/93‑1 (MQ=255)
                                           aGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACtt                                                >  2:3287027/1‑93 (MQ=255)
                                           aGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACtt                                                >  2:1562336/1‑93 (MQ=255)
                                           aGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACtt                                                >  2:3648940/1‑93 (MQ=255)
                                            gCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGCGCAGAATACatt                                               >  1:5996084/1‑93 (MQ=255)
                                            gCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACttt                                               <  2:3124435/93‑1 (MQ=255)
                                            gCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACttt                                               >  1:2765341/1‑93 (MQ=255)
                                            gCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACttt                                               >  2:1876880/1‑93 (MQ=255)
                                            gCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACttt                                               <  1:1602330/93‑1 (MQ=255)
                                             cGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACtttt                                              >  1:2565948/1‑93 (MQ=255)
                                               gAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTaa                                            <  1:3965114/93‑1 (MQ=255)
                                                aaCGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAAt                                           <  2:1324925/93‑1 (MQ=255)
                                                aaCGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAAt                                           <  1:4950854/93‑1 (MQ=255)
                                                aaCGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAAt                                           <  1:6349879/93‑1 (MQ=255)
                                                 aCGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATg                                          <  1:2804616/93‑1 (MQ=255)
                                                   ggNACGTTTCGGGTTTTTAACCNTACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGAt                                        <  2:4951190/93‑1 (MQ=255)
                                                     tACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCg                                      <  1:3288823/93‑1 (MQ=255)
                                                     tACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCg                                      >  1:856937/1‑93 (MQ=255)
                                                     tACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCg                                      >  2:1138333/1‑93 (MQ=255)
                                                      aCGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGc                                     >  2:4595476/1‑93 (MQ=255)
                                                         tttCGGGTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATg                                  >  1:2249200/1‑93 (MQ=255)
                                                               gTTTTTAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATc                            >  1:4648888/1‑93 (MQ=255)
                                                                tttttAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCa                           <  1:3518166/93‑1 (MQ=255)
                                                                tttttAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCa                           <  1:4679853/93‑1 (MQ=255)
                                                                tttttAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCa                           >  1:5267708/1‑93 (MQ=255)
                                                                tttttAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGACCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCa                           >  2:2953736/1‑93 (MQ=255)
                                                                gtgttAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATTCTTTTAATGATCGCATGACCATCa                           <  2:5533354/90‑1 (MQ=255)
                                                                  tttAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGt                         <  2:1943639/93‑1 (MQ=255)
                                                                    tAACCATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTgg                       <  2:4747313/93‑1 (MQ=255)
                                                                        cATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATc                   >  2:3934663/1‑93 (MQ=255)
                                                                        cATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATc                   >  2:4169462/1‑93 (MQ=255)
                                                                        cATACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATc                   >  2:452359/1‑93 (MQ=255)
                                                                         aTACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATcc                  >  2:4411542/1‑93 (MQ=255)
                                                                         aTACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATcc                  <  1:5743471/93‑1 (MQ=255)
                                                                         aTACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATcc                  >  1:2410245/1‑93 (MQ=255)
                                                                          tACCAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCg                 >  1:4131788/1‑93 (MQ=255)
                                                                            ccAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCa               >  2:1025365/1‑93 (MQ=255)
                                                                            ccAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCa               >  1:3601006/1‑93 (MQ=255)
                                                                             cAGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAg              <  2:4208537/93‑1 (MQ=255)
                                                                              aGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGt             <  1:5077775/93‑1 (MQ=255)
                                                                              aGTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGt             <  2:5605380/93‑1 (MQ=255)
                                                                               gTACTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTa            >  1:6691355/1‑93 (MQ=255)
                                                                                 aCTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTAtt          >  2:6542217/1‑93 (MQ=255)
                                                                                 aCTTACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTAtt          >  2:277598/1‑93 (MQ=255)
                                                                                   ttACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTcc        >  1:5168893/1‑93 (MQ=255)
                                                                                   ttACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTcc        >  1:6101020/1‑93 (MQ=255)
                                                                                   ttACAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTcc        >  1:2157439/1‑93 (MQ=255)
                                                                                     aCAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAg      <  2:886482/93‑1 (MQ=255)
                                                                                      cAGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGt     >  1:1379081/1‑93 (MQ=255)
                                                                                       aGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGtt    <  2:1245479/93‑1 (MQ=255)
                                                                                       aGCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGtt    <  2:1693102/93‑1 (MQ=255)
                                                                                        gCTGTGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTa   <  2:696648/93‑1 (MQ=255)
                                                                                            |                                                                                         
GATGTAAACAATACCTGCTAAACCAGTACCCAGCATGGTTGCCAGCGGAACGGTACGTTTCGGGTTTTTAACCATACCAGTACTTACAGCTGCGGATTCAACACCCACGAAGGCCCACAGGCAGAGCAGAATACTTTTAATGATCGCATGACCATCAGTGGTATCCGCAGTATTCCAGTTAG  >  NC_000913/4356902‑4357083

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