Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,519,222 T→G 15.2% 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,2220TG15.2% 312.2 / 51.5 130intergenic (+17/+35)yfeD/gltXDUF1323 family putative DNA‑binding protein/glutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (54/56);  new base G (15/5);  total (69/61)
Fisher's exact test for biased strand distribution p-value = 4.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

Only 100 of 144 total aligned reads displayed.

AGCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCGGG  >  NC_000913/2519123‑2519321
                                                                                                   |                                                                                                   
 gCAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGata                                                                                                    >  1:757576/1‑100 (MQ=255)
  cAGGAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGatat                                                                                                   <  1:448311/100‑1 (MQ=255)
    gtaaaaaGTCGCGCAATTTTTGTAACGTGAAGGAATTCGCAATTTCATCGTCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatta                                                                                                 >  2:549414/3‑100 (MQ=255)
    ggAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatta                                                                                                 >  2:2059284/1‑100 (MQ=255)
    ggAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtatta                                                                                                 >  2:617778/1‑100 (MQ=255)
     gAAAAAGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCTTCGTCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAtattat                                                                                                <  1:1024555/100‑1 (MQ=255)
      aaaaaGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtc                                                                                               >  2:391304/1‑100 (MQ=255)
      aaaaaGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtc                                                                                               <  1:571623/100‑1 (MQ=255)
      aaaaaGTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTAtc                                                                                               <  1:2827637/100‑1 (MQ=255)
           gTCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCt                                                                                           >  1:2924821/1‑99 (MQ=255)
            tCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGc                                                                                         <  2:1089354/100‑1 (MQ=255)
            tCGCGCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGc                                                                                         <  1:135125/100‑1 (MQ=255)
             cgcgCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCGGTCTCGATATTGAGGAAGCAGCATAAATAAAAAACGCCAGGATATTAGGTCCTga                                                                                         >  2:3931802/1‑98 (MQ=255)
             cgcgCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGcc                                                                                        >  2:2451925/1‑100 (MQ=255)
             cgcgCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGcc                                                                                        <  1:400111/100‑1 (MQ=255)
             cgcgCAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGc                                                                                         >  2:3660943/1‑99 (MQ=255)
              gcgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCg                                                                                       >  1:1672816/1‑100 (MQ=255)
               cgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGcc                                                                                        <  2:1398414/98‑1 (MQ=255)
               cgcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGt                                                                                      <  1:3532470/100‑1 (MQ=39)
               cacAATTTTTGTCACGTTAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGt                                                                                      <  1:2725710/98‑1 (MQ=38)
                gcAATTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATAAAAAAAAAAAAAAAACGGCAGGATATTATCTCCTGCCGtt                                                                                     <  2:2813860/100‑1 (MQ=255)
                   aTTTTTGTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTAt                                                                                  <  2:480163/100‑1 (MQ=255)
                         gTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttttt                                                                            <  2:1157858/100‑1 (MQ=255)
                         gTCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCttttt                                                                            <  1:834634/100‑1 (MQ=255)
                          tCACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTa                                                                           <  1:3460056/100‑1 (MQ=255)
                           cACGTGAAGGAATTCGCAATTTGCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttttatttatg                                                                          <  2:2257620/100‑11 (MQ=25)
                           cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTTGATATTTACGAATCAGCATAAATAAAAAACGGCAGGAGATAATATCCTGCCGttttttatttatg                                                                          <  1:3605721/100‑11 (MQ=21)
                           cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTac                                                                          >  2:2263867/1‑100 (MQ=255)
                           cACGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTac                                                                          <  1:3319408/100‑1 (MQ=255)
                             cGTGAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTacac                                                                        >  1:1034520/1‑100 (MQ=255)
                                gAAGGAATTCGCAATTTCCTCGCCCGTCTCGATATTGAAGAATCAGAATAAATAAAAAACGGTAGGATATTATCTCCTGCCGTTTATCTTTTTAAACGCt                                                                     >  2:179941/1‑100 (MQ=255)
                                 aaGGAATTCGCAATTACATCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAATCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTacac                                                                        <  2:867434/96‑1 (MQ=255)
                                     aaTTCGCAATTTCCTCGCCCGTCTCGATATTGACGCATCAGCATAAATAAACAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTa                                                                >  1:2871637/1‑100 (MQ=255)
                                      aTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAc                                                               <  1:1485165/100‑1 (MQ=255)
                                      aTTCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAc                                                               >  1:309134/1‑100 (MQ=255)
                                       ttCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATTTTTTTACACGCTAATTAct                                                              <  2:393886/100‑1 (MQ=255)
                                       ttCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTAct                                                              >  2:348354/1‑100 (MQ=255)
                                        tCGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActg                                                             >  2:262769/1‑100 (MQ=255)
                                         cGCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                            >  2:3840235/1‑100 (MQ=255)
                                          gCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgct                                                           <  1:201035/100‑1 (MQ=255)
                                          gCAATTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTActgc                                                            >  2:2725494/1‑99 (MQ=255)
                                            aaTTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGa                                                         <  2:2991453/100‑1 (MQ=255)
                                             aTTTCCTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCATGAAGATAATCTTTTTACACGCTAATTACTGCTGAt                                                        >  1:115537/1‑100 (MQ=255)
                                                 ccTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTc                                                    >  1:2307167/1‑100 (MQ=255)
                                                 ccTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTc                                                    >  1:979601/1‑100 (MQ=255)
                                                  cTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCTTGCCGTTTAACTTTTTACACGATAATTACAACTGAtt                                                       >  1:2863899/1‑96 (MQ=255)
                                                  cTCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcg                                                   <  2:2321643/100‑1 (MQ=255)
                                                   tCGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgc                                                  <  1:391310/100‑1 (MQ=255)
                                                    cGCCCGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTcgcg                                                 <  2:2099205/100‑1 (MQ=255)
                                                        cGTCTCGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCa                                             <  1:2331095/100‑1 (MQ=255)
                                                          tctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGc                                           <  1:1067656/100‑1 (MQ=255)
                                                          tctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGc                                           <  1:2009137/100‑1 (MQ=255)
                                                           ctcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCa                                          >  2:3683927/1‑100 (MQ=255)
                                                            tcGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCaa                                         >  1:2599277/1‑100 (MQ=255)
                                                             cGATATTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAAt                                        >  1:63370/1‑100 (MQ=255)
                                                               atatTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATTTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATa                                       <  1:3106616/99‑1 (MQ=255)
                                                                tatTGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaaa                                     >  2:149622/1‑100 (MQ=255)
                                                                tatTGAAGAATCAGTATAAATAAAAATCGGTAGTTTATTATTTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATaa                                      <  1:360911/99‑1 (MQ=255)
                                                                  ttGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAAt                                   >  1:2400805/1‑100 (MQ=255)
                                                                   tGACGAATCAGCATAAATAAAAAACTGCTGTATTTTATTTCCTGCCGTTTATCTTTTAACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATc                                  <  1:1247773/100‑1 (MQ=255)
                                                                   tGACGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATc                                  <  1:464580/100‑1 (MQ=255)
                                                                      cGAATCAGCATAAATAAAAAACGGCAGGATATTATTTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAg                               <  1:578680/100‑1 (MQ=255)
                                                                      cGAATCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAg                               <  1:937867/100‑1 (MQ=255)
                                                                        aaTCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgcg                             >  1:3549047/1‑100 (MQ=255)
                                                                        aaTCAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAgcg                             >  1:3335064/1‑100 (MQ=255)
                                                                           cAGCATAAATAAAAAACTGCAGGATATTATCTCATGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTGAGCAATAAAATCCAGCGCt                           >  2:437766/1‑99 (MQ=255)
                                                                           cAGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCtt                          <  1:3350183/100‑1 (MQ=255)
                                                                            aGCATAAATAAAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCttt                         >  2:3434416/1‑100 (MQ=255)
                                                                               tctaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                      <  2:388259/97‑1 (MQ=255)
                                                                               taaaaagataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                      >  2:3452223/10‑100 (MQ=25)
                                                                               ataaataaAAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgtt                      <  2:2939432/100‑1 (MQ=255)
                                                                                taaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATACAGCGCTttgttg                     >  1:3452547/1‑100 (MQ=255)
                                                                                taaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCACGTTCAGAAATAAAATCCAGCGCTttgtt                      >  1:1348356/1‑99 (MQ=255)
                                                                                taaataaaAAACGGCAGGATATTATCTACTGCCGTTAATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATACAGCGCTttgtt                      >  2:3519319/1‑99 (MQ=255)
                                                                                 aaataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGa                    <  1:1196268/100‑1 (MQ=255)
                                                                                  aataaaAAACGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTttgttaga                   >  1:388265/1‑97 (MQ=255)
                                                                                  aataaaAAACGGAAGAATATTATCTACTGCCGTTTATCTTTTTACACGATAATTACTTCTGATTTTCGGGTTTAACAATAAAATCCAACGCTTTGTTGAt                   >  2:1427780/1‑100 (MQ=255)
                                                                                  aaagataaaCGGCAGGAGATAATTTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAt                   >  2:2238248/7‑100 (MQ=37)
                                                                                  aaagataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTTATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGAt                   >  2:1011373/7‑100 (MQ=37)
                                                                                   ataaaAAACGGCAGGATATTATCTCCTGCTGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                  <  1:1525091/100‑1 (MQ=255)
                                                                                   agaaaaaaCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                  >  1:2653019/3‑100 (MQ=255)
                                                                                   aagataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATa                  >  1:318941/6‑100 (MQ=37)
                                                                                      ataaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGc               >  1:3189229/3‑100 (MQ=38)
                                                                                       taaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCt              >  2:1368122/2‑100 (MQ=38)
                                                                                        aaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTc             >  2:2624535/1‑100 (MQ=39)
                                                                                        aaaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTc             >  1:2541151/1‑100 (MQ=38)
                                                                                        aaaCAGCAGGTTATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTGCAGCAATAAAATCCAGCGCTTTGTTGATACGCTc             <  2:2658844/100‑1 (MQ=255)
                                                                                         aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg            >  1:817800/1‑100 (MQ=39)
                                                                                         aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCg            >  2:3247384/1‑100 (MQ=39)
                                                                                         aaCGGCAGGAGATAATATCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTAGCAATAAGAAATAAAAACCAGAGATTTGTTGATACGCTCg            >  1:3129574/1‑100 (MQ=25)
                                                                                         aaCGGCAGGAGACCATATCCTGCAGCTTATCTTTTTACACGATAATTAATGCTGATGTTCGCGTTCAGCAATAGAATGACGCGCTTTGTTGATACACTc             >  2:163420/1‑99 (MQ=255)
                                                                                          aCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGa           <  1:2969130/100‑1 (MQ=255)
                                                                                          aCGGCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGa           >  1:2524684/1‑100 (MQ=255)
                                                                                             gCAGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCACGATCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATAc        >  1:1039044/1‑100 (MQ=255)
                                                                                              cAGGATATTATCTCCTGCCGTTTATCTTTTTATTTTTTTTTTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAACGCTTTGTTGATACGCTCGATACt       >  2:246888/1‑100 (MQ=255)
                                                                                               aGGATATTATCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTg      >  2:898652/1‑100 (MQ=255)
                                                                                                ggATATTATCTCCTGCCGTTTATCTTTTTACACTCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGc     >  2:3397704/1‑100 (MQ=255)
                                                                                                 aatatTATCTCCTGCCGTTTATCTTATTACACACTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTTATACGCTCGATACTGCg    >  2:230946/2‑100 (MQ=255)
                                                                                                   tattatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGCAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCggg  <  2:309696/100‑1 (MQ=255)
                                                                                                   tattatCTCCTGCCGTTTATCTTTTTACACGCTAATTACTGCTGATTTTCGCGTTCAGAAATAAAATCCAGCGCTTTGTTGATACGCTCGATACTGCggg  <  2:1363741/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 < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: