Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,519,228 C→A 15.3% intergenic (+23/+29) yfeD → / ← gltX DUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,519,2280CA15.3% 396.9 / 57.7 144intergenic (+23/+29)yfeD/gltXDUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (58/64);  new base A (18/4);  total (76/68)
Fisher's exact test for biased strand distribution p-value = 4.52e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-02

Only 100 of 148 total aligned reads displayed.

AAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAC  >  NC_000913/2519130‑2519327
                                                                                                  |                                                                                                   
aaaaGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtct                                                                                                    >  2:1784439/1‑100 (MQ=255)
aaaaGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtct                                                                                                    <  1:3248673/100‑1 (MQ=255)
      cgcgCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGcc                                                                                              <  1:1149093/100‑1 (MQ=255)
      cgcgCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGcc                                                                                              <  1:2140559/100‑1 (MQ=255)
        cgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGt                                                                                            >  1:587401/1‑100 (MQ=255)
         gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGtt                                                                                           <  1:3609360/100‑1 (MQ=255)
         gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGtt                                                                                           >  1:3737619/1‑100 (MQ=255)
         gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtt                                                                                           <  2:1685258/100‑1 (MQ=39)
          cAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGttt                                                                                          <  1:1852455/100‑1 (MQ=255)
              ttttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGTTTTTca                                                                                      <  1:2578095/100‑2 (MQ=38)
               tttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCtt                                                                                     <  1:899567/100‑1 (MQ=255)
                  gTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttttt                                                                                  <  1:3232149/100‑1 (MQ=255)
                  gTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttttt                                                                                  <  1:2567339/100‑1 (MQ=255)
                   taaCGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGACAGGATATTATCTCCTGCCGTTTATCTTTTTa                                                                                 <  2:1399854/98‑1 (MQ=255)
                    cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGGAGGATATTATCTCCTGCCGTTTATCTTTTTac                                                                                >  2:593719/1‑100 (MQ=255)
                    cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTac                                                                                >  1:1840619/1‑100 (MQ=255)
                    cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTac                                                                                >  2:480094/1‑100 (MQ=255)
                     acttGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTaca                                                                               <  2:3321207/97‑1 (MQ=255)
                     aCGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTaca                                                                               >  1:259837/1‑100 (MQ=255)
                      cGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTacac                                                                              <  2:293869/100‑1 (MQ=255)
                      cGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTacac                                                                              <  1:2236474/100‑1 (MQ=255)
                       gTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGTATATTATCTCCTGCCGTTTATCTTTTTACACg                                                                             <  2:896081/100‑1 (MQ=255)
                       gTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACg                                                                             <  1:921188/100‑1 (MQ=255)
                        tGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGc                                                                            <  2:3468324/100‑1 (MQ=255)
                        tGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGc                                                                            >  1:2332734/1‑100 (MQ=255)
                         gAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCt                                                                           <  2:3023410/100‑1 (MQ=255)
                          aaGGAATTCGCAATTTCCTCGCCCGTCTCTATATTGACCAATCAGCATAAATAAAAAACGGCAGGATAATATCTCCTGCCGTTTATCTTTTTACACGCt                                                                           <  2:1588480/99‑1 (MQ=255)
                           aGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTaa                                                                         <  2:2830271/100‑1 (MQ=255)
                            ggAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGAATATTATCTCCTGCCGTTTATCTTTTTACACGCTAAt                                                                        >  1:3396465/1‑100 (MQ=255)
                               aTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAc                                                                     <  1:1251838/100‑1 (MQ=255)
                                 tCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActg                                                                   >  1:1779290/1‑100 (MQ=255)
                                  tgCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                                  <  1:2946611/99‑1 (MQ=255)
                                  cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                                  >  2:2987921/1‑100 (MQ=255)
                                   gCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAACAACGGCCGGATCTTATCTCCTGCCGTTTATCTTTTTACACGCTAATTCcttct                                                                 >  2:2341749/1‑100 (MQ=255)
                                   gCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgct                                                                 >  1:738437/1‑100 (MQ=255)
                                     aaTTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACGGCTGa                                                               <  2:543545/100‑1 (MQ=255)
                                        ttCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAttt                                                            >  1:1495505/1‑100 (MQ=255)
                                          ccTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTc                                                          >  2:3257214/1‑100 (MQ=255)
                                             cGCCCGTCTCCATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTACCTGTTTACACGCTAATTACTGCTGATTTTcgcg                                                       <  2:2511301/100‑1 (MQ=255)
                                              gCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGt                                                      >  1:1392924/1‑100 (MQ=255)
                                              gCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGt                                                      <  2:1052009/100‑1 (MQ=255)
                                               cccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGtt                                                     >  1:1593464/1‑100 (MQ=255)
                                                ccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTc                                                    <  2:39640/100‑1 (MQ=255)
                                                  gTCTCGATATTGACGAATCAGAATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAg                                                  <  2:2077001/100‑1 (MQ=255)
                                                     tcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCaa                                               <  1:1527254/100‑1 (MQ=255)
                                                         tatTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaa                                           >  1:1166974/1‑100 (MQ=255)
                                                          atTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaaa                                          >  1:1239646/1‑100 (MQ=255)
                                                          atTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCACCAATaaaa                                          >  2:1707510/1‑100 (MQ=255)
                                                           ttGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAAt                                         <  1:2946743/100‑1 (MQ=255)
                                                           ttGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAAt                                         >  2:2138480/1‑100 (MQ=255)
                                                           ttGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAAt                                         >  1:1344906/1‑100 (MQ=255)
                                                            tGACGAATAAGTACATATAAAGAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTGACTCGTTAATGACTGCTGATGTTCGCGTTCa                                                   <  2:123313/89‑1 (MQ=255)
                                                             gACGAATCAGCATAAATAAAAAACGGCATGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                       <  2:2584324/100‑1 (MQ=255)
                                                             gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTCCACCCTAATCACTGCTGATTTTCGCGGTCAGCAATAAACTcc                                       >  1:3756842/1‑100 (MQ=255)
                                                             gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                       <  2:417364/100‑1 (MQ=255)
                                                             gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                       <  2:2791125/100‑1 (MQ=255)
                                                               cGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAg                                     >  1:3576479/1‑100 (MQ=255)
                                                               cGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAg                                     <  1:649361/100‑1 (MQ=255)
                                                                gAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgc                                    <  2:3200019/100‑1 (MQ=255)
                                                                  aTCAGTATAAATAAAAAACGGCAGGATATTAGCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAg                                     <  2:3500978/97‑1 (MQ=255)
                                                                      gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttg                              >  2:2755110/1‑100 (MQ=255)
                                                                       cATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgt                             >  2:2753283/1‑100 (MQ=255)
                                                                        ataaataaAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                            <  2:1941434/100‑1 (MQ=255)
                                                                        ataaataaAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                            <  1:417922/100‑1 (MQ=255)
                                                                         taaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgttg                           >  1:2673291/1‑100 (MQ=255)
                                                                         taaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgttg                           >  1:3652710/1‑100 (MQ=255)
                                                                           aataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAt                         <  1:2300613/100‑1 (MQ=255)
                                                                            ataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTGCTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                        >  2:3593645/1‑100 (MQ=255)
                                                                            ataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                        <  2:1105700/100‑1 (MQ=255)
                                                                            aagataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                        >  1:447643/6‑100 (MQ=37)
                                                                             taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATAc                       <  2:1970475/100‑1 (MQ=255)
                                                                             agataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATAc                       >  1:2237845/5‑100 (MQ=37)
                                                                             agataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATAc                       >  1:2238074/5‑100 (MQ=37)
                                                                              gataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                        >  2:3278858/4‑98 (MQ=37)
                                                                              ataaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACg                      >  1:3535523/3‑100 (MQ=255)
                                                                               ataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTTTTTTTTTTTTACTGATGATATTCGCGGTTAGCAATAAAATCTAGCGTTttgttcgtgagc                     >  2:2568828/3‑93 (MQ=9)
                                                                               ataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGc                     >  1:2355778/3‑100 (MQ=38)
                                                                                taaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTGACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGc                     >  2:2034094/2‑99 (MQ=38)
                                                                                aaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCt                    <  2:2425481/100‑1 (MQ=255)
                                                                                aaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCt                    >  1:3660960/1‑100 (MQ=255)
                                                                                  aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg                  >  2:3709904/1‑100 (MQ=39)
                                                                                   aCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATCCGCTCGa                 >  2:834289/1‑100 (MQ=38)
                                                                                   aCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGa                 >  1:3374987/1‑100 (MQ=39)
                                                                                   aCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGa                 >  1:2672965/1‑100 (MQ=39)
                                                                                    cGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt                >  2:1344910/1‑100 (MQ=255)
                                                                                    cGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt                <  2:206626/100‑1 (MQ=255)
                                                                                    cGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCTTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt                >  2:2686142/1‑100 (MQ=39)
                                                                                    cGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt                >  1:1106237/1‑100 (MQ=39)
                                                                                    cGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACCCTAATTACTGCTGCTTTTCGCGTACGGCAATAAAGTCCAGCGCTTTGTTGATACGCTCTAt                >  2:2354783/1‑100 (MQ=37)
                                                                                     ggCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATa               >  2:2913454/1‑100 (MQ=255)
                                                                                      gCAGGATATTATCTCCTGCCGTTTATCTTTATACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATCGCCTCGATAc              >  2:2510995/1‑100 (MQ=255)
                                                                                       cAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACt             >  1:2602270/1‑100 (MQ=255)
                                                                                        aGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATAccg            >  2:2307607/1‑98 (MQ=255)
                                                                                         ggATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCGCGATACTGc           >  1:2355991/1‑100 (MQ=255)
                                                                                           atatTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCgg         <  1:2683811/100‑1 (MQ=255)
                                                                                              ttatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTc      <  2:163929/100‑1 (MQ=255)
                                                                                              ttatCTCCTGCCGTGTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTc      <  2:112859/100‑1 (MQ=255)
                                                                                               tatCTCCTGCCCTTTATCGTTTTACACGCTAATTACTGCTGATTTTCGCGTTAAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTc      <  1:3551562/99‑1 (MQ=255)
                                                                                                atCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCtt    <  2:2374716/100‑1 (MQ=255)
                                                                                                 tctcCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTa   >  2:932287/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
AAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAC  >  NC_000913/2519130‑2519327

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: