Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,519,225 T→A 15.5% intergenic (+20/+32) 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,2250TA15.5% 368.6 / 53.5 142intergenic (+20/+32)yfeD/gltXDUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (59/61);  new base A (18/4);  total (77/65)
Fisher's exact test for biased strand distribution p-value = 5.05e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

Only 100 of 149 total aligned reads displayed.

GGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTC  >  NC_000913/2519127‑2519323
                                                                                                  |                                                                                                  
ggAAAAAGTCGCGCCATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatta                                                                                                   <  1:3634187/100‑1 (MQ=255)
ggAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatta                                                                                                   >  2:1251829/1‑100 (MQ=255)
   aaaaGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtct                                                                                                >  2:1784439/1‑100 (MQ=255)
         cgcgCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGcc                                                                                          <  1:1149093/100‑1 (MQ=255)
         cgcgCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGcc                                                                                          <  1:2140559/100‑1 (MQ=255)
            gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGtt                                                                                       <  1:37232/100‑1 (MQ=255)
            gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGtt                                                                                       <  1:514878/100‑1 (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:2567339/100‑1 (MQ=255)
                     gTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttttt                                                                              <  1:3232149/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                                                                          <  1:2236474/100‑1 (MQ=255)
                         cGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTacac                                                                          <  2:293869/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)
                               ggAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAAt                                                                    >  2:2555568/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)
                                         aTTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAt                                                          >  1:3100537/1‑100 (MQ=255)
                                           ttCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAttt                                                        >  1:1495505/1‑100 (MQ=255)
                                             ccTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTc                                                      >  2:3257214/1‑100 (MQ=255)
                                                cGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgcg                                                   >  1:3625613/1‑100 (MQ=255)
                                                cGCCCGTCTCCATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTACCTGTTTACACGCTAATTACTGCTGATTTTcgcg                                                   <  2:2511301/100‑1 (MQ=255)
                                                 gCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGt                                                  <  2:1052009/100‑1 (MQ=255)
                                                 gCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGt                                                  >  1:1392924/1‑100 (MQ=255)
                                                  cccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGtt                                                 >  1:1593464/1‑100 (MQ=255)
                                                   ccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTc                                                <  2:39640/100‑1 (MQ=255)
                                                     gTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAg                                              <  1:3315513/100‑1 (MQ=255)
                                                     gTCTCGATATTGACGAATCAGAATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAg                                              <  2:2077001/100‑1 (MQ=255)
                                                        tcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCaa                                           <  1:1527254/100‑1 (MQ=255)
                                                         cGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAAt                                          >  1:3554728/1‑100 (MQ=255)
                                                         cGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAAt                                          >  2:862959/1‑100 (MQ=255)
                                                           atatTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaa                                        <  2:259837/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                                     >  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:2791125/100‑1 (MQ=255)
                                                                  cGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAg                                 <  1:649361/100‑1 (MQ=255)
                                                                  cGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAg                                 >  1:3576479/1‑100 (MQ=255)
                                                                   gAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgc                                <  2:3200019/100‑1 (MQ=255)
                                                                     aTCAGTATAAATAAAAAACGGCAGGATATTAGCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAg                                 <  2:3500978/97‑1 (MQ=255)
                                                                        aGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCttt                           >  1:3777284/1‑100 (MQ=255)
                                                                         gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttg                          >  2:2755110/1‑100 (MQ=255)
                                                                          cATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgt                         >  2:2753283/1‑100 (MQ=255)
                                                                           ataaataaAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                        <  2:1941434/100‑1 (MQ=255)
                                                                            taaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgttg                       >  1:2673291/1‑100 (MQ=255)
                                                                              aataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAt                     <  1:2300613/100‑1 (MQ=255)
                                                                               ataacAAACGGCGGGAGATTTTCTCCTGCAGTTTATCTTTTTACACGCTAAGTACTGATGATTTTCGACATCATCAATAAAATCCATAGCTTTGTTGATa                    >  2:3739392/1‑100 (MQ=255)
                                                                               ataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTGCTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                    >  2:3593645/1‑100 (MQ=255)
                                                                               ataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                    <  1:3521941/100‑1 (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:2238074/5‑100 (MQ=37)
                                                                                agataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATAc                   >  1:2237845/5‑100 (MQ=37)
                                                                                 gataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                    >  2:3278858/4‑98 (MQ=37)
                                                                                 gataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACg                  >  1:684707/4‑100 (MQ=38)
                                                                                 ataaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACg                  >  1:3535523/3‑100 (MQ=255)
                                                                                  ataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGc                 >  1:2355778/3‑100 (MQ=38)
                                                                                   taaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTGACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGc                 >  2:2034094/2‑99 (MQ=38)
                                                                                   aaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCt                >  1:3660960/1‑100 (MQ=255)
                                                                                   aaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCt                <  2:2425481/100‑1 (MQ=255)
                                                                                    aaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTc               <  2:2430743/100‑1 (MQ=255)
                                                                                    aaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGCTACGCTc               >  2:448184/1‑100 (MQ=38)
                                                                                     aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg              >  2:3709904/1‑100 (MQ=39)
                                                                                      aCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGa             >  1:3374987/1‑100 (MQ=39)
                                                                                      aCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGa             >  1:2672965/1‑100 (MQ=39)
                                                                                      aCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGa             >  2:569710/1‑100 (MQ=39)
                                                                                       cGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt            <  2:206626/100‑1 (MQ=255)
                                                                                       cGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt            >  2:1344910/1‑100 (MQ=255)
                                                                                       cGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt            <  1:2889515/100‑1 (MQ=255)
                                                                                       cGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt            >  1:1106237/1‑100 (MQ=39)
                                                                                       cGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACCCTAATTACTGCTGCTTTTCGCGTACGGCAATAAAGTCCAGCGCTTTGTTGATACGCTCTAt            >  2:2354783/1‑100 (MQ=37)
                                                                                        ggCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATa           >  2:2913454/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)
                                                                                                 ttatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTc  <  2:163929/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
GGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTC  >  NC_000913/2519127‑2519323

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: