Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,519,225 T→A 18.1% 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,2250TA18.1% 294.0 / 51.7 126intergenic (+20/+32)yfeD/gltXDUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (47/56);  new base A (12/11);  total (59/67)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 134 total aligned reads displayed.

GGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCT  >  NC_000913/2519127‑2519324
                                                                                                  |                                                                                                   
ggAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatta                                                                                                    >  1:1682982/1‑100 (MQ=255)
   aaaaGTCGCGCAAGTTTAGTCACGTGAACGAAATCGAAAAAAAATCGCCCGTCTCGATATTGACGAAAAAAAAAAAAAAAAAAACGGCAGGATATTAtc                                                                                                  <  2:1698003/99‑1 (MQ=255)
     aaGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTcc                                                                                               >  2:2561437/1‑100 (MQ=255)
        tCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGc                                                                                            <  2:10242/100‑1 (MQ=255)
        tCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGc                                                                                            <  2:1165054/100‑1 (MQ=255)
         cgcgCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGcc                                                                                           >  1:2440770/1‑100 (MQ=255)
         cgcgCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGcc                                                                                           >  2:1039630/1‑100 (MQ=255)
         cgcgCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGcc                                                                                           <  2:1372971/100‑1 (MQ=255)
          gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCg                                                                                          <  1:3496101/100‑1 (MQ=255)
            gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtt                                                                                        <  2:2991318/100‑1 (MQ=39)
            gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtt                                                                                        <  1:2980190/100‑1 (MQ=39)
             cAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGttt                                                                                       >  2:3455134/1‑100 (MQ=255)
             cAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttt                                                                                       <  2:141002/100‑1 (MQ=38)
             cAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttt                                                                                       <  1:3198160/100‑1 (MQ=38)
             cAATTTTGGTCACGTGAAGGAATTCGCAACTTCCGCGCCCGTGTGGATATTGACAAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGtt                                                                                        <  2:1722588/99‑1 (MQ=255)
              aaTTTTGGTCACGTGAAGGAATTCGCAATGTCCTCGCGTGGCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGtgt                                                                                       <  2:2612325/99‑3 (MQ=255)
               aTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTAt                                                                                     >  1:2252942/1‑100 (MQ=255)
               aTTTTTGTCACGTGAAGGAATTCGCAAGTTACTCGACTGTCTCGATATCGACGAATCAGCATAAATAAAAAACGGCAGGACATAATATCCTGCCGtt                                                                                        <  2:1488413/97‑1 (MQ=37)
                tttttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGtttt                                                                                      <  2:2971985/98‑2 (MQ=38)
                  tttGTCACGTGAAGGATTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttttat                                                                                  <  2:1004087/100‑6 (MQ=37)
                  tttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttttat                                                                                  <  1:782715/100‑6 (MQ=37)
                     gTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttttt                                                                               <  1:2956971/100‑1 (MQ=255)
                     gTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttttt                                                                               >  2:751815/1‑100 (MQ=255)
                      tCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTa                                                                              <  1:2014268/100‑1 (MQ=255)
                       cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTan                                                                             <  1:2838378/100‑2 (MQ=255)
                       cACGTGAAGGAATTAGCAATTTCCTCGCCCGTCTAGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTac                                                                             >  1:742843/1‑100 (MQ=255)
                         cGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTacac                                                                           >  2:1813077/1‑100 (MQ=255)
                         cGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTaca                                                                            <  1:3642652/99‑1 (MQ=255)
                           tGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGc                                                                         >  1:2054937/1‑100 (MQ=255)
                           tGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGc                                                                         >  1:3353038/1‑100 (MQ=255)
                           tGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGc                                                                         >  2:185538/1‑100 (MQ=255)
                           tGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGc                                                                         >  1:1348429/1‑100 (MQ=255)
                            gAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCt                                                                        >  1:241481/1‑100 (MQ=255)
                             aaGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTa                                                                       <  1:3470680/100‑1 (MQ=255)
                                gAATTCGCAATTTGCTCGCCCGTCGCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAAt                                                                     <  2:1580023/99‑1 (MQ=255)
                                 aaTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTa                                                                   <  1:2884303/100‑1 (MQ=255)
                                  aTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAc                                                                  >  1:1283602/1‑100 (MQ=255)
                                   ttCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAct                                                                 >  1:124547/1‑100 (MQ=255)
                                   ttCGCAATTTCCTCACCCGTCTCCATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAc                                                                  <  2:3409512/99‑1 (MQ=255)
                                    tCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActg                                                                <  1:1385043/100‑1 (MQ=255)
                                    tCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActg                                                                <  2:2589522/100‑1 (MQ=255)
                                    tCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActg                                                                <  2:1897289/100‑1 (MQ=255)
                                    tCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActg                                                                <  1:1790454/100‑1 (MQ=255)
                                     cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                               >  1:867065/1‑100 (MQ=255)
                                     cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACCGCAGGATATTATCTCCTGCCGTATATCTTTTTACAAGCTAATTActg                                                                >  2:1984084/1‑99 (MQ=255)
                                      gCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgct                                                              >  2:3116814/1‑100 (MQ=255)
                                       cAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgctg                                                             <  2:400899/100‑1 (MQ=255)
                                       cAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgctg                                                             >  2:737491/1‑100 (MQ=255)
                                        aaTTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGa                                                            <  1:3326901/100‑1 (MQ=255)
                                           ttCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAttt                                                         <  2:3286004/100‑1 (MQ=255)
                                            tCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAtttt                                                        <  2:1283594/100‑1 (MQ=255)
                                             ccTCGCCCGTCTGGATATAGACGAATCAGCATAAAGAAAAACCGGCCGGGTATTATATGCTGCCGCGTATATTTTTTCGCGGTAATAActtcggtatttc                                                       >  2:2105483/1‑90 (MQ=255)
                                             ccTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAttttn                                                       <  2:3171753/100‑2 (MQ=255)
                                              cTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcg                                                      <  1:2115534/100‑1 (MQ=255)
                                               tagCCCGTCTAGATATTGACGAATCAACATAAATAAATAACTGCAGGATATTATCTCCTGCCGTTTATCTTTTTACTCGCTAATTACTGCTGATTTTcg                                                      <  2:2973805/97‑1 (MQ=255)
                                                cGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgcg                                                    <  2:2002873/100‑1 (MQ=255)
                                                cGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCGTTTTACACGCTAATTACTGCTGATTTTcg                                                      <  2:3341363/98‑1 (MQ=255)
                                                  cccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGtt                                                  >  1:1071259/1‑100 (MQ=255)
                                                    cGTCTCTATAGTGACTAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCa                                                <  1:1576032/100‑1 (MQ=255)
                                                     gTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAg                                               <  2:460516/100‑1 (MQ=255)
                                                      tctcGATATTGACGAATCAGCATAAATAAAAAACGGCCGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAgn                                              <  1:1730025/100‑2 (MQ=255)
                                                        tcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTCTCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGc                                              <  1:1103444/98‑1 (MQ=255)
                                                          gagatTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCATGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATa                                          <  2:500811/97‑1 (MQ=255)
                                                            tatTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaa                                        >  2:654862/1‑100 (MQ=255)
                                                              ttGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAAt                                      >  2:2051511/1‑100 (MQ=255)
                                                                gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                    <  1:2181033/100‑1 (MQ=255)
                                                                 aCGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCa                                   <  1:2656469/100‑1 (MQ=255)
                                                                   gAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgc                                 <  1:2712392/100‑1 (MQ=255)
                                                                      tCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTATTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgcgc                               <  2:2324208/99‑1 (MQ=255)
                                                                        aGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCtt                             >  2:1191065/1‑99 (MQ=255)
                                                                         gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttg                           >  1:288885/1‑100 (MQ=255)
                                                                          cATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgt                          >  1:536536/1‑100 (MQ=255)
                                                                           ataaataaAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTTGCGTTCAGCAATAAAATATAGCGCTttgtt                         >  2:3303804/1‑100 (MQ=255)
                                                                           ataaataaAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                         <  1:3210491/100‑1 (MQ=255)
                                                                           ataaataaAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                         <  1:732285/100‑1 (MQ=255)
                                                                            taaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgttg                        >  1:499138/1‑100 (MQ=255)
                                                                             aaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGa                       >  2:2128341/1‑100 (MQ=255)
                                                                              aataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAt                      >  2:2131172/1‑100 (MQ=255)
                                                                               ataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                     <  1:1039638/100‑1 (MQ=255)
                                                                               ataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                     <  1:562303/100‑1 (MQ=255)
                                                                               aagataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                     >  1:2538841/6‑100 (MQ=37)
                                                                                taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATAc                    >  1:2683959/1‑100 (MQ=255)
                                                                                  aaaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGc                  >  2:1457961/1‑100 (MQ=255)
                                                                                   aaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAt                      >  2:2204632/1‑95 (MQ=255)
                                                                                     aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg               >  1:1544748/1‑100 (MQ=39)
                                                                                     aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg               >  1:1856923/1‑100 (MQ=39)
                                                                                     aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg               >  1:3184177/1‑100 (MQ=39)
                                                                                     aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg               >  2:3490642/1‑100 (MQ=39)
                                                                                      aCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGCTTTTCGCGTTCAGCAATAACATCCAGCGCTTTGTTGATACGCTc                >  2:3101274/1‑98 (MQ=38)
                                                                                      aCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATaagttc                >  2:1146005/1‑93 (MQ=38)
                                                                                       cGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt             >  1:3147542/1‑100 (MQ=39)
                                                                                         gCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATAc           >  2:2370844/1‑100 (MQ=255)
                                                                                          cAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACt          >  2:3636047/1‑100 (MQ=255)
                                                                                             gATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCg       <  1:737494/100‑1 (MQ=255)
                                                                                               tattatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCg       <  1:1420321/98‑1 (MQ=255)
                                                                                               tattatCTCCGGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCggg     <  1:3311691/100‑1 (MQ=255)
                                                                                                attatCTCCTGCCGTTTATCTTTTTACTCGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGt    <  1:2341189/100‑1 (MQ=255)
                                                                                                attatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGt    >  1:1339561/1‑100 (MQ=255)
                                                                                                 ttatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTc   <  2:145749/100‑1 (MQ=255)
                                                                                                  tatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAACTCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCt  <  2:3458755/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGGTCT  >  NC_000913/2519127‑2519324

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: