Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,519,228 C→A 15.0% 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.0% 445.6 / 56.5 165intergenic (+23/+29)yfeD/gltXDUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (71/69);  new base A (9/16);  total (80/85)
Fisher's exact test for biased strand distribution p-value = 1.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 173 total aligned reads displayed.

AAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAC  >  NC_000913/2519129‑2519327
                                                                                                   |                                                                                                   
aaaaaGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtc                                                                                                     <  2:978500/100‑1 (MQ=255)
   aaGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTcc                                                                                                  >  2:2056037/1‑100 (MQ=255)
     gTCGCCCAATTTTTGTCACGTGAAGTGATTCTCAATTTCCTCGCCCTGCTGGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATAGTAtctc                                                                                                   >  1:2518706/1‑97 (MQ=255)
      tCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGc                                                                                               <  2:3369383/100‑1 (MQ=255)
        gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCg                                                                                             <  2:404791/100‑1 (MQ=255)
        gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCg                                                                                             <  2:3860720/100‑1 (MQ=255)
          gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGtt                                                                                           <  2:2502543/100‑1 (MQ=255)
          gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGtt                                                                                           <  1:3170447/100‑1 (MQ=255)
          gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtt                                                                                           <  2:3495985/100‑1 (MQ=255)
          gcAATTTTTGTCACGTGAAGGAATTCGCAATGTCCGCGCCCGGCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGt                                                                                            <  1:277151/99‑1 (MQ=38)
           cAATTTTTGTCTCGTGAAGGAATTGGCAATGTCCGCCCACGTCTCGATATTGACGAATCAGCTTCAATAAAAAACGGCAGGAGATTATCTCGTGCCGttt                                                                                          >  1:658668/1‑100 (MQ=255)
           cAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGttt                                                                                          >  1:1025263/1‑100 (MQ=255)
            aaTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtttt                                                                                         <  2:3660932/100‑2 (MQ=38)
             aTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTAt                                                                                        >  2:2724002/1‑100 (MQ=255)
             aTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttt                                                                                        <  1:4402830/100‑3 (MQ=38)
              tttttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtttttt                                                                                       <  2:787861/100‑4 (MQ=38)
              tttttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtttttt                                                                                       <  2:2039020/100‑4 (MQ=38)
              tttttGTCACGTGAAAGAATTCGCAATTTCCTCGCCCGACTCGATATTGACAAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtttttt                                                                                       <  1:719279/100‑4 (MQ=37)
                tttGTCACGAAAAGTAATTCGCAATTTACTTGCCAGTCTCGATATTGACGAATCAGAATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCt                                                                                      <  1:4395777/99‑1 (MQ=255)
                 ttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttttatt                                                                                    <  2:3008346/100‑7 (MQ=37)
                  tGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCtttt                                                                                   >  2:3227483/1‑100 (MQ=255)
                   gTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttttt                                                                                  <  2:712795/100‑1 (MQ=255)
                      acggggaGGAATTCGCAATTTCCTCGCCCGTCGCGATATTGACGAATCAGCATAAATGAAAAACGGCAGGATGTGATCTCCTGCCGTTTATCTTTTTaaa                                                                               >  1:3082843/7‑98 (MQ=255)
                      aCGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTaca                                                                               <  2:536326/100‑1 (MQ=255)
                                aTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAc                                                                     <  2:420513/100‑1 (MQ=255)
                                 ttCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAct                                                                    >  1:1247129/1‑100 (MQ=255)
                                 ttCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAct                                                                    >  1:2733538/1‑100 (MQ=255)
                                   cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                                  >  2:1241917/1‑100 (MQ=255)
                                   cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                                  >  2:443640/1‑100 (MQ=255)
                                    gCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgct                                                                 <  2:2297298/100‑1 (MQ=255)
                                    gCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgct                                                                 >  2:3254938/1‑100 (MQ=255)
                                    gCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgct                                                                 >  1:4201136/1‑100 (MQ=255)
                                     cAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgctg                                                                >  2:2578828/1‑100 (MQ=255)
                                       aTTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAt                                                              <  1:300191/100‑1 (MQ=255)
                                        tttCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAtt                                                             >  1:1211605/1‑100 (MQ=255)
                                            cTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcg                                                         <  2:4400671/100‑1 (MQ=255)
                                             tCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcg                                                         <  1:2056029/99‑1 (MQ=255)
                                              cgacggTCTCGATGTTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgcg                                                       <  1:855456/95‑1 (MQ=255)
                                              cGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgcg                                                       <  1:2433965/100‑1 (MQ=255)
                                              cGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGGTTATCTTTTTACACGCTAATTACTGCTGATTTTcgcg                                                       <  1:2347204/100‑1 (MQ=255)
                                               gCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGt                                                      >  1:397838/1‑100 (MQ=255)
                                                cccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGtt                                                     >  1:2207842/1‑100 (MQ=255)
                                                 ccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTc                                                    <  2:3890371/100‑1 (MQ=255)
                                                    tctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGc                                                 <  1:14458/100‑1 (MQ=255)
                                                    tctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGc                                                 <  1:2918887/100‑1 (MQ=255)
                                                    tctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGc                                                 <  2:3201083/100‑1 (MQ=255)
                                                       cgctatTGCCGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAAt                                              >  1:1134822/4‑100 (MQ=255)
                                                          tatTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaa                                           >  1:1629767/1‑100 (MQ=255)
                                                          tatTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaa                                           >  2:2476547/1‑100 (MQ=255)
                                                           atTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaaa                                          >  1:1916405/1‑100 (MQ=255)
                                                            ttGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAAt                                         <  2:4302787/100‑1 (MQ=255)
                                                              gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                       <  1:783415/100‑1 (MQ=255)
                                                              gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                       >  2:1822616/1‑100 (MQ=255)
                                                                cGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAg                                     <  2:1718227/100‑1 (MQ=255)
                                                                 gAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgc                                    <  2:4137163/100‑1 (MQ=255)
                                                                  aaTCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgcg                                   >  1:4335893/1‑100 (MQ=255)
                                                                   aTCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCTCGTTCAGCAATAAAATCCAgcgc                                  >  1:357958/1‑100 (MQ=255)
                                                                    tCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCt                                 >  2:3701990/1‑100 (MQ=255)
                                                                      aGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTTCTGATTTTCGCGATCAGCAATAAAATCCAGCGCtt                                >  1:4147927/1‑99 (MQ=255)
                                                                      aGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTATCGCGTTCAGCAATAAAATCCAGCGCttt                               >  1:2652371/1‑100 (MQ=255)
                                                                       gtataaataaAAAACGGCAGCATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTGTCGCGTTCAGCAAGAAAATCCAGCGCtt                                <  1:3381460/96‑1 (MQ=255)
                                                                       gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttg                              >  2:3612405/1‑100 (MQ=255)
                                                                        cATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgt                             <  1:3156245/100‑1 (MQ=255)
                                                                         ataaataaAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                            <  2:2652384/100‑1 (MQ=255)
                                                                         ataaataaAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                            >  2:3293800/1‑100 (MQ=255)
                                                                           aaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATTCAGCGCTttgttg                           >  1:3636091/1‑99 (MQ=255)
                                                                           aaaagataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGa                          >  2:1288135/8‑100 (MQ=37)
                                                                            aataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAt                         >  2:3140350/1‑100 (MQ=255)
                                                                              taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAtgc                       >  1:980868/1‑98 (MQ=255)
                                                                              taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATAc                       >  2:2173582/1‑100 (MQ=255)
                                                                              agataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATAc                       >  2:1223209/5‑100 (MQ=37)
                                                                               gataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGGTTTTTGCGTTTGGAACTAAAATCCAGCGCTTTGTTGATa                        >  1:3793120/4‑98 (MQ=21)
                                                                               aaaaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACg                      >  1:1951690/1‑100 (MQ=255)
                                                                                aaaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATTCAGCGCTTTGTTGATCCGc                     >  1:946578/1‑100 (MQ=255)
                                                                                aaaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGc                     >  1:3896678/1‑100 (MQ=255)
                                                                                 aaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCt                    <  2:1688865/100‑1 (MQ=255)
                                                                                   aaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTc                   <  1:2883481/99‑1 (MQ=255)
                                                                                   aaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg                  >  1:2182877/1‑100 (MQ=255)
                                                                                   aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg                  >  2:2740661/1‑100 (MQ=39)
                                                                                    ncGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGCTTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACTCTCGa                 >  1:2071043/2‑100 (MQ=38)
                                                                                    aCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGa                 >  2:4063269/1‑100 (MQ=39)
                                                                                     cGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt                >  2:662112/1‑100 (MQ=255)
                                                                                     cGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt                <  1:3006673/100‑1 (MQ=255)
                                                                                     cGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt                >  1:669909/1‑100 (MQ=255)
                                                                                       gCAGGATATTATCTCCTGCCGTTTCTCTTTTTACACGCTAATTACTGCTGATTTTCGCTTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATa               >  1:238905/1‑99 (MQ=255)
                                                                                       gCAGGATAGTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCTTTCAGCAATAAAATCCAGTGCTTTGTTGATACTTTCTATAc              >  1:1454622/1‑100 (MQ=255)
                                                                                        cAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACt             >  1:1582335/1‑100 (MQ=255)
                                                                                           gATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCg          >  1:3374104/1‑100 (MQ=255)
                                                                                           gATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCg          <  1:1648718/100‑1 (MQ=255)
                                                                                             tattatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCggg        <  2:2382228/100‑1 (MQ=255)
                                                                                             tattatCTCCTGCCGTTTATCTTTTTACACGCCAATTACGGCCGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCggg        <  1:3085964/100‑1 (MQ=255)
                                                                                              attatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCgggn       <  1:2071162/100‑2 (MQ=255)
                                                                                               ttatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTc      <  1:1193044/100‑1 (MQ=255)
                                                                                               ctatCTCCTGCCGTTTATCTTTCTACACGCTAATTACTGCTTATTTTCGCTCTCAGCTATACAATCCACTGCTTTGTTTATACGCTCGATACTGCGGGTc      <  1:4318750/99‑1 (MQ=255)
                                                                                                tatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCt     >  2:1496518/1‑100 (MQ=255)
                                                                                                tatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCt     <  1:2611621/100‑1 (MQ=255)
                                                                                                 atCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCtt    >  2:2569695/1‑100 (MQ=255)
                                                                                                  tctcCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCt     <  1:1692515/98‑1 (MQ=255)
                                                                                                   ctcCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTGGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAc  <  1:2823310/100‑1 (MQ=255)
                                                                                                   ctcCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAc  <  2:924078/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCTTAC  >  NC_000913/2519129‑2519327

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: