Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,840,547 T→C intergenic (+39/+2) yicS → / ← nepI putative periplasmic protein/putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,840,5470TC97.0% 453.7 / ‑1.3 134intergenic (+39/+2)yicS/nepIputative periplasmic protein/putative transporter
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (64/66);  total (66/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

TAATAATCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAATA  >  NC_000913/3840451‑3840642
                                                                                                |                                                                                               
taataaTCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATc                                                                                              >  1:1848728/1‑100 (MQ=255)
        agagCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttc                                                                                      >  2:2808454/1‑100 (MQ=255)
          agCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttctt                                                                                    >  2:1120732/1‑100 (MQ=255)
          agCAATTATTGGGAAGCATTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttctt                                                                                    <  1:773935/100‑1 (MQ=255)
           gCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttcttc                                                                                   <  1:3002512/100‑1 (MQ=255)
           gCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttcttc                                                                                   >  1:1175570/1‑100 (MQ=255)
            cAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCa                                                                                  <  1:2465861/100‑1 (MQ=255)
            cAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCa                                                                                  <  2:693374/100‑1 (MQ=255)
             aaTTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCAt                                                                                 <  1:1139104/100‑1 (MQ=255)
              attattGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCAtt                                                                                >  1:2251751/1‑100 (MQ=255)
              attattGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCAtt                                                                                >  2:2982561/1‑100 (MQ=255)
               ttataGGCAAGCGATGGGTAAGGTAGAATGCCTGGATATGTAAGAACAGCAACAACTACTGACGAACGTGTCTTTACCTGACATCAGCATTTCTTCAtt                                                                                >  1:2558900/1‑99 (MQ=255)
                tattGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCAtttt                                                                              <  1:2813490/100‑1 (MQ=255)
                 attGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTc                                                                             <  1:1268838/100‑1 (MQ=255)
                  ttGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCa                                                                            >  2:192038/1‑100 (MQ=255)
                  ttGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCa                                                                            <  1:2191417/100‑1 (MQ=255)
                   ttggAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCAc                                                                           <  2:2025581/98‑1 (MQ=255)
                   tGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCAc                                                                           <  2:169975/100‑1 (MQ=255)
                   tGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCAc                                                                           <  2:1179232/100‑1 (MQ=255)
                    gggCAGCGTTGGGTAAGGTAGAATGTCAGGATATGTAAGAACAGAAACCACTACTGACCACCGCTTCTTTACATGATATCCGGATTTATTCATTTTCa                                                                            >  2:695822/1‑98 (MQ=255)
                    gggAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCAcc                                                                          >  1:3080810/1‑100 (MQ=255)
                    gggAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCAcc                                                                          <  2:1555003/100‑1 (MQ=255)
                    gggAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCAcc                                                                          <  2:662324/100‑1 (MQ=255)
                    gggAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCAcc                                                                          >  2:2991301/1‑100 (MQ=255)
                    gggAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCAca                                                                          >  1:1179217/1‑99 (MQ=255)
                         gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGATCAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGc                                                                     <  2:1945779/100‑1 (MQ=255)
                         gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTTCTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGc                                                                     <  1:1338549/100‑1 (MQ=255)
                         gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGc                                                                     <  2:1583095/100‑1 (MQ=255)
                         gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGc                                                                     <  1:1140210/100‑1 (MQ=255)
                         gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGc                                                                     <  1:1642364/100‑1 (MQ=255)
                          cGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCa                                                                    <  2:1544948/100‑1 (MQ=255)
                           gTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAg                                                                   <  1:2177996/100‑1 (MQ=255)
                           gTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAg                                                                   >  2:21830/1‑100 (MQ=255)
                           gTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAg                                                                   <  1:437746/100‑1 (MQ=255)
                              gggTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTaa                                                                >  2:2374088/1‑100 (MQ=255)
                              gggTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTaa                                                                >  2:2373571/1‑100 (MQ=255)
                                gTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAAcc                                                              <  1:3226953/100‑1 (MQ=255)
                                 tAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACca                                                             >  1:3450474/1‑100 (MQ=255)
                                   aGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaac                                                           <  2:3146801/100‑1 (MQ=255)
                                   aGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaac                                                           >  2:136449/1‑100 (MQ=255)
                                   aGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaac                                                           >  1:5496/1‑100 (MQ=255)
                                   aGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaac                                                           <  1:669219/100‑1 (MQ=255)
                                    ggTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaaca                                                          >  2:2990734/1‑100 (MQ=255)
                                     gTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaacaa                                                         >  1:1840038/1‑100 (MQ=255)
                                     gTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaacaa                                                         >  2:2075594/1‑100 (MQ=255)
                                        gAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGc                                                      >  2:2173389/1‑100 (MQ=255)
                                        gAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGc                                                      <  2:1688259/100‑1 (MQ=255)
                                            gCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTc                                                  <  1:159833/100‑1 (MQ=255)
                                             cccGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCa                                                 <  2:3353594/100‑1 (MQ=255)
                                              ccgatcttGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAg                                                <  1:1571659/93‑1 (MQ=255)
                                              ccGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAg                                                <  2:2194847/100‑1 (MQ=255)
                                                  atatGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAAc                                            <  1:3045302/100‑1 (MQ=255)
                                                  atatGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAAc                                            <  2:3056578/100‑1 (MQ=255)
                                                   tatGTAAGAACAGCCACAACTACTGACCAACGCTTCTTTTCCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAAc                                            <  1:203162/99‑1 (MQ=255)
                                                     tGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATc                                         >  2:2653526/1‑100 (MQ=255)
                                                     tGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATc                                         >  2:2713215/1‑100 (MQ=255)
                                                     tGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATc                                         >  1:1638816/1‑100 (MQ=255)
                                                      gTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCa                                        <  2:2465862/100‑1 (MQ=255)
                                                      gTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTAACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCa                                        >  2:1154700/1‑100 (MQ=255)
                                                       tAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCaa                                       >  1:2044259/1‑100 (MQ=255)
                                                        aaGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAAt                                      <  2:1314342/100‑1 (MQ=255)
                                                         aGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCaa                                       >  1:1032413/1‑98 (MQ=255)
                                                           aaCAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAAtgtn                                   <  2:1428846/100‑2 (MQ=255)
                                                             cAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGcc                                 >  1:705072/1‑100 (MQ=255)
                                                             cAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTAGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGc                                  >  1:1437412/1‑99 (MQ=255)
                                                              aGCAACAATTACTGACCAACGCTTCTTTACCTGACAGCAGGATTTCTTCATTTTCACCTTTGCAGTAACAAACAATGGAGTCAGCAAGATCatggggcc                                 <  1:152393/99‑8 (MQ=255)
                                                               gCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCgg                               >  2:2110487/1‑100 (MQ=255)
                                                               gCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCgg                               >  2:535594/1‑100 (MQ=255)
                                                               gCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCgg                               >  1:551978/1‑100 (MQ=255)
                                                                 aacaacTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGac                             <  1:210024/100‑1 (MQ=255)
                                                                   caacTACTGACCAACGCTTCTTTTCCTGACGTCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGaca                            <  1:460510/99‑1 (MQ=255)
                                                                   caacTACTGACCAACGCTTCTTTCCCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGacaa                           >  2:1199729/1‑100 (MQ=255)
                                                                     actactGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGacaaca                         <  1:1106769/100‑1 (MQ=255)
                                                                      ctactGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAGCAt                        >  2:2070338/1‑100 (MQ=255)
                                                                      ctactGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACAt                        >  2:2605120/1‑100 (MQ=255)
                                                                      ctactGACCAACGCTTCTTTACCTGACACCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGacaaca                         <  1:314006/99‑1 (MQ=255)
                                                                       tactGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATc                       <  2:685222/100‑1 (MQ=255)
                                                                         ntGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCaa                     >  2:129539/2‑100 (MQ=255)
                                                                         cttaCCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCa                      <  1:2602455/96‑1 (MQ=255)
                                                                         ctGCCCAACGCTGCTTTACCGGACATCAGGATTTCTTCATTTGCACCTTTGCAGGAAGAAACAATGCAGGCAGCAACATCAATGTGCGGGacaac                          <  1:3143333/95‑1 (MQ=255)
                                                                         ctGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCaa                     >  2:402654/1‑100 (MQ=255)
                                                                         ctGACCAACGCTTCTTTACCTGACAGCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATc                       <  1:1707521/98‑1 (MQ=255)
                                                                           gACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACg                   >  2:2804925/1‑100 (MQ=255)
                                                                            aCCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACgg                  >  2:624086/1‑100 (MQ=255)
                                                                             ccAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGc                 <  2:2633920/100‑1 (MQ=255)
                                                                              cAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCg                >  1:912516/1‑100 (MQ=255)
                                                                              cAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCg                <  2:330243/100‑1 (MQ=255)
                                                                                 cGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAg             <  1:70688/100‑1 (MQ=255)
                                                                                   cttcttTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTc           <  2:531220/100‑1 (MQ=255)
                                                                                     tcttTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAg         >  2:1219257/1‑100 (MQ=255)
                                                                                     tcttTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAg         >  2:3239004/1‑100 (MQ=255)
                                                                                     tcttTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAg         <  1:2796231/100‑1 (MQ=255)
                                                                                      cttTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAgn        <  2:3136580/100‑2 (MQ=255)
                                                                                      cttTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGa        <  2:803542/100‑1 (MQ=255)
                                                                                      ccttACCTGATATCATTATTTCTGCAGTTGCACCTTTTCATTAACCAACAATTCAGGCAGCAACATCAATGTTCCGGACAACATCAACGGCGAAGTCAGa        <  1:695813/98‑1 (MQ=255)
                                                                                        ttACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAcc      >  2:828880/1‑100 (MQ=255)
                                                                                        ttACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAcc      <  1:953262/100‑1 (MQ=255)
                                                                                            cTGACATCAGGATTTCTTCATTTTTATTTTTGTTGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAata  >  1:1071963/1‑100 (MQ=255)
                                                                                            cTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAata  >  2:582013/1‑100 (MQ=255)
                                                                                            cTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAata  >  1:1980689/1‑100 (MQ=255)
                                                                                                |                                                                                               
TAATAATCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAATA  >  NC_000913/3840451‑3840642

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