Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,840,547 T→C 33.5% 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,5470TC33.5% 179.0 / 157.0 176intergenic (+39/+2)yicS/nepIputative periplasmic protein/putative transporter
Reads supporting (aligned to +/- strand):  ref base T (46/71);  new base C (25/34);  total (71/105)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.80e-01

Only 100 of 184 total aligned reads displayed.

ACGTAATAATCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTAT  >  NC_000913/3840448‑3840646
                                                                                                   |                                                                                                   
 cGTAATAATCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGata                                                                                                    <  2:347636/100‑1 (MQ=255)
  gTAATAATCAGAGCAATTAATGGGAAGCGTTGGGTAAGGTAGAGTGCCCGGATATGTAATAACAGCAACAACTACTGACCAACGCTTCTTTACCTGatat                                                                                                   <  2:2805439/100‑1 (MQ=255)
   taataaTCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATCTGTACCAACAGCCACAACTACTAACCAACGCTTCTTTACCTGATATc                                                                                                  >  2:3952503/1‑100 (MQ=255)
   taataaTCAGAGCAATTATTGGGAAGCGTTGGGTAAGATAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATc                                                                                                  >  2:1010718/1‑100 (MQ=255)
   taataaTCAGAGCAATTATTGGGAAGCGTTAGGTAAGGTAGAATGACCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGata                                                                                                    <  2:499120/98‑1 (MQ=255)
        atCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGAt                                                                                             >  1:689209/1‑100 (MQ=255)
         tCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGAtt                                                                                            <  1:3547319/100‑1 (MQ=255)
           agagCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATttc                                                                                          <  2:3408870/100‑1 (MQ=255)
           agagCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATttc                                                                                          >  1:1132400/1‑100 (MQ=255)
            gagCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTTCTGACCAACGCTTCTTTACCTGATATCAGGATttct                                                                                         <  1:2914159/100‑1 (MQ=255)
             agCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATttctt                                                                                        <  1:730017/100‑1 (MQ=255)
             agCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATAACAGGATttctt                                                                                        <  2:1389623/100‑1 (MQ=255)
             agCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttctt                                                                                        <  1:1289749/100‑1 (MQ=255)
              gCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGTCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATGAGGATttcttc                                                                                       <  2:2169155/100‑1 (MQ=255)
               cAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTAATGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCa                                                                                      >  1:618234/1‑100 (MQ=255)
                acttattGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTAACCAACGCTTCTTTACCTGATAATAGGATctcttc                                                                                       >  2:619135/3‑98 (MQ=255)
                 attattGGGAAGCGTTGGGTAAGGTAGAATACCCTGTTATGTAATAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttcttc                                                                                       <  2:946226/97‑1 (MQ=255)
                  ttattGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCAttt                                                                                   <  1:2338757/100‑1 (MQ=255)
                   tattGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCAtttt                                                                                  <  2:2291162/100‑1 (MQ=255)
                     ttGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCa                                                                                <  1:4047095/100‑1 (MQ=255)
                       gggtcgaggtGGGTAAGGTAGAATGCCCGGTTATGTAAGAACTGCAAAAACTACTCACCCACGCTTCTATCCATGATAGCACGCATTCTCCTTTTTAAc                                                                               >  2:2000630/10‑99 (MQ=255)
                       gggAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCAcc                                                                              <  2:3570944/100‑1 (MQ=255)
                       gggAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCAcc                                                                              >  1:468668/1‑100 (MQ=255)
                       gggAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCCTCTTTACCTGATATCAGGATTTCTTCATTTTCAcc                                                                              >  2:1171049/1‑100 (MQ=255)
                           aGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTg                                                                          <  2:4067159/100‑1 (MQ=255)
                           aGCGCTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTg                                                                          <  1:3507850/100‑1 (MQ=255)
                            gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGc                                                                         <  2:1383222/100‑1 (MQ=255)
                            gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGc                                                                         <  1:1468413/100‑1 (MQ=255)
                            gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGc                                                                         <  1:316909/100‑1 (MQ=255)
                            gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGc                                                                         <  1:2506622/100‑1 (MQ=255)
                              gTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAg                                                                       <  2:2244964/100‑1 (MQ=255)
                              gTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAg                                                                       <  2:2360650/100‑1 (MQ=255)
                              gTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAg                                                                       <  1:3074892/100‑1 (MQ=255)
                              gTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAg                                                                       <  2:194041/100‑1 (MQ=255)
                                  ggTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAAc                                                                   <  1:3428448/100‑1 (MQ=255)
                                  ggTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAAc                                                                   <  1:209393/100‑1 (MQ=255)
                                   gTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAAcc                                                                  >  2:3365622/1‑100 (MQ=255)
                                     aaGGTAGAATGCCCGGATATGTAAGAACAGCCACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaa                                                                <  1:1698162/100‑1 (MQ=255)
                                     aaGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaa                                                                >  2:562819/1‑100 (MQ=255)
                                       ggTAGAACGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaaca                                                              >  1:2160742/1‑100 (MQ=255)
                                         tataaTGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTTCAGTAACcaacaa                                                             <  2:2650247/96‑1 (MQ=255)
                                         tAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaacaa                                                             >  1:2093657/1‑99 (MQ=255)
                                          aGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATg                                                           <  2:412246/100‑1 (MQ=255)
                                           gAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGc                                                          <  2:3835192/100‑1 (MQ=255)
                                             aTGCCCGGATATGTAAGAACAGCAACAACTACTCAACAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAg                                                        <  2:1716985/100‑1 (MQ=255)
                                              tGCCCGGATATGTAAGAACAGCAACACCTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTCCAGTACCccacaa                                                             >  2:2852814/1‑94 (MQ=255)
                                               gCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTc                                                      >  1:466380/1‑100 (MQ=255)
                                               gCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTc                                                      >  2:1321442/1‑100 (MQ=255)
                                               gCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTc                                                      >  1:2244229/1‑100 (MQ=255)
                                                cccGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCa                                                     <  1:3691592/100‑1 (MQ=255)
                                                cccGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCa                                                     <  1:3365642/100‑1 (MQ=255)
                                                   ggATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCa                                                  <  2:1145384/100‑1 (MQ=255)
                                                   ggATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCa                                                  >  2:2337571/1‑100 (MQ=255)
                                                     atatGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAAc                                                <  2:468664/100‑1 (MQ=255)
                                                     atatGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAAc                                                >  2:3215559/1‑100 (MQ=255)
                                                      tatGTAAGAACAGCAACAACTACTGCCCAACGCTTCTTTACCTGATACCAGCATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAAc                                                <  2:18677/99‑1 (MQ=255)
                                                      tatGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACa                                               >  2:303713/1‑100 (MQ=255)
                                                        tGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATc                                             >  2:3248868/1‑100 (MQ=255)
                                                          tAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCaa                                           >  1:446580/1‑100 (MQ=255)
                                                           aaGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAAt                                          >  1:1143369/1‑100 (MQ=255)
                                                           aaGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAAt                                          >  2:1011155/1‑100 (MQ=255)
                                                           aaGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAAt                                          <  1:3270610/100‑1 (MQ=255)
                                                             gAACAGCAACAAGTACTCACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAAtgt                                        <  2:1296794/100‑1 (MQ=255)
                                                                  gCAACAGCTCCTGACCAACGCTTCTTTTCCTGATATCATGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACCTCAATGTGc                                      <  2:2744252/97‑1 (MQ=255)
                                                                  gCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCgg                                   <  2:3782051/100‑1 (MQ=255)
                                                                  gCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCgg                                   <  1:1551656/100‑1 (MQ=255)
                                                                  gCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCgg                                   >  1:2655997/1‑100 (MQ=255)
                                                                  gCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCgg                                   >  1:1769693/1‑100 (MQ=255)
                                                                  gCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCgg                                   <  2:4127016/100‑1 (MQ=255)
                                                                     acaacTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGaca                                <  1:1005838/100‑1 (MQ=255)
                                                                      caacTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGacaa                               <  2:2373884/100‑1 (MQ=255)
                                                                       aacTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGacaac                              <  2:3666015/100‑1 (MQ=255)
                                                                        actactGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGacaaca                             <  1:3208311/100‑1 (MQ=255)
                                                                          tactGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATc                           >  2:2049358/1‑100 (MQ=255)
                                                                           actGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCa                          <  2:293396/100‑1 (MQ=255)
                                                                             tgtccAACGCTTCTTTACATGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACCATGCAGTCAGCAACATCAATGTGCCGGACAACATCa                          <  1:1581675/95‑1 (MQ=255)
                                                                              gACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACg                       <  1:1896544/100‑1 (MQ=255)
                                                                              gACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACg                       >  2:772648/1‑100 (MQ=255)
                                                                               actaaCGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAATgg                      <  1:2464575/97‑1 (MQ=255)
                                                                                ccAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACg                       >  2:1687738/1‑98 (MQ=255)
                                                                                ccAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGc                     <  1:1063758/100‑1 (MQ=255)
                                                                                ccAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGc                     <  2:327690/100‑1 (MQ=255)
                                                                                ccAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGc                     >  2:3223346/1‑100 (MQ=255)
                                                                                 cAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACCTCAATGTGCCTGACAACATCAACGGCg                    >  2:680486/1‑100 (MQ=255)
                                                                                  aaCGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGc                     >  2:3254167/1‑98 (MQ=255)
                                                                                  aaCGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGa                   >  1:2374518/1‑100 (MQ=255)
                                                                                  aaCGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGa                   >  1:2493151/1‑100 (MQ=255)
                                                                                  aaCGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGa                   <  2:1562319/100‑1 (MQ=255)
                                                                                   aCGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCg                    <  2:2508375/98‑1 (MQ=255)
                                                                                    cGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCAGGACAACATCAACGGCGAAg                 <  2:822378/100‑1 (MQ=255)
                                                                                    cGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAg                 <  2:4066840/100‑1 (MQ=255)
                                                                                          tttACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAc           <  1:1948727/100‑1 (MQ=255)
                                                                                           ttACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAcc          <  2:3615369/100‑1 (MQ=255)
                                                                                            tatcTGATATCAGGATTTCTTCCTTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGTACAACATCAACGGCGAAGTCAGAcc          <  2:2577539/96‑1 (MQ=255)
                                                                                            tACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCa         >  1:956256/1‑100 (MQ=255)
                                                                                             aCCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGCCGTCAGCCCaa        >  1:3866714/1‑100 (MQ=255)
                                                                                             aCCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCaa        >  2:1192303/1‑100 (MQ=255)
                                                                                                tGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAatat     <  1:2551720/100‑1 (MQ=255)
                                                                                                 gATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAGGTCAGACCAAtatt    >  1:3528202/1‑100 (MQ=255)
                                                                                                   tatCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAAtattat  >  2:690459/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ACGTAATAATCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATTAT  >  NC_000913/3840448‑3840646

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: