Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,704,865 T→C K157R (AAG→AGG)  dppA ← dipeptide transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,704,8650TC98.0% 466.9 / ‑1.0 151K157R (AAG→AGG) dppAdipeptide transport protein
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (81/66);  total (84/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.44e-01

Only 100 of 164 total aligned reads displayed.

GACAGAATAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCTTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGTTTTGCGCGTT  >  NC_000913/3704775‑3704956
                                                                                          |                                                                                           
gACAGAATAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCctt                                                                                           >  1:3357182/1‑93 (MQ=255)
gACAGAATAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCctt                                                                                           >  1:4052138/1‑93 (MQ=255)
 aCAGAATAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCcttt                                                                                          >  1:1216436/1‑93 (MQ=255)
   aGAATAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGTGGGAGGCAGCCAAAACTTAACGGTGCTGGCGTCCACCCTTttc                                                                                        >  1:2313452/1‑93 (MQ=255)
    gAATAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTttca                                                                                       >  1:3348214/1‑93 (MQ=255)
     aaTAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTttcac                                                                                      <  2:2492027/93‑1 (MQ=255)
     aaTAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTttcac                                                                                      <  2:2492026/93‑1 (MQ=255)
       tAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTttcactt                                                                                    >  1:4162270/1‑93 (MQ=255)
        agagGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCCCCCTTttcacttc                                                                                   >  1:2764934/1‑93 (MQ=255)
        agagGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTttcacttc                                                                                   >  1:1545019/1‑93 (MQ=255)
         gagGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTttcacttca                                                                                  >  1:1721373/1‑93 (MQ=255)
           ggCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTacacttcact                                                                                >  2:4535657/1‑93 (MQ=255)
            gCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTg                                                                               <  1:4102673/93‑1 (MQ=255)
            gCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTg                                                                               >  2:4456170/1‑93 (MQ=255)
            gCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTg                                                                               <  2:3868319/93‑1 (MQ=255)
              gAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGAt                                                                             >  2:1434127/1‑93 (MQ=255)
               aaGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATc                                                                            >  2:510179/1‑93 (MQ=255)
               aaGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATc                                                                            <  2:2683054/93‑1 (MQ=255)
               aaGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATc                                                                            >  1:195913/1‑93 (MQ=255)
               aaGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATc                                                                            >  2:1862921/1‑93 (MQ=255)
               aaGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATc                                                                            >  2:1427827/1‑93 (MQ=255)
               aaGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATc                                                                            >  2:1783848/1‑93 (MQ=255)
                aGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCa                                                                           >  2:2177975/1‑93 (MQ=255)
                aGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCa                                                                           <  2:3434748/93‑1 (MQ=255)
                 gTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCTTTTTCACTTCACTGATCAg                                                                          >  2:1580110/1‑93 (MQ=255)
                 gTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAg                                                                          >  1:3675830/1‑93 (MQ=255)
                  tCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGc                                                                         <  2:3453388/93‑1 (MQ=255)
                   ccATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGct                                                                        >  2:786714/1‑93 (MQ=255)
                    cATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGctc                                                                       <  1:3814941/93‑1 (MQ=255)
                    cATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGctc                                                                       <  2:4300894/93‑1 (MQ=255)
                     aTTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGctct                                                                      >  1:2971287/1‑93 (MQ=255)
                      ttGCCAGGTCAGCGAGGCACGGCGCTACCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTg                                                                     <  1:4216638/93‑1 (MQ=255)
                      ttGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTg                                                                     >  1:4663672/1‑93 (MQ=255)
                      ttGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTg                                                                     <  1:107148/93‑1 (MQ=255)
                      ttGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTg                                                                     <  2:2929583/93‑1 (MQ=255)
                       tGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTgg                                                                    <  2:4613701/93‑1 (MQ=255)
                       tGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCATTTTCACTTCACTGATCAGCTCTgg                                                                    >  1:4490292/1‑93 (MQ=255)
                          cAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCaa                                                                 <  2:1725202/93‑1 (MQ=255)
                          cAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCaa                                                                 <  1:1429479/93‑1 (MQ=255)
                           aGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAg                                                                >  1:3128691/1‑93 (MQ=255)
                            ggTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGc                                                               >  2:3949163/1‑93 (MQ=255)
                             gTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGcc                                                              <  2:3104524/93‑1 (MQ=255)
                             gTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTGTTCACTTCACTGATCAGCTCGGGCAAGcc                                                              >  2:4209748/1‑93 (MQ=255)
                               cAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCa                                                            <  1:1929822/93‑1 (MQ=255)
                               cAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCa                                                            <  1:1630073/93‑1 (MQ=255)
                                 gCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATg                                                          <  1:321028/93‑1 (MQ=255)
                                   gAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGACCATGcc                                                        >  2:590934/1‑93 (MQ=255)
                                    aGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCt                                                       >  1:178225/1‑93 (MQ=255)
                                    aGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCt                                                       >  2:1905742/1‑93 (MQ=255)
                                    aGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCt                                                       >  1:89985/1‑93 (MQ=255)
                                     ggAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCtt                                                      <  1:4184693/93‑1 (MQ=255)
                                      gAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTc                                                     >  2:150926/1‑93 (MQ=255)
                                      gAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTc                                                     >  1:4588003/1‑93 (MQ=255)
                                        aCGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGAAAGCCCATGCCTTCGa                                                   >  1:1269449/1‑93 (MQ=255)
                                         cGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGaa                                                  <  2:2205504/93‑1 (MQ=255)
                                          ggCGCTTCCGGGCTAGTCAGCACAAACTGGACGGTGTTGTCGTCCACCCTTTTCACTTCACTGGTCGGCTTTGGCCAGGACATGGCTTCGacg                                                 >  2:1633299/1‑91 (MQ=255)
                                             gcTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTAt                                              <  2:3285095/93‑1 (MQ=255)
                                             gcTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTAt                                              >  2:3081939/1‑93 (MQ=255)
                                              cTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTAtt                                             <  1:1192762/93‑1 (MQ=255)
                                              cTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTAtt                                             <  1:2697907/93‑1 (MQ=255)
                                               ttCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGa                                                   <  1:3793730/86‑1 (MQ=255)
                                               ttCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTc                                            <  2:3982759/93‑1 (MQ=255)
                                               ttCCGGGCGAGTCAGCAAAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTc                                            <  2:2011522/93‑1 (MQ=255)
                                                tCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCTCTGCAATGATCAGCGCTGGAAAGCCCATGCCCTCGCattattat                                           >  2:3788003/1‑91 (MQ=255)
                                                  cGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTa                                         >  2:1814792/1‑93 (MQ=255)
                                                  cGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTa                                         <  1:2776551/93‑1 (MQ=255)
                                                  cGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTa                                         <  2:2417362/93‑1 (MQ=255)
                                                   gggCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAg                                        <  2:3647756/93‑1 (MQ=255)
                                                    ggCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGGAgt                                       >  2:2655414/1‑92 (MQ=255)
                                                     gCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGTGAAGCCCATGCCTTCGAAGGATTCGTAGCt                                      >  2:2555129/1‑93 (MQ=255)
                                                       gAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTgc                                    >  2:272345/1‑93 (MQ=255)
                                                       gAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTgc                                    >  1:507194/1‑93 (MQ=255)
                                                       gAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTgc                                    <  1:4638112/93‑1 (MQ=255)
                                                       gAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTgc                                    >  2:1371893/1‑93 (MQ=255)
                                                         gTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTgccg                                  >  2:3091588/1‑93 (MQ=255)
                                                         gTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTgccg                                  >  1:4521535/1‑93 (MQ=255)
                                                         gTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTgccg                                  >  2:1234603/1‑93 (MQ=255)
                                                            aGCACAAACTGAACGGTGCTGTCGTCCACCCTTTTCACTTCGCTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCa                               <  2:1300980/93‑1 (MQ=255)
                                                              cacaAACGGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATGCGTAGCTGCCGCCAGa                             <  2:611784/93‑1 (MQ=255)
                                                               acaAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAACTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGaa                            >  2:4133302/1‑93 (MQ=255)
                                                                caAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGaaa                           <  1:581263/93‑1 (MQ=255)
                                                                 aaaCTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAAc                          <  1:3233095/93‑1 (MQ=255)
                                                                   aCTGAATGGTGTTGTCGTCCCCCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACtt                        <  1:3450718/93‑1 (MQ=255)
                                                                    cTGAACGGTGTTGTCGTCCACCCTTTTCCCTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACttt                       <  1:1334599/93‑1 (MQ=255)
                                                                     tGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTa                      <  1:2517484/93‑1 (MQ=255)
                                                                      gAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTAt                     <  1:3419685/93‑1 (MQ=255)
                                                                       aaCGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATg                    <  2:2289878/93‑1 (MQ=255)
                                                                            tgttgtCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACg               <  1:4137456/93‑1 (MQ=255)
                                                                             gttgtCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACgg              >  2:4134241/1‑93 (MQ=255)
                                                                              ttgtCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACggg             >  1:1869855/1‑93 (MQ=255)
                                                                              ttgtCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCCCCAGAAACTTTATGGTACggg             >  2:1199554/1‑93 (MQ=255)
                                                                               tgtCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGt            <  2:3769499/93‑1 (MQ=255)
                                                                                gtcgtcCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGtt           >  1:2260752/1‑93 (MQ=255)
                                                                                 tcgtcCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGttt          <  1:57194/93‑1 (MQ=255)
                                                                                 tcgtcCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGttt          <  2:4354919/93‑1 (MQ=255)
                                                                                    tcCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGTTTTgc       >  2:1297596/1‑93 (MQ=255)
                                                                                    tcCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGTTTTgc       >  2:1904632/1‑93 (MQ=255)
                                                                                       aCCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGTTTTGcgcg    >  1:3271540/1‑93 (MQ=255)
                                                                                       aCCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGTTTTGcgcg    >  1:3008432/1‑93 (MQ=255)
                                                                                       aCCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGTTTTGcgcg    >  2:2025233/1‑93 (MQ=255)
                                                                                          |                                                                                           
GACAGAATAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCTTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGTTTTGCGCGTT  >  NC_000913/3704775‑3704956

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