Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,702,780 T→C T214A (ACT→GCT)  dppB ← dipeptide transport system permease protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,702,7800TC99.6% 804.8 / ‑5.3 256T214A (ACT→GCT) dppBdipeptide transport system permease protein 1
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (113/141);  total (115/141)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 260 total aligned reads displayed.

CTTTGGCGCGCGCGGTGCGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGTACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTCGATTAGCATAA  >  NC_000913/3702690‑3702869
                                                                                          |                                                                                         
cTTTGGCGCGCGCGGTGCGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCAc                                                                                         >  1:1182708/1‑93 (MQ=255)
cTTTGGCGCGCGCGGTGCGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCAc                                                                                         >  1:1591394/1‑93 (MQ=255)
 tttGGCGCGCGCGGTGCGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCAcc                                                                                        <  1:1832885/93‑1 (MQ=255)
 tttGGCGCGCGCGGTGCGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCAcc                                                                                        >  1:487478/1‑93 (MQ=255)
   tGGCGCGCGCGGTGCGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCa                                                                                      <  1:41987/93‑1 (MQ=255)
    ggCGCGCGCGGTGCGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAg                                                                                     >  1:353073/1‑93 (MQ=255)
    ggCGCGCGCGGTGCGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAg                                                                                     <  2:4192866/93‑1 (MQ=255)
    ggCGCGCGCGGGGCGGATGTAATCCTCGCCCAGCACTTACAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAg                                                                                     >  2:734816/1‑93 (MQ=255)
      cgcgcgcgGTGCGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGca                                                                                   <  2:491772/93‑1 (MQ=255)
          cgcgGTGCGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAAt                                                                               >  1:3116844/1‑93 (MQ=255)
          cgcgGTGCGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAAt                                                                               >  2:1865496/1‑93 (MQ=255)
           gcgGTGCGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATg                                                                              <  2:2821548/93‑1 (MQ=255)
           gcgGTGCGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATg                                                                              <  2:1915459/93‑1 (MQ=255)
            cgGTGCGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATgg                                                                             <  2:2311263/93‑1 (MQ=255)
             ggTGCGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGACTAGCACCCAGCACAATggc                                                                            >  1:1333965/1‑93 (MQ=255)
             ggTGCGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATggc                                                                            <  1:2757081/93‑1 (MQ=255)
               tGCGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATggcag                                                                          <  2:4149793/93‑1 (MQ=255)
               tGCGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGAACCCAGCACAATggcag                                                                          >  1:654408/1‑93 (MQ=255)
                 cGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATggcaggc                                                                        <  1:2374355/93‑1 (MQ=255)
                   gATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCaa                                                                      <  1:950706/93‑1 (MQ=255)
                     tGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGa                                                                    <  2:4171350/93‑1 (MQ=255)
                     tGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGa                                                                    <  1:2362421/93‑1 (MQ=255)
                     tGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGa                                                                    <  1:105665/93‑1 (MQ=255)
                        aaTCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCa                                                                 <  2:4552327/93‑1 (MQ=255)
                        aaTCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCa                                                                 <  2:1596162/93‑1 (MQ=255)
                        aaTCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCa                                                                 <  1:2297264/93‑1 (MQ=255)
                          tCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCata                                                               >  2:1955567/1‑93 (MQ=255)
                            cTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATg                                                             >  2:3352335/1‑93 (MQ=255)
                            cTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATg                                                             <  2:2053275/93‑1 (MQ=255)
                            cTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATg                                                             >  2:1003389/1‑93 (MQ=255)
                             tCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATgg                                                            <  1:379090/93‑1 (MQ=255)
                              cGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAACGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATggg                                                           <  1:317655/93‑1 (MQ=255)
                               gCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATgggg                                                          >  2:585043/1‑92 (MQ=255)
                                  cAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGAc                                                       <  2:4630225/93‑1 (MQ=255)
                                  cAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGAc                                                       >  1:3053673/1‑93 (MQ=255)
                                   aGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACg                                                      <  2:240420/93‑1 (MQ=255)
                                   aGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACg                                                      >  1:2623711/1‑93 (MQ=255)
                                    gCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACgg                                                     >  1:2267970/1‑93 (MQ=255)
                                     cACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGc                                                    >  2:954271/1‑93 (MQ=255)
                                     cACTTCAAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGc                                                    >  1:1534386/1‑93 (MQ=255)
                                      aCTTCCAGCATCGAGGAGCGTGTCACACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCa                                                   <  1:1680962/93‑1 (MQ=255)
                                       cTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCAt                                                  <  1:300105/93‑1 (MQ=255)
                                       cTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGAAt                                                  >  1:3877024/1‑93 (MQ=255)
                                        ttCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATc                                                 <  1:4464429/93‑1 (MQ=255)
                                          ccAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGa                                               <  1:3290353/93‑1 (MQ=255)
                                           cAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGAt                                              >  1:4328918/1‑93 (MQ=255)
                                            aGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATa                                             >  1:2175081/1‑93 (MQ=255)
                                            aGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATa                                             <  1:563147/93‑1 (MQ=255)
                                             gCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATaa                                            >  2:2709479/1‑93 (MQ=255)
                                               aTCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAg                                          <  2:468579/93‑1 (MQ=255)
                                               aTCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAg                                          >  2:898511/1‑93 (MQ=255)
                                                tCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGt                                         <  2:3924790/93‑1 (MQ=255)
                                                tCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGt                                         >  1:75229/1‑93 (MQ=255)
                                                tCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGt                                         <  1:2642988/93‑1 (MQ=255)
                                                 cGAGGAGCGTGTCATACGCACAAGGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGtt                                        <  2:1175840/93‑1 (MQ=255)
                                                  gaggagCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTg                                       >  2:1405295/1‑93 (MQ=255)
                                                  gaggagCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTg                                       <  1:3833538/93‑1 (MQ=255)
                                                   aggagCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGc                                      >  2:767438/1‑93 (MQ=255)
                                                   aggagCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGc                                      <  1:1659234/93‑1 (MQ=255)
                                                    ggagCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGcc                                     >  2:3405503/1‑93 (MQ=255)
                                                       gCGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTc                                  >  1:3679977/1‑93 (MQ=255)
                                                        cGTGTCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCt                                 >  1:1447238/1‑93 (MQ=255)
                                                          tgtCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCtta                               >  1:2905373/1‑92 (MQ=255)
                                                          tgtCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTc                               >  2:3789917/1‑93 (MQ=255)
                                                           gtCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCa                              <  2:1256671/93‑1 (MQ=255)
                                                            tCATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCAc                             >  1:4262081/1‑93 (MQ=255)
                                                             cATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCAcc                            >  1:2871113/1‑93 (MQ=255)
                                                             cATACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCAcc                            <  1:2045260/93‑1 (MQ=255)
                                                              aTACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCAccc                           >  2:3532258/1‑93 (MQ=255)
                                                               tACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCAcccc                          <  1:2465225/93‑1 (MQ=255)
                                                               tACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCAcccc                          >  2:655502/1‑93 (MQ=255)
                                                               tACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCAcccc                          <  1:123829/93‑1 (MQ=255)
                                                               tACGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCAcccc                          <  2:3479266/93‑1 (MQ=255)
                                                                aCGCACAATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCa                         <  2:1699355/93‑1 (MQ=255)
                                                                    acaATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATg                     <  2:256454/93‑1 (MQ=255)
                                                                    acaATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATg                     >  2:320460/1‑93 (MQ=255)
                                                                    acaATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATg                     >  1:8808/1‑93 (MQ=255)
                                                                     caATGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATgg                    <  1:2528071/93‑1 (MQ=255)
                                                                       aTGACCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGGGACGGCATCGATAAAGTTGCCGTCTTGACCCCAGATcgcg                  >  2:2945572/1‑93 (MQ=255)
                                                                          aCCGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGgtg               <  1:804841/93‑1 (MQ=255)
                                                                           ccGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGgtgt              <  1:4536477/93‑1 (MQ=255)
                                                                            cGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTc             <  1:2618563/93‑1 (MQ=255)
                                                                            cGCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTc             >  2:2783137/1‑93 (MQ=255)
                                                                             gCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTCg            >  2:3249145/1‑93 (MQ=255)
                                                                             gCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTCg            >  1:1519144/1‑93 (MQ=255)
                                                                             gCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTCg            <  1:3259258/93‑1 (MQ=255)
                                                                             gCCAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTCg            >  2:4255668/1‑93 (MQ=255)
                                                                              caaGCGGAATAGCACCCAGCACAATTGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTCGa           <  2:432967/91‑1 (MQ=255)
                                                                               cAGCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTCGAt          >  2:37131/1‑93 (MQ=255)
                                                                                 gCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTCGATTa        <  2:1793506/93‑1 (MQ=255)
                                                                                 gCGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGAGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTCGATTa        >  1:2859789/1‑93 (MQ=255)
                                                                                  cGGAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTCGATTAg       <  2:3140600/93‑1 (MQ=255)
                                                                                   ggAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGTCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTCGATTAGc      <  1:4017408/93‑1 (MQ=255)
                                                                                   ggAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTCGATTAGc      >  1:3717782/1‑93 (MQ=255)
                                                                                   ggAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTCGATTAGc      <  2:388841/93‑1 (MQ=255)
                                                                                   ggAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTCGATTAGc      >  1:1125903/1‑93 (MQ=255)
                                                                                    gAATAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTCGATTAGCa     >  1:4553295/1‑93 (MQ=255)
                                                                                      aTAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTCGATTAGCATa   >  2:4505158/1‑93 (MQ=255)
                                                                                      aTAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTCGATTAGCATa   <  1:1314138/93‑1 (MQ=255)
                                                                                       tAGCACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGCCTTCACCCCAGATGGTGGTGTCGATTATCATaa  >  2:213792/1‑93 (MQ=255)
                                                                                          |                                                                                         
CTTTGGCGCGCGCGGTGCGGATGTAATCCTCGCCCAGCACTTCCAGCATCGAGGAGCGTGTCATACGCACAATGACCGCCAGCGGAATAGTACCCAGCACAATGGCAGGCAAGATCATATGGGCGACGGCATCGATAAAGTTGCCGTCTTCACCCCAGATGGCGGTGTCGATTAGCATAA  >  NC_000913/3702690‑3702869

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