Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,519,222 T→G 16.3% intergenic (+17/+35) 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,2220TG16.3% 361.5 / 62.0 140intergenic (+17/+35)yfeD/gltXDUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (57/60);  new base G (19/4);  total (76/64)
Fisher's exact test for biased strand distribution p-value = 2.85e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

Only 100 of 151 total aligned reads displayed.

GCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGG  >  NC_000913/2519124‑2519320
                                                                                                  |                                                                                                  
gaaGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGata                                                                                                   >  1:2575072/3‑100 (MQ=255)
 cAGGAAAAAGTCGCGCAATTTTTTTTTTTTTTTTTTATTCGCAATTTCCTCGCGCGTCTTTATGTTGTCGATTCAGCATAAATAAAAAGTTGTAGGatat                                                                                                  >  2:1243153/1‑100 (MQ=255)
 cAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGata                                                                                                   >  2:1531797/1‑99 (MQ=255)
   ggAAAAAGTCGCGCCATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatta                                                                                                <  1:3634187/100‑1 (MQ=255)
   ggAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatta                                                                                                <  2:3734391/100‑1 (MQ=255)
   ggAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatta                                                                                                >  1:708367/1‑100 (MQ=255)
      aaaaGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtct                                                                                             >  2:1784439/1‑100 (MQ=255)
            cgcgCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGcc                                                                                       <  1:1149093/100‑1 (MQ=255)
             gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCg                                                                                      <  2:3583502/100‑1 (MQ=255)
             gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACCGCAGGATATTATCTCCTGCCg                                                                                      >  2:1562151/1‑100 (MQ=255)
              cgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGt                                                                                     >  1:587401/1‑100 (MQ=255)
              cgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGt                                                                                     <  1:2342535/100‑1 (MQ=39)
               gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGtt                                                                                    <  1:3609360/100‑1 (MQ=255)
               gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGtt                                                                                    <  1:37232/100‑1 (MQ=255)
               gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGtt                                                                                    >  1:3737619/1‑100 (MQ=255)
               gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGtt                                                                                    <  1:514878/100‑1 (MQ=255)
               gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGtt                                                                                    <  1:1869116/100‑1 (MQ=255)
               gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtt                                                                                    <  2:1685258/100‑1 (MQ=39)
                cAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGttt                                                                                   <  1:1852455/100‑1 (MQ=255)
                cAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGttt                                                                                   >  2:1947837/1‑100 (MQ=255)
                 aaTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtt                                                                                    <  2:3263136/98‑1 (MQ=255)
                   tttttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGATGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATc                                                                                <  1:2243605/100‑1 (MQ=255)
                   tttttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATc                                                                                >  2:2771078/1‑100 (MQ=255)
                    ttttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGTTTTTca                                                                               <  1:2578095/100‑2 (MQ=38)
                     tttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCtt                                                                              <  1:899567/100‑1 (MQ=255)
                     tttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCt                                                                               >  2:2365862/1‑99 (MQ=255)
                       tGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCtttt                                                                            <  2:3058351/100‑1 (MQ=255)
                        gTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttttt                                                                           <  1:3232149/100‑1 (MQ=255)
                          cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGGAGGATATTATCTCCTGCCGTTTATCTTTTTac                                                                         >  2:593719/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)
                             gTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAGAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACg                                                                      <  2:2082159/100‑1 (MQ=255)
                             gTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGTATATTATCTCCTGCCGTTTATCTTTTTACACg                                                                      <  2:896081/100‑1 (MQ=255)
                              tGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGc                                                                     >  1:2332734/1‑100 (MQ=255)
                                   gAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAAtt                                                                >  2:1768665/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                                                           >  1:1578728/1‑100 (MQ=255)
                                        cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                           <  1:1218863/100‑1 (MQ=255)
                                        cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                           >  2:2987921/1‑100 (MQ=255)
                                         gCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAACAACGGCCGGATCTTATCTCCTGCCGTTTATCTTTTTACACGCTAATTCcttct                                                          >  2:2341749/1‑100 (MQ=255)
                                         gCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgca                                                          >  2:2946590/1‑99 (MQ=255)
                                           aaTTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACGGCTGa                                                        <  2:543545/100‑1 (MQ=255)
                                            aTTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAt                                                       >  1:3100537/1‑100 (MQ=255)
                                              ttCCTCGCCCGTCTCGATATTGACGTATCAGCAGAAATAAACAACGGCAGGATATTATCGCCTGCCGTTTAGCTTTTTACACGCTAATTActgctg                                                         <  1:1811085/96‑1 (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)
                                                      ccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTc                                             <  2:39640/100‑1 (MQ=255)
                                                        gTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAg                                           <  1:3315513/100‑1 (MQ=255)
                                                           tcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCaa                                        >  1:900995/1‑100 (MQ=255)
                                                           tcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCaa                                        <  1:1527254/100‑1 (MQ=255)
                                                           tcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCa                                         >  2:3323807/1‑99 (MQ=255)
                                                               tatTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaa                                    >  1:2366105/1‑100 (MQ=255)
                                                                atTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaaa                                   >  1:1239646/1‑100 (MQ=255)
                                                                 ttGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAAt                                  >  2:2138480/1‑100 (MQ=255)
                                                                  tGACGAATAAGTACATATAAAGAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTGACTCGTTAATGACTGCTGATGTTCGCGTTCa                                            <  2:123313/89‑1 (MQ=255)
                                                                   gACGAATCAGCATAAATAAAAAACGGCATGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                <  2:2584324/100‑1 (MQ=255)
                                                                   gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                <  2:417364/100‑1 (MQ=255)
                                                                   gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                <  2:2791125/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)
                                                                          cAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCtt                         <  2:1332268/100‑1 (MQ=255)
                                                                           aGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCttt                        >  1:3777284/1‑100 (MQ=255)
                                                                             cATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgt                      >  2:2753283/1‑100 (MQ=255)
                                                                              ataaataaAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAACTCCAGCGCTttgtt                     <  2:2833133/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)
                                                                                  ataaaAAACGGCAGGATATTGTCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                 >  2:2331029/1‑100 (MQ=255)
                                                                                  ataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                 <  1:3521941/100‑1 (MQ=255)
                                                                                   taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATAc                <  2:1970475/100‑1 (MQ=255)
                                                                                   taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTGTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTCTTGATAc                >  2:2047203/1‑100 (MQ=255)
                                                                                   agataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATAc                >  1:2238074/5‑100 (MQ=37)
                                                                                    gataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACg               >  1:684707/4‑100 (MQ=38)
                                                                                    ataaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACg               >  1:3535523/3‑100 (MQ=255)
                                                                                    aaaaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACg               <  1:1488975/100‑1 (MQ=255)
                                                                                     ataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTTTTTTTTTTTTACTGATGATATTCGCGGTTAGCAATAAAATCTAGCGTTttgttcgtgagc              >  2:2568828/3‑93 (MQ=9)
                                                                                     aaaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGc              >  1:1332131/1‑100 (MQ=255)
                                                                                      aaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCt             <  2:2425481/100‑1 (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:206626/100‑1 (MQ=255)
                                                                                          cGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt         <  1:2889515/100‑1 (MQ=255)
                                                                                          cGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt         >  2:1344910/1‑100 (MQ=255)
                                                                                          cGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCTTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt         >  2:2686142/1‑100 (MQ=39)
                                                                                          cGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt         >  1:1106237/1‑100 (MQ=39)
                                                                                           ggCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATa        >  1:223722/1‑100 (MQ=255)
                                                                                           ggCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATa        >  2:2913454/1‑100 (MQ=255)
                                                                                            gCAGGATATTATCTCCTGCCGTTTATCTTTATACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATCGCCTCGATAc       >  2:2510995/1‑100 (MQ=255)
                                                                                             cAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACt      >  1:2602270/1‑100 (MQ=255)
                                                                                               ggATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCGCGATACTGc    >  1:2355991/1‑100 (MQ=255)
                                                                                                 atatTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCgg  >  2:1683221/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
GCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGG  >  NC_000913/2519124‑2519320

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: