Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,840,547 T→C 100% 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,5470TC98.7% 537.8 / ‑3.2 156intergenic (+39/+2)yicS/nepIputative periplasmic protein/putative transporter
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (79/75);  total (79/77)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.84e-01

Only 100 of 157 total aligned reads displayed.

TAATAATCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAATAT  >  NC_000913/3840451‑3840643
                                                                                                |                                                                                                
taataaTCTGAGCAATTATTGGGAAGCGTTGGGTAATATAGAATGCCCGGATATGTAAGACCAGCAACAACTACTGACCAACGCTTCTTTACTTGACATc                                                                                               <  1:1590388/100‑1 (MQ=255)
taataaTCAGAGCAATTATTGGGAATCGTTTGGTAAGGTAGAATGTCCGGTTATGTAAGAACAGCAACAACTGCTGACCAACGCTTCTTTACCTGACATc                                                                                               <  2:552135/100‑1 (MQ=255)
  ataatCAGAGCAATTATTGTGAAGCGTTGGGTATTGTAGTATTTCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAg                                                                                             <  1:202265/100‑1 (MQ=255)
  ataatCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAg                                                                                             <  1:2802467/100‑1 (MQ=255)
    aatCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATGTGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAgg                                                                                            >  2:2273605/1‑99 (MQ=255)
    aatCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGa                                                                                           <  2:2876146/100‑1 (MQ=255)
      tCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGAtt                                                                                         <  1:2365663/100‑1 (MQ=255)
          agCAATTATTGGGAAGCGTTGGGTAGGGTAGAATGCCCGGATATGTAAGAACATCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttctt                                                                                     <  2:2470444/100‑1 (MQ=255)
          agCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGCCCAACGCTTCTTTACCTGACATCAGGATttctt                                                                                     <  2:2007223/100‑1 (MQ=255)
           gCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttcttc                                                                                    >  2:53855/1‑100 (MQ=255)
            cAATTATTGGGAAGCGTTGGGTAAGGTAGAATGACCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCa                                                                                   >  2:1789857/1‑100 (MQ=255)
             aaTTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCAt                                                                                  <  1:1528271/100‑1 (MQ=255)
                 attGGGAAGCGTTGGGTAAGGTAGAATGCCAGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTc                                                                              <  1:442963/100‑1 (MQ=255)
                   tGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCAc                                                                            <  2:2195383/100‑1 (MQ=255)
                    gggAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCAcc                                                                           >  1:3034902/1‑100 (MQ=255)
                      gAAGCGTTGGGTAAGGTAGAATGCCCGGGTATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATAAGAATTTCTTCATTTTCACCtt                                                                         >  1:745314/1‑100 (MQ=255)
                         gCGTTGTGTAAGGTAGAATGCCTGGATATGTAAGAACAGCAACCTCTACTGACCAACGCTTTTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGc                                                                      <  2:1344589/100‑1 (MQ=255)
                         gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGc                                                                      <  1:3027994/100‑1 (MQ=255)
                         gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGc                                                                      <  2:2820027/100‑1 (MQ=255)
                         gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGc                                                                      <  2:2317443/100‑1 (MQ=255)
                         gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGc                                                                      <  1:193961/100‑1 (MQ=255)
                          cGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCa                                                                     <  1:527222/100‑1 (MQ=255)
                          cGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCa                                                                     >  1:2552583/1‑100 (MQ=255)
                           gTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAg                                                                    <  1:293955/100‑1 (MQ=255)
                           gTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAAAAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTATTCATTTTCACCTTTGCAg                                                                    >  2:740068/1‑100 (MQ=255)
                             tGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTa                                                                  >  1:1557918/1‑100 (MQ=255)
                              gggTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTaa                                                                 <  2:931171/100‑1 (MQ=255)
                              gggTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTaa                                                                 >  2:211057/1‑100 (MQ=255)
                               ggTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAAc                                                                >  1:1750761/1‑100 (MQ=255)
                               ggTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAAc                                                                <  2:788946/100‑1 (MQ=255)
                                gTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAAcc                                                               <  2:967408/100‑1 (MQ=255)
                                gTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAAcc                                                               >  2:79927/1‑100 (MQ=255)
                                 tAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACca                                                              <  2:2966539/100‑1 (MQ=255)
                                 tAAGGTAGAATGCCCGAATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAAcaa                                                              <  2:1196366/100‑3 (MQ=255)
                                  aaGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaa                                                             >  1:2615390/1‑100 (MQ=255)
                                   aGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaac                                                            <  1:1177129/100‑1 (MQ=255)
                                   aGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaac                                                            >  1:1789309/1‑100 (MQ=255)
                                   aGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTAACTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaac                                                            >  2:427400/1‑100 (MQ=255)
                                    ggTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCCCCTTTGCAGTAACcaaca                                                           <  2:1892080/100‑1 (MQ=255)
                                    ggTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaaca                                                           >  2:1375041/1‑100 (MQ=255)
                                     gTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaacaa                                                          >  1:210726/1‑100 (MQ=255)
                                       aGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGGCATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATg                                                        >  2:1043999/1‑100 (MQ=255)
                                       aGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATg                                                        <  2:2997268/100‑1 (MQ=255)
                                        gAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGc                                                       <  1:1294314/100‑1 (MQ=255)
                                        gAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGc                                                       <  1:2957037/100‑1 (MQ=255)
                                          aTGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAg                                                     <  2:1014752/100‑1 (MQ=255)
                                           tGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGt                                                    >  2:1616073/1‑100 (MQ=255)
                                            gCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTc                                                   >  2:2869744/1‑100 (MQ=255)
                                              ccGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAg                                                 >  1:347920/1‑100 (MQ=255)
                                               cGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAg                                                 >  2:1668878/1‑99 (MQ=255)
                                                   tatGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGAATTTCTTCATTTTCAACTTTGCAGTAACCAACAATGCAGTCAGCAACa                                            >  1:115209/1‑100 (MQ=255)
                                                    atGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACa                                            >  1:1984324/1‑99 (MQ=255)
                                                     tGTAATAACAGCAACAACTATTAACCAACGATTATTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACAAACAATGCAGTAAGAAACAt                                           >  2:1093679/1‑99 (MQ=255)
                                                      gTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCa                                         >  2:2430904/1‑100 (MQ=255)
                                                       tAAGAACAGCAACAACTACTGACTAACGCTTCTTTACCTGACCTCAGCATTTCTTCATTTTCACCTTTCCACTACCCAACAATGCAGTCAGCAACATCaa                                        >  2:1024886/1‑100 (MQ=255)
                                                         aGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAAtg                                      <  1:89579/100‑1 (MQ=255)
                                                         aGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAAtg                                      >  2:2936941/1‑100 (MQ=255)
                                                          gAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTACAGTAACCAACAATGAAGTCAGCAACATCAAtgt                                     >  1:580869/1‑100 (MQ=255)
                                                            aCAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGc                                   <  1:1048332/100‑1 (MQ=255)
                                                             cAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCAACTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGcc                                  >  1:878317/1‑100 (MQ=255)
                                                              aGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCg                                 <  1:393332/100‑1 (MQ=255)
                                                              aGCAACAACTACTGACCAACGCTTCGTTACCGGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCg                                 <  2:232132/100‑1 (MQ=255)
                                                               gCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCgg                                <  2:59122/100‑1 (MQ=255)
                                                                 aacaacTACTTACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGac                              <  2:59240/100‑1 (MQ=255)
                                                                 aacaacTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAAAAATGCAGTCAGCAACATCAATGTGCCGGac                              >  1:995040/1‑100 (MQ=255)
                                                                     actactGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGacaaca                          <  1:406099/100‑1 (MQ=255)
                                                                     actactGACCAACGCTTCTTTACCTGACAACAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCATTCCGCAACATCAATGTGCCGGacaaca                          >  2:123882/1‑100 (MQ=255)
                                                                      ctactGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACAt                         <  2:1839872/100‑1 (MQ=255)
                                                                      ctactGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAATAACAAACAATACAGTCAGCAAAATCAATGTGCCGAac                              >  1:2530034/1‑95 (MQ=255)
                                                                        actGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCa                       >  1:2221250/1‑100 (MQ=255)
                                                                         ctGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCaa                      >  2:72359/1‑100 (MQ=255)
                                                                           gACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACg                    >  2:2837668/1‑100 (MQ=255)
                                                                           gACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACg                    <  2:1941138/100‑1 (MQ=255)
                                                                             ccAACGCTTCTTTACCTGACATCAGGCTTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGc                  >  1:11763/1‑100 (MQ=255)
                                                                             ccAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCATCAACATCAATGTGCCGGTCAACATCAACgtc                  <  1:2047925/100‑3 (MQ=255)
                                                                             ccAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGc                  >  2:2219489/1‑100 (MQ=255)
                                                                             ccAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGc                  >  2:2711005/1‑100 (MQ=255)
                                                                             ccAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGCACAACATCAACGGc                  <  2:2083458/100‑1 (MQ=255)
                                                                             ccAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATACAGTCAGCAACATCAATGTGCCGGACAACATCAACGGc                  >  2:611459/1‑100 (MQ=255)
                                                                              cAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCg                 >  2:2660025/1‑100 (MQ=255)
                                                                              cAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCg                 >  2:1327120/1‑100 (MQ=255)
                                                                               aaCGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGc                  <  1:780309/98‑1 (MQ=255)
                                                                                 cGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAg              <  1:1048683/100‑1 (MQ=255)
                                                                                 cGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAg              >  1:2175780/1‑100 (MQ=255)
                                                                                 cGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAg              <  1:2501350/100‑1 (MQ=255)
                                                                                  gCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAg              <  2:48013/99‑1 (MQ=255)
                                                                                   cttcttTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTc            <  1:339462/100‑1 (MQ=255)
                                                                                   cttcttTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTc            >  1:1094427/1‑100 (MQ=255)
                                                                                    ttcttTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCa           <  2:196131/100‑1 (MQ=255)
                                                                                      cttTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGa         <  2:2584747/100‑1 (MQ=255)
                                                                                      cttTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGa         >  1:58459/1‑100 (MQ=255)
                                                                                       tttACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAc        >  1:1969799/1‑100 (MQ=255)
                                                                                        ttACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAcc       <  1:1534137/100‑1 (MQ=255)
                                                                                        ttACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAcc       <  2:695129/100‑1 (MQ=255)
                                                                                        ttACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAcc       <  1:1439391/100‑1 (MQ=255)
                                                                                        ttACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAcc       <  1:2668323/100‑1 (MQ=255)
                                                                                           ccTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAat    >  1:52173/1‑100 (MQ=255)
                                                                                           ccTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAAGGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCaa     <  1:625273/99‑1 (MQ=255)
                                                                                            cTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAata   >  2:1901435/1‑100 (MQ=255)
                                                                                            cTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAata   <  1:1242525/100‑1 (MQ=255)
                                                                                                |                                                                                                
TAATAATCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAATAT  >  NC_000913/3840451‑3840643

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