Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,840,547 T→C 64.3% 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,5470TC64.3% 155.5 / 126.4 142intergenic (+39/+2)yicS/nepIputative periplasmic protein/putative transporter
Reads supporting (aligned to +/- strand):  ref base T (20/31);  new base C (47/44);  total (67/75)
Fisher's exact test for biased strand distribution p-value = 1.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

Only 100 of 144 total aligned reads displayed.

CGTAATAATCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTA  >  NC_000913/3840449‑3840645
                                                                                                  |                                                                                                  
  taataaTCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATc                                                                                                 >  2:1604482/1‑100 (MQ=255)
   aataatCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCa                                                                                                >  1:803996/1‑100 (MQ=255)
   aataatCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCa                                                                                                <  1:2178867/100‑1 (MQ=255)
      aatCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGa                                                                                             >  1:782763/1‑100 (MQ=255)
       atCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGAt                                                                                            >  2:772992/1‑100 (MQ=255)
          agagCAATTATTAGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttc                                                                                         <  2:2215867/100‑1 (MQ=255)
           gagCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATttct                                                                                        <  2:2800825/100‑1 (MQ=255)
            agCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATttctt                                                                                       >  1:1988575/1‑100 (MQ=255)
            agCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATttctt                                                                                       <  2:694280/100‑1 (MQ=255)
            agCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttctt                                                                                       <  2:130315/100‑1 (MQ=255)
              cAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCa                                                                                     <  2:1824232/100‑1 (MQ=255)
               aaTTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCAt                                                                                    <  2:1397230/100‑1 (MQ=255)
                attattGGGCAGCGTTAGGTCAGGTAGAATGCCCGGATATGTAAGAACAGCAACACCTACTGACCAACGCTTCTATACCTGCCATCAGGATttctac                                                                                      >  1:2897389/1‑95 (MQ=255)
                attattGGGAGGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCACCGCTTCTTTACCTGATATCAGGATTTCTTCAtt                                                                                   <  1:2229867/100‑1 (MQ=255)
                attattGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCAtt                                                                                   <  1:463323/100‑1 (MQ=255)
                   attGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTc                                                                                >  2:2768636/1‑100 (MQ=255)
                   attGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTc                                                                                <  1:2754546/100‑1 (MQ=255)
                     tGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCAc                                                                              >  1:3022490/1‑100 (MQ=255)
                     tGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCAc                                                                              <  2:612855/100‑1 (MQ=255)
                          aGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTg                                                                         <  1:2768613/100‑1 (MQ=255)
                           gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGc                                                                        <  2:2943782/100‑1 (MQ=255)
                           gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGc                                                                        <  2:812332/100‑1 (MQ=255)
                           gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGc                                                                        <  1:1364446/100‑1 (MQ=255)
                             gTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAg                                                                      <  1:1924837/100‑1 (MQ=255)
                             gTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAg                                                                      <  1:165528/100‑1 (MQ=255)
                               tGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTa                                                                    <  2:2576279/100‑1 (MQ=255)
                                 ggTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACGACTACTGACCAACGCTTCTTTAACTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAAc                                                                  >  1:1202692/1‑100 (MQ=255)
                                 ggTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAAc                                                                  >  2:2693159/1‑100 (MQ=255)
                                  gTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAAcc                                                                 <  2:1707278/100‑1 (MQ=255)
                                  gTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAAcc                                                                 <  1:1657112/100‑1 (MQ=255)
                                  gTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAAcc                                                                 <  1:1268178/100‑1 (MQ=255)
                                  gTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAAcc                                                                 <  2:1446952/100‑1 (MQ=255)
                                   tAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACca                                                                <  2:2538116/100‑1 (MQ=255)
                                      ggTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTAACATGACAACAGGATTTATTCATTTTCACCTTTGCAGTAACcaaca                                                             >  2:1125529/1‑100 (MQ=255)
                                       gTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaacaa                                                            >  1:1869485/1‑100 (MQ=255)
                                       gTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaacaa                                                            >  1:207464/1‑100 (MQ=255)
                                        tAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAAt                                                           >  1:1856127/1‑100 (MQ=255)
                                        tAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAAt                                                           >  1:1562974/1‑100 (MQ=255)
                                          gAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGc                                                         <  2:3022516/100‑1 (MQ=255)
                                           aaTGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGc                                                         <  1:2356009/99‑1 (MQ=255)
                                             tGCCCGGATATGTAAGAACAGCAACAACTACTTACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACACTGCa                                                        >  1:3251594/1‑98 (MQ=255)
                                              gCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCACCAATGCAgcc                                                     >  1:71865/1‑98 (MQ=255)
                                              gCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTGGCAGTAACCAACAATGCAGTc                                                     <  1:1352807/100‑1 (MQ=255)
                                               cccGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCa                                                    >  2:2036260/1‑100 (MQ=255)
                                                ccGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAg                                                   <  2:72681/100‑1 (MQ=255)
                                                 tccgatcttAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGc                                                  <  1:2885162/92‑1 (MQ=255)
                                                  ggATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCa                                                 <  2:1184642/100‑1 (MQ=255)
                                                     tatGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACa                                              <  1:1582027/100‑1 (MQ=255)
                                                      atGTAAGAACAGCAACAACTACTGACCAACGCTTGTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACAt                                             <  1:1787891/100‑1 (MQ=255)
                                                       tGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTACCCAACAATGCAGTCAGCAACATc                                            >  2:1648512/1‑100 (MQ=255)
                                                       tGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATc                                            >  1:2664849/1‑100 (MQ=255)
                                                        gTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCa                                           >  2:2355353/1‑100 (MQ=255)
                                                        gTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCa                                           >  2:2862084/1‑100 (MQ=255)
                                                         tAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCaa                                          >  1:384206/1‑100 (MQ=255)
                                                         tAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCagatcgga                                          >  2:2885185/1‑93 (MQ=255)
                                                          aaGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAAt                                         >  2:1869365/1‑100 (MQ=255)
                                                          aaGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAAt                                         >  1:266317/1‑100 (MQ=255)
                                                           aGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAAtg                                        <  2:531559/100‑1 (MQ=255)
                                                              aCAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTGACCAACAATGCAGTCAGCAACATCAAggt                                       >  1:2832404/1‑98 (MQ=255)
                                                              aCAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGc                                     >  2:3128827/1‑100 (MQ=255)
                                                               cAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGcc                                    >  2:910102/1‑100 (MQ=255)
                                                               cAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGcc                                    >  1:602078/1‑100 (MQ=255)
                                                                 gCAACAACTCCTGACCAACGCTTCTTTACCTGACATCGGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCgg                                  <  1:2280893/100‑1 (MQ=255)
                                                                 gCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCgg                                  >  1:2211158/1‑100 (MQ=255)
                                                                 gCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCgg                                  <  2:896867/100‑1 (MQ=255)
                                                                 gCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCgg                                  >  1:90755/1‑100 (MQ=255)
                                                                 gCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCgg                                  <  2:3365317/100‑1 (MQ=255)
                                                                 gCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCgg                                  >  2:3252685/1‑100 (MQ=255)
                                                                    acaacTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACACCAATGTGCCGGaca                               >  1:529653/1‑100 (MQ=255)
                                                                    acaacTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGaca                               <  2:1856139/100‑1 (MQ=255)
                                                                    acaacTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGaca                               <  1:2206775/100‑1 (MQ=255)
                                                                     caacTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGGCGGacaa                              <  1:335142/100‑1 (MQ=255)
                                                                     caacTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGacaa                              <  1:3185723/100‑1 (MQ=255)
                                                                      aacTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGacaac                             <  2:1102433/100‑1 (MQ=255)
                                                                          actGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCa                         >  2:1541728/1‑100 (MQ=255)
                                                                           ctGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCaa                        >  2:206766/1‑100 (MQ=255)
                                                                           ctGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCaa                        >  2:421928/1‑100 (MQ=255)
                                                                             gatcaccGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCATCAACATCAATGTGCCGGACAACATCAACg                      <  1:836156/94‑1 (MQ=255)
                                                                             gACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACg                      >  2:1089296/1‑100 (MQ=255)
                                                                             gACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACg                      <  2:316508/100‑1 (MQ=255)
                                                                              aCCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACgg                     <  1:1450214/100‑1 (MQ=255)
                                                                              aCCAACGCGTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACgg                     <  2:605690/100‑1 (MQ=255)
                                                                                cAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCg                   >  1:2135997/1‑100 (MQ=255)
                                                                                 aaCGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCg                   >  1:113643/1‑99 (MQ=255)
                                                                                 aaCGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGa                  <  2:585757/100‑1 (MQ=255)
                                                                                  aCGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGaa                 >  1:1098763/1‑100 (MQ=255)
                                                                                   cGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAg                <  2:2211173/100‑1 (MQ=255)
                                                                                   cGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAg                <  1:1016843/100‑1 (MQ=255)
                                                                                   cGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAg                <  1:1286576/100‑1 (MQ=255)
                                                                                    gCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGt               >  1:902957/1‑100 (MQ=255)
                                                                                     cttcttTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTc              <  1:1491106/100‑1 (MQ=255)
                                                                                     cttcttTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTc              <  1:2079804/100‑1 (MQ=255)
                                                                                     cttcttTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTc              <  1:2738618/100‑1 (MQ=255)
                                                                                     cttcttTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTc              <  1:3252657/100‑1 (MQ=255)
                                                                                      ttcttTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCATTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCa             >  1:1094282/1‑100 (MQ=255)
                                                                                         tttACTTGATATCATTATTTCGTCCTGTTCACCTTTGCAGTATCCAACAATGTAGTCAACAACATCAATGGGCCGTACAACATCAACGGCGAACTCATAc          <  1:2790824/100‑1 (MQ=255)
                                                                                         tttACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAc          <  1:331119/100‑1 (MQ=255)
                                                                                          ttACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAcc         <  2:2960517/100‑1 (MQ=255)
                                                                                             ccTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAat      >  1:1384650/1‑100 (MQ=255)
                                                                                                 atatCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGCCAGACCAAtatta  <  2:2278307/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
CGTAATAATCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTA  >  NC_000913/3840449‑3840645

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