Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,386,329 C→T R207H (CGC→CAC)  yjeP ← putative periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,386,3290CT75.8% 623.3 / 179.2 311R207H (CGC→CAC) yjePputative periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (44/31);  new base T (120/115);  total (165/146)
Fisher's exact test for biased strand distribution p-value = 2.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.40e-02

Only 100 of 313 total aligned reads displayed.

CCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGC  >  NC_000913/4386237‑4386421
                                                                                            |                                                                                            
    cTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTaa                                                                                          <  1:1067206/93‑1 (MQ=255)
    cTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTaa                                                                                          <  1:2725386/93‑1 (MQ=255)
     tCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGTCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAg                                                                                         <  2:1953021/93‑1 (MQ=255)
     tCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAg                                                                                         <  1:6820855/93‑1 (MQ=255)
      cAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAgc                                                                                        <  1:6317054/93‑1 (MQ=255)
      cAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAgc                                                                                        <  1:907538/93‑1 (MQ=255)
        aGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAgcgc                                                                                      <  2:4364248/93‑1 (MQ=255)
           tGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCt                                                                                   >  1:4421758/1‑93 (MQ=255)
            gACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTa                                                                                  >  2:703081/1‑93 (MQ=255)
             aCGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTaa                                                                                 <  2:6627868/93‑1 (MQ=255)
             aCGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTaa                                                                                 <  2:4026435/93‑1 (MQ=255)
             aCGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTaa                                                                                 <  2:391478/93‑1 (MQ=255)
               gTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAAtt                                                                               >  1:140957/1‑93 (MQ=255)
               gTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAAtt                                                                               >  1:4113217/1‑93 (MQ=255)
                ttGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTc                                                                              <  1:1334387/93‑1 (MQ=255)
                 tGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTcc                                                                             <  1:7159709/93‑1 (MQ=255)
                 tGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTcc                                                                             >  2:2980955/1‑93 (MQ=255)
                 tGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTcc                                                                             <  2:7198661/93‑1 (MQ=255)
                 tGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTcc                                                                             <  1:4515322/93‑1 (MQ=255)
                  ggCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCt                                                                            <  2:1556834/93‑1 (MQ=255)
                  ggCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCt                                                                            <  2:2192954/93‑1 (MQ=255)
                    cTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTgg                                                                          <  2:905635/93‑1 (MQ=255)
                       tttAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCgg                                                                       <  2:3419803/93‑1 (MQ=255)
                       tttAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCgg                                                                       >  1:5334368/1‑93 (MQ=255)
                        ttAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGt                                                                      >  1:2175834/1‑93 (MQ=255)
                        ttAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGt                                                                      >  1:5113169/1‑93 (MQ=255)
                          aaTTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTa                                                                    <  1:3650311/93‑1 (MQ=255)
                          aaTTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTa                                                                    >  1:313058/1‑93 (MQ=255)
                          aaTTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTa                                                                    <  1:1967867/93‑1 (MQ=255)
                          aaTTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTa                                                                    <  1:601769/93‑1 (MQ=255)
                              gATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTgg                                                                <  2:5661635/93‑1 (MQ=255)
                               aTTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGc                                                               <  2:55111/93‑1 (MQ=255)
                                  aCGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGa                                                            <  2:796597/93‑1 (MQ=255)
                                   cGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGAc                                                           <  1:2905681/93‑1 (MQ=255)
                                   cGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGAc                                                           <  2:4997463/93‑1 (MQ=255)
                                    gCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACa                                                          <  2:1413871/93‑1 (MQ=255)
                                    gCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACa                                                          >  2:1718586/1‑93 (MQ=255)
                                    gCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACa                                                          <  2:1789364/93‑1 (MQ=255)
                                     cAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAg                                                         <  1:4279218/93‑1 (MQ=255)
                                        ggCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTg                                                      <  1:7216677/93‑1 (MQ=255)
                                        ggCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTg                                                      >  1:2994697/1‑93 (MQ=255)
                                         gCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTgc                                                     >  1:7106499/1‑93 (MQ=255)
                                           cTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTgcgc                                                   <  1:3570905/93‑1 (MQ=255)
                                            tGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGcc                                                  <  2:4931326/93‑1 (MQ=255)
                                            tGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGcc                                                  >  2:5809732/1‑93 (MQ=255)
                                             gCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCa                                                 >  1:6604608/1‑93 (MQ=255)
                                                aaTACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCt                                              >  2:1039713/1‑93 (MQ=255)
                                                aaTACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCt                                              >  2:2815062/1‑93 (MQ=255)
                                                aaTACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCt                                              >  2:5175896/1‑93 (MQ=255)
                                                  tACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTcc                                            >  1:5805940/1‑93 (MQ=255)
                                                   aCGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCa                                           <  1:1825636/93‑1 (MQ=255)
                                                    cGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAg                                          >  2:5395881/1‑93 (MQ=255)
                                                     gCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGt                                         <  1:3258427/93‑1 (MQ=255)
                                                     gCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGt                                         >  2:666044/1‑93 (MQ=255)
                                                      cATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGtt                                        <  2:73352/93‑1 (MQ=255)
                                                        tCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTtca                                      >  1:4624506/1‑93 (MQ=255)
                                                          cAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTtcatc                                    >  1:3720937/1‑93 (MQ=255)
                                                          cAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTtcatc                                    >  1:6367038/1‑93 (MQ=255)
                                                          cAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTtcatc                                    <  1:93500/93‑1 (MQ=255)
                                                           aGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTtcatca                                   >  1:1644470/1‑93 (MQ=255)
                                                             ttGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAAc                                 >  1:2824237/1‑93 (MQ=255)
                                                             ttGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAAc                                 >  1:1037415/1‑93 (MQ=255)
                                                             ttGCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAAc                                 >  1:2056389/1‑93 (MQ=255)
                                                               gCTGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGa                               <  2:4483423/93‑1 (MQ=255)
                                                               gCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGa                               >  1:781040/1‑93 (MQ=255)
                                                                 tGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAgc                             >  2:510797/1‑93 (MQ=255)
                                                                 tGGCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAgc                             >  1:3360934/1‑93 (MQ=255)
                                                                  ggCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAgcg                            >  1:3963238/1‑93 (MQ=255)
                                                                   gCTCTCTTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAgcgc                           >  2:4106944/1‑93 (MQ=255)
                                                                    ctctctTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGcc                          >  1:7429093/1‑93 (MQ=255)
                                                                     tctctTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCTCCt                         <  1:6980094/93‑1 (MQ=255)
                                                                      ctctTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCtt                        <  2:5915867/93‑1 (MQ=255)
                                                                      ctctTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCtt                        >  1:2671914/1‑93 (MQ=255)
                                                                      ctctTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCtt                        >  2:6814171/1‑93 (MQ=255)
                                                                       tctTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTa                       <  1:2574085/93‑1 (MQ=255)
                                                                        ctTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTaa                      <  1:487972/93‑1 (MQ=255)
                                                                        ctTTTTCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTaa                      <  2:3769520/93‑1 (MQ=255)
                                                                         tttttCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAg                     >  1:6253588/1‑93 (MQ=255)
                                                                         tttttCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAg                     <  1:5556686/93‑1 (MQ=255)
                                                                            ttCCGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACg                  <  2:3161460/93‑1 (MQ=255)
                                                                              ccGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTg                <  2:6020629/93‑1 (MQ=255)
                                                                              ccGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAATTCATCAACGAGCGCCTTAAGACGTg                >  1:7530443/1‑93 (MQ=255)
                                                                               cGCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGc               >  2:265067/1‑93 (MQ=255)
                                                                                gCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCa              >  1:3195934/1‑93 (MQ=255)
                                                                                gCCAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCa              >  1:5443729/1‑93 (MQ=255)
                                                                                  cAGCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCaga            <  1:2390856/93‑1 (MQ=255)
                                                                                  cAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCCACGAGCGCCTTAAGACGTGCaga            >  1:6682328/1‑93 (MQ=255)
                                                                                    gCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGt          >  2:6379744/1‑93 (MQ=255)
                                                                                    gCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGt          <  2:7120743/93‑1 (MQ=255)
                                                                                    gCTCTGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGt          >  2:3592850/1‑93 (MQ=255)
                                                                                     ctctGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTc         >  1:1909730/1‑93 (MQ=255)
                                                                                     ctctGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTc         <  2:570766/93‑1 (MQ=255)
                                                                                      tctGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCa        <  2:1217100/93‑1 (MQ=255)
                                                                                      tctGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCa        >  2:2488570/1‑93 (MQ=255)
                                                                                      tctGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCa        >  1:5734257/1‑93 (MQ=255)
                                                                                        tGAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGa      >  2:4781107/1‑93 (MQ=255)
                                                                                         gAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGAc     <  1:5147658/93‑1 (MQ=255)
                                                                                         gAGTGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGAc     <  1:651507/93‑1 (MQ=255)
                                                                                         gAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGAc     >  1:5506535/1‑93 (MQ=255)
                                                                                            cGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGc  <  2:6352812/93‑1 (MQ=255)
                                                                                            |                                                                                            
CCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGTTATTGGCAGACAGCTGCGCCAGCTCCAGTTCATCAACGAGCGCCTTAAGACGTGCAGAGTCAGACTGC  >  NC_000913/4386237‑4386421

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