Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,519,222 T→G 17.7% 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,2220TG17.7% 286.8 / 61.2 125intergenic (+17/+35)yfeD/gltXDUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (47/55);  new base G (12/10);  total (59/66)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.95e-01

Only 100 of 131 total aligned reads displayed.

AGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGG  >  NC_000913/2519123‑2519321
                                                                                                   |                                                                                                   
aGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGat                                                                                                     >  1:3171776/1‑100 (MQ=255)
 gCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGata                                                                                                    <  1:2453301/100‑1 (MQ=255)
    ggAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatta                                                                                                 >  1:1682982/1‑100 (MQ=255)
       aaaaGTCGCGCAAGTTTAGTCACGTGAACGAAATCGAAAAAAAATCGCCCGTCTCGATATTGACGAAAAAAAAAAAAAAAAAAACGGCAGGATATTAtc                                                                                               <  2:1698003/99‑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)
                 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)
                      tttGTCACGTGAAGGATTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttttat                                                                               <  2:1004087/100‑6 (MQ=37)
                      tttGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttttat                                                                               <  1:778858/100‑6 (MQ=37)
                         gTCACGTGAAGGGATTAGCAATTGAAGAGCGGGTCTCGATATTGACGAATCAGCATGAATAAAAAACGGCAGGATGTTATCTCCTGCCGTTTATc                                                                                 <  2:3015993/95‑1 (MQ=255)
                         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)
                            aCGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTaca                                                                         >  1:145749/1‑100 (MQ=255)
                             cGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTacac                                                                        >  2:1813077/1‑100 (MQ=255)
                             cGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTacac                                                                        <  1:1705687/100‑1 (MQ=255)
                               tGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGc                                                                      >  1:2054937/1‑100 (MQ=255)
                               tGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGc                                                                      >  1:1348429/1‑100 (MQ=255)
                               tGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGc                                                                      >  2:185538/1‑100 (MQ=255)
                                gAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCt                                                                     >  1:241481/1‑100 (MQ=255)
                                 aaGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTa                                                                    <  1:3470680/100‑1 (MQ=255)
                                 aaGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTa                                                                    <  1:2807101/100‑1 (MQ=255)
                                    gAATTCGCAATTTGCTCGCCCGTCGCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAAt                                                                  <  2:1580023/99‑1 (MQ=255)
                                     aaTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTa                                                                <  1:2884303/100‑1 (MQ=255)
                                      aTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAc                                                               >  2:3487237/1‑100 (MQ=255)
                                      aTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAc                                                               >  1:1283602/1‑100 (MQ=255)
                                       ttCGCAATTTCCTCACCCGTCTCCATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAc                                                               <  2:3409512/99‑1 (MQ=255)
                                        tCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActg                                                             <  2:1897289/100‑1 (MQ=255)
                                        tCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActg                                                             <  2:2589522/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:737491/1‑100 (MQ=255)
                                            aaTTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGa                                                         <  1:3326901/100‑1 (MQ=255)
                                               ttCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGAttt                                                      <  2:3286004/100‑1 (MQ=255)
                                               ttCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTCACACGCTAATTACTGCTGAttt                                                      <  2:1050070/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)
                                                 ccTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTc                                                    >  1:1341388/1‑100 (MQ=255)
                                                   tagCCCGTCTAGATATTGACGAATCAACATAAATAAATAACTGCAGGATATTATCTCCTGCCGTTTATCTTTTTACTCGCTAATTACTGCTGATTTTcg                                                   <  2:2973805/97‑1 (MQ=255)
                                                   tCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgc                                                  <  1:965590/100‑1 (MQ=255)
                                                    cGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgcg                                                 <  2:2002873/100‑1 (MQ=255)
                                                    cGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgcg                                                 >  2:66170/1‑100 (MQ=255)
                                                    cGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCGTTTTACACGCTAATTACTGCTGATTTTcg                                                   <  2:3341363/98‑1 (MQ=255)
                                                      cccGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGtt                                               >  1:1071259/1‑100 (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:3175777/100‑1 (MQ=255)
                                                                tatTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaa                                     >  1:2160021/1‑100 (MQ=255)
                                                                  ttGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAAt                                   >  2:2051511/1‑100 (MQ=255)
                                                                  ttGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAAt                                   >  1:495027/1‑100 (MQ=255)
                                                                    gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                 <  1:3360967/100‑1 (MQ=255)
                                                                    gACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATcc                                 <  1:2181033/100‑1 (MQ=255)
                                                                     aCGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCa                                <  1:2656469/100‑1 (MQ=255)
                                                                       gAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgc                              <  1:1217944/100‑1 (MQ=255)
                                                                          tCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTATTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgcgc                            <  2:2324208/99‑1 (MQ=255)
                                                                            aGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCttt                         <  1:2575774/100‑1 (MQ=255)
                                                                            aGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCtt                          >  2:1191065/1‑99 (MQ=255)
                                                                             gCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttg                        >  1:288885/1‑100 (MQ=255)
                                                                              cATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgt                       >  1:536536/1‑100 (MQ=255)
                                                                              cATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgt                       >  2:2769026/1‑100 (MQ=255)
                                                                               ataaataaAAAACTGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                      <  2:403202/100‑1 (MQ=255)
                                                                                taaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgttg                     >  2:2491247/1‑100 (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)
                                                                                    taaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATAc                 >  1:2683959/1‑100 (MQ=255)
                                                                                      ataaaCGGCAGGAGATAATATCCTGCCGTTTATCTCTTTACACGCTAATTACTGTTGATTTTCGCGTTCTGCAATCAAATCCAGCGCTTTGTTGCTACGc               >  2:3363155/3‑100 (MQ=25)
                                                                                      aaaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGc               >  2:1457961/1‑100 (MQ=255)
                                                                                      aaaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGAAATAAAATCCAGCGCTTTGTTGATACGc               >  2:1412738/1‑100 (MQ=255)
                                                                                       aaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACgat              >  2:2203759/1‑98 (MQ=255)
                                                                                        aaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTc             >  1:1873008/1‑100 (MQ=39)
                                                                                         aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg            >  1:3184177/1‑100 (MQ=39)
                                                                                         aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg            >  1:1856923/1‑100 (MQ=39)
                                                                                         aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg            >  1:1544748/1‑100 (MQ=39)
                                                                                         aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg            >  2:3490642/1‑100 (MQ=39)
                                                                                          aCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATaagttc             >  2:1146005/1‑93 (MQ=38)
                                                                                           cGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt          >  1:3147542/1‑100 (MQ=39)
                                                                                           cGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGAt          >  2:1224605/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)
                                                                                                  atatTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCgg   <  1:3309106/100‑1 (MQ=255)
                                                                                                   tattatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCg    <  1:1420321/98‑1 (MQ=255)
                                                                                                   tattatCTCCGGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCggg  <  1:3311691/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGG  >  NC_000913/2519123‑2519321

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: