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,8650TC99.8% 1382.5 / ‑6.3 425K157R (AAG→AGG) dppAdipeptide transport protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (228/196);  total (229/196)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 450 total aligned reads displayed.

GACAGAATAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCTTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGTTTTGCGCG  >  NC_000913/3704775‑3704954
                                                                                          |                                                                                         
taCAGAATAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCctt                                                                                         <  1:2268550/92‑1 (MQ=255)
 aCAGAATAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCcttt                                                                                        <  2:8553025/93‑1 (MQ=255)
  cAGAATAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCctttt                                                                                       >  1:8594291/1‑93 (MQ=255)
    gAATAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTttca                                                                                     <  1:2876727/93‑1 (MQ=255)
      aTAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTttcact                                                                                   >  2:9170113/1‑93 (MQ=255)
      aTAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTAGTCCACCCTTttcact                                                                                   <  2:7135177/93‑1 (MQ=255)
       tAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTttcgctt                                                                                  >  1:1587624/1‑93 (MQ=255)
       tAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTttcactt                                                                                  >  2:1071767/1‑93 (MQ=255)
         gagGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTttcacttca                                                                                <  1:7177565/93‑1 (MQ=255)
         gagGCGAAGTCCATTGACAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTttcacttca                                                                                <  1:8985238/93‑1 (MQ=255)
            gCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTg                                                                             <  1:9553873/93‑1 (MQ=255)
             ggAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGa                                                                            <  1:6587019/92‑1 (MQ=255)
             cGAAGTCCATTGCCATGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGa                                                                            <  2:508296/93‑1 (MQ=255)
             cGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGa                                                                            <  2:635431/93‑1 (MQ=255)
             cGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGa                                                                            >  2:6324679/1‑93 (MQ=255)
                aGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGCTCa                                                                         >  1:2633374/1‑93 (MQ=255)
                aGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCAACCCTTTTCACTTCACTGATCa                                                                         <  1:1811089/93‑1 (MQ=255)
                 gTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAg                                                                        >  1:7606691/1‑93 (MQ=255)
                 gTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAg                                                                        <  1:6010503/93‑1 (MQ=255)
                  tCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGc                                                                       <  1:9804137/93‑1 (MQ=255)
                  tCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGc                                                                       >  1:2204279/1‑93 (MQ=255)
                  tCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGc                                                                       >  1:4636310/1‑93 (MQ=255)
                  tCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGc                                                                       <  2:4900375/93‑1 (MQ=255)
                    cATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGctc                                                                     >  1:4756161/1‑93 (MQ=255)
                    cATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGctc                                                                     >  2:9369430/1‑93 (MQ=255)
                      ttGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTg                                                                   >  2:7060637/1‑93 (MQ=255)
                      ttGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTg                                                                   >  1:5582351/1‑93 (MQ=255)
                      ttGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTg                                                                   >  2:6701057/1‑93 (MQ=255)
                      ttGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTg                                                                   <  1:6713398/93‑1 (MQ=255)
                           aGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAg                                                              >  1:4037200/1‑93 (MQ=255)
                            ggTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGc                                                             >  2:4504338/1‑93 (MQ=255)
                             gTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGcc                                                            >  1:3625974/1‑93 (MQ=255)
                               cAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCa                                                          <  1:6185004/93‑1 (MQ=255)
                               cAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCa                                                          <  2:422814/93‑1 (MQ=255)
                               cAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCa                                                          <  2:5826943/93‑1 (MQ=255)
                               cAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCa                                                          <  2:8262128/93‑1 (MQ=255)
                                aGCGAGGAACGGCGCTTCCGTGCGAGTCAGCACAAACTGAACGGTGGTGTCGTCCACCCTTTTCACTTAACTGATCAGCTCTGGCAAGACCac                                                         >  2:7589644/1‑92 (MQ=255)
                                aGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCAt                                                         >  1:8740191/1‑93 (MQ=255)
                                  cGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGc                                                       <  2:2843087/93‑1 (MQ=255)
                                    aGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCt                                                     >  1:203068/1‑93 (MQ=255)
                                     ggAACCGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCACCACTTTACACTTCACTGATCAGCTCTGGCAAGCCCATGCCtt                                                    <  2:7719796/93‑1 (MQ=255)
                                      gAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTc                                                   >  2:3753213/1‑93 (MQ=255)
                                       aaCGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGGTCAGCTCTGGCAAGCCCATGCCTTCg                                                  <  2:6508135/93‑1 (MQ=255)
                                        aCGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGa                                                 <  1:4187714/93‑1 (MQ=255)
                                        aCGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGa                                                 >  2:946824/1‑93 (MQ=255)
                                         cGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGaa                                                <  2:5955546/93‑1 (MQ=255)
                                           gcgcTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTTATCAGCTCTGGCAAGCCCATGCCTTCGAAGt                                              >  1:6798123/1‑93 (MQ=255)
                                            cgcTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTa                                             <  2:5311510/93‑1 (MQ=255)
                                            cgcTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTa                                             <  2:7683102/93‑1 (MQ=255)
                                            cgcTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGGCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTa                                             <  1:7648724/93‑1 (MQ=255)
                                             gcTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCGAGCCCATGCGTTCGAAGTAt                                            >  2:4937537/1‑93 (MQ=255)
                                             gcTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTAt                                            <  2:7487106/93‑1 (MQ=255)
                                              cTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTAtt                                           <  1:6412864/93‑1 (MQ=255)
                                              cTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTAtt                                           <  2:7276681/93‑1 (MQ=255)
                                              cTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTAtt                                           <  2:790693/93‑1 (MQ=255)
                                              cTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTAtt                                           <  1:8019742/93‑1 (MQ=255)
                                               ttCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTc                                          >  1:1702969/1‑93 (MQ=255)
                                                tCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCg                                         >  1:3521551/1‑93 (MQ=255)
                                                    ggCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGc                                     <  2:2669848/93‑1 (MQ=255)
                                                    ggCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGc                                     >  2:6839170/1‑93 (MQ=255)
                                                    ggCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGc                                     <  1:8570619/93‑1 (MQ=255)
                                                     gCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCt                                    >  2:5179682/1‑93 (MQ=255)
                                                     gCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCt                                    >  1:2986290/1‑93 (MQ=255)
                                                        aGTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTgcc                                 >  1:8096854/1‑93 (MQ=255)
                                                         gTCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTgccg                                <  1:1333969/93‑1 (MQ=255)
                                                          tCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGTAGTATTCGTAGCTgccgc                               >  1:9403664/1‑93 (MQ=255)
                                                          tCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTgccgc                               >  2:3837272/1‑93 (MQ=255)
                                                          tCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTgccgc                               >  2:1023448/1‑93 (MQ=255)
                                                          tCAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTgccgc                               >  1:3329237/1‑93 (MQ=255)
                                                           cAGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTgccgcc                              <  2:1660961/93‑1 (MQ=255)
                                                            aGCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCa                             >  2:1352966/1‑93 (MQ=255)
                                                             gCACAAACTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAg                            >  1:6072726/1‑93 (MQ=255)
                                                                  aaCTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACt                       <  2:2088561/93‑1 (MQ=255)
                                                                   aCTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACtt                      <  1:3752629/93‑1 (MQ=255)
                                                                    cTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACttt                     <  1:4973088/93‑1 (MQ=255)
                                                                    cTGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACttt                     <  1:10031559/93‑1 (MQ=255)
                                                                     tGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTTGGAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTa                    >  2:6063360/1‑93 (MQ=255)
                                                                     tGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTa                    >  2:4560287/1‑93 (MQ=255)
                                                                     tGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTa                    <  2:9722536/93‑1 (MQ=255)
                                                                     tGAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTa                    <  1:5159726/93‑1 (MQ=255)
                                                                      gAACGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTAt                   >  2:4306781/1‑93 (MQ=255)
                                                                       aaCGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATg                  <  1:5298461/93‑1 (MQ=255)
                                                                        aCGGTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATgg                 <  2:8738911/93‑1 (MQ=255)
                                                                         cGGGGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGt                >  2:9615989/1‑93 (MQ=255)
                                                                          ggTGTTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTa               <  2:3946402/93‑1 (MQ=255)
                                                                          cgtgtTGTCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTa               <  1:1538157/92‑1 (MQ=255)
                                                                             gttgtCGTCCACCCTTTTCACTTCACTGATCAGCTCAGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACgg            >  1:569602/1‑93 (MQ=255)
                                                                             gttgaCGTCCCTCCTTTTCCCTTCGCTGATCAGCTCTGGCAAGCTCATGCCTTCGAAGTATTAGTAGCTGCCGCCAGAAACTTTATGGTACgg            <  2:1849635/93‑1 (MQ=255)
                                                                              ttgtCGTCCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACggg           >  1:4476163/1‑93 (MQ=255)
                                                                                gtcgtcCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGtt         <  2:3094707/93‑1 (MQ=255)
                                                                                gtcgtcCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGtt         >  1:2118542/1‑93 (MQ=255)
                                                                                gtcgtcCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGtt         >  1:6928230/1‑93 (MQ=255)
                                                                                gtcgtcCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGtt         >  1:9200963/1‑93 (MQ=255)
                                                                                gtcgtcCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGtt         >  2:3298135/1‑93 (MQ=255)
                                                                                 tcgtcCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGttt        <  2:1487310/93‑1 (MQ=255)
                                                                                  cgtcCCCCCTTTTCACTTCACTGTGCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGtttt       <  2:5502615/93‑1 (MQ=255)
                                                                                    tcCACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGTTTTgc     >  2:5005501/1‑93 (MQ=255)
                                                                                      cACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGTTTTgcgc   >  2:8624091/1‑93 (MQ=255)
                                                                                      cACCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGTTTTgcgc   >  2:2227241/1‑93 (MQ=255)
                                                                                       aCCCTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGTTTTGcgcg  >  2:3433362/1‑93 (MQ=255)
                                                                                          |                                                                                         
GACAGAATAGAGGCGAAGTCCATTGCCAGGTCAGCGAGGAACGGCGCTTCCGGGCGAGTCAGCACAAACTGAACGGTGTTGTCGTCCACCTTTTTCACTTCACTGATCAGCTCTGGCAAGCCCATGCCTTCGAAGTATTCGTAGCTGCCGCCAGAAACTTTATGGTACGGGTTTTGCGCG  >  NC_000913/3704775‑3704954

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