Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 451,881 A→G intergenic (‑271/+189) cyoA ← / ← ampG cytochrome o ubiquinol oxidase subunit II/muropeptide transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913451,8810AG96.2% 587.3 / 7.0 185intergenic (‑271/+189)cyoA/ampGcytochrome o ubiquinol oxidase subunit II/muropeptide transporter
Reads supporting (aligned to +/- strand):  ref base A (4/3);  new base G (91/86);  total (96/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

Only 100 of 185 total aligned reads displayed.

ATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGACAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCG  >  NC_000913/451785‑451976
                                                                                                |                                                                                               
  aaCAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAAc                                                                                             <  1:2897662/99‑1 (MQ=255)
  aaCAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAAc                                                                                             >  1:180629/1‑99 (MQ=255)
   aCAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACt                                                                                            >  1:76670/1‑99 (MQ=255)
   aCAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACt                                                                                            >  2:1225208/1‑99 (MQ=255)
     aaaaGCGTAGTGCAGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAGTAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTgg                                                                                          <  2:3975454/99‑1 (MQ=255)
       aaGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGt                                                                                        <  1:2816629/99‑1 (MQ=255)
       aaGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGt                                                                                        <  1:1607339/99‑1 (MQ=255)
        aGGGGAGTGCTGAAGGAGTCTGTGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGtg                                                                                       <  2:875189/99‑1 (MQ=255)
          gggAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGtgat                                                                                     <  2:3711794/99‑1 (MQ=255)
          gggAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGtgat                                                                                     <  2:1329258/99‑1 (MQ=255)
           ggAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGtgatg                                                                                    <  1:754574/99‑1 (MQ=255)
              gTGCTGAAGGAGACTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATgaga                                                                                 >  1:2539605/1‑99 (MQ=255)
               tGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGaa                                                                                >  2:1119553/1‑99 (MQ=255)
                gCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGaaa                                                                               >  2:1690657/1‑99 (MQ=255)
                 cTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGaaaa                                                                              <  2:1036243/99‑1 (MQ=255)
                  tGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAAt                                                                             >  1:2859930/1‑99 (MQ=255)
                   gAAGGAGACGGGGCGGGCAATTGGTATAACCAATGCGAAATAAAACACTCATAAAACAAAGAATTATCATGAGAGTGCCAACTGGAAGATGAAAAAATa                                                                            >  2:839765/1‑99 (MQ=255)
                    aaGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATACTTATCATTAGAGTGGCAACTGGGTGATGTGCAAATaa                                                                           >  2:3502667/1‑99 (MQ=255)
                    aaGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATaa                                                                           >  2:3236198/1‑99 (MQ=255)
                     aGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATaaa                                                                          <  2:924635/99‑1 (MQ=255)
                     aGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATaaa                                                                          <  1:2147955/99‑1 (MQ=255)
                       gAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGa                                                                        <  1:1253802/99‑1 (MQ=255)
                        aGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGAc                                                                       >  2:4151763/1‑99 (MQ=255)
                         gTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTACCAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACt                                                                      <  2:4206634/99‑1 (MQ=255)
                         gTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACt                                                                      <  2:3188221/99‑1 (MQ=255)
                              ggcgggcAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTa                                                                 <  1:273765/99‑1 (MQ=255)
                              ggcgggcAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTa                                                                 >  1:1870527/1‑99 (MQ=255)
                               gcgggcAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAg                                                                >  2:260925/1‑99 (MQ=255)
                               gcgggcAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAg                                                                <  2:4031922/99‑1 (MQ=255)
                               gcgggcAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAg                                                                >  1:919234/1‑99 (MQ=255)
                                cgggcAATTGGTATAACCAATGTGAAATAGAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGc                                                               <  1:3969255/99‑1 (MQ=255)
                                 gggcAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCt                                                              <  2:2988547/99‑1 (MQ=255)
                                   gccaTTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGa                                                            >  1:3342247/4‑99 (MQ=255)
                                     aaTTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAAt                                                          >  2:391220/1‑99 (MQ=255)
                                     aaTTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAAt                                                          >  1:1173823/1‑99 (MQ=255)
                                     aaTTGGTATAACCAATGTGAAAGAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAAGAAAGACTAATTAGCTGAAt                                                          >  2:549693/1‑99 (MQ=255)
                                      aTTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATc                                                         >  2:2750681/1‑99 (MQ=255)
                                       ttGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCg                                                        >  2:1192797/1‑99 (MQ=255)
                                       ttGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGACAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAAtag                                                        >  1:307516/1‑97 (MQ=255)
                                        tGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCgt                                                       <  1:3908278/99‑1 (MQ=255)
                                          gTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCgtgt                                                     <  1:2438191/99‑1 (MQ=255)
                                          gTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCgtgt                                                     >  2:3218101/1‑99 (MQ=255)
                                          gTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCgtgt                                                     <  2:2737914/99‑1 (MQ=255)
                                          gTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGACAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCgtgt                                                     >  1:359934/1‑99 (MQ=255)
                                            ataACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTaa                                                   <  2:244517/99‑1 (MQ=255)
                                            ataACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTaa                                                   <  1:2333774/99‑1 (MQ=255)
                                            ataACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTaa                                                   <  1:4148674/99‑1 (MQ=255)
                                            ataACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTaa                                                   <  1:3436685/99‑1 (MQ=255)
                                             taACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGTAACGAGTAAg                                                  >  2:389761/1‑99 (MQ=255)
                                             taACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAg                                                  >  2:2113777/1‑99 (MQ=255)
                                             taACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAg                                                  >  2:2295253/1‑99 (MQ=255)
                                             taACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAg                                                  >  2:3197404/1‑99 (MQ=255)
                                              aaCCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGa                                                 >  2:474766/1‑99 (MQ=255)
                                              aaCCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGa                                                 >  1:2877108/1‑99 (MQ=255)
                                               aCCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGaa                                                >  2:2466519/1‑99 (MQ=255)
                                                ccAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGaaa                                               >  2:3023734/1‑99 (MQ=255)
                                                  aaTGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGaaaaa                                             <  1:2193222/99‑1 (MQ=255)
                                                  aaTGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGaaaaa                                             >  2:3386249/1‑99 (MQ=255)
                                                  aaTGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGaaaaa                                             <  1:248624/99‑1 (MQ=255)
                                                  aaTGTGAAATAAAACAATAATTAAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGaaaaa                                             <  2:3256144/99‑1 (MQ=255)
                                                    tgtgAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATa                                           <  2:1517821/99‑1 (MQ=255)
                                                     gtgAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAg                                          >  2:135391/1‑99 (MQ=255)
                                                      tgAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGa                                         >  2:3576845/1‑99 (MQ=255)
                                                      tgAAATAAAACAATTATTTAACAAATAATCATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGaaaaagatc                                         >  1:1410267/1‑95 (MQ=255)
                                                        aaaTAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAAt                                       >  2:3862951/1‑99 (MQ=255)
                                                          aTAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAAttt                                     >  1:356640/1‑99 (MQ=255)
                                                            aaaaCAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTtat                                   >  2:2805702/1‑99 (MQ=255)
                                                            aaaaCAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTtat                                   <  1:4094874/99‑1 (MQ=255)
                                                             aaaCAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTAAAAAAAATAAAAAAAATAGAATTtata                                  <  2:779344/99‑1 (MQ=255)
                                                                cAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAAt                               <  1:2607059/99‑1 (MQ=255)
                                                                 aaTTATTTAACAAATAATTATCATTAGAGTGGGAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATg                              <  2:4271973/99‑1 (MQ=255)
                                                                 aaTTATTTAACAAATAATTATCATTAGAGTGACAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATg                              <  2:354556/99‑1 (MQ=255)
                                                                 aaTGATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATGTATAAATg                              >  2:3623969/1‑99 (MQ=255)
                                                                     attTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGgcgc                          <  1:3259575/99‑1 (MQ=255)
                                                                     attTAACAAATAATTATCATTAGAGTGGCAACTGGGTGAGGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGgcgc                          <  1:2499120/99‑1 (MQ=255)
                                                                           cAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCaataa                    >  1:475184/1‑99 (MQ=255)
                                                                           cAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCaataa                    >  2:3795740/1‑99 (MQ=255)
                                                                            aaaTAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCaataat                   >  2:2919363/1‑99 (MQ=255)
                                                                            aaaTAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCaataat                   >  1:3063277/1‑99 (MQ=255)
                                                                            aaaTAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAAAAAAATAAAAAAAATAGAATTTATAAATGGCGCCaataat                   <  2:3071069/99‑1 (MQ=255)
                                                                              ataattataattAGAGTGGCAACTGGGTGAGGATAAAATAAAGACTAATTAGCTGAATCTTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAAttt                 <  2:3281729/99‑1 (MQ=255)
                                                                              ataatTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAAttt                 <  2:3527638/99‑1 (MQ=255)
                                                                              ataatTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAAttt                 <  2:2525195/99‑1 (MQ=255)
                                                                               taatTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTc                <  2:2502358/99‑1 (MQ=255)
                                                                               taatTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTc                <  1:962807/99‑1 (MQ=255)
                                                                               taatTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTc                <  2:3008404/99‑1 (MQ=255)
                                                                               taatTATCATTAGAGTGACAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTc                <  1:310130/99‑1 (MQ=255)
                                                                                aatTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTcc               <  1:61657/99‑1 (MQ=255)
                                                                                aatTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTcc               <  2:467518/99‑1 (MQ=255)
                                                                                   tATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCaaa            <  1:2978496/99‑1 (MQ=255)
                                                                                    aTCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAAc           >  1:1104533/1‑99 (MQ=255)
                                                                                      cATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGc         >  1:1581847/1‑99 (MQ=255)
                                                                                      cATTAGAGTGACAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGc         <  2:3092085/99‑1 (MQ=255)
                                                                                        ttAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTa       >  2:3138352/1‑99 (MQ=255)
                                                                                        tgagagTTGCAACTGGTTGATGAGTAAATAAAGACTCATTATCTGAATCGGGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTa       >  2:198409/3‑99 (MQ=255)
                                                                                          agagTGGCAACTGGGTGATGAGAAAATAAAGACTCATTGGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAg     <  1:1378437/99‑1 (MQ=255)
                                                                                           gagTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGc    >  2:2364323/1‑99 (MQ=255)
                                                                                            agTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGcc   <  2:2948429/99‑1 (MQ=255)
                                                                                            agTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGcc   <  1:744179/99‑1 (MQ=255)
                                                                                             gTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCg  >  2:1779089/1‑99 (MQ=255)
                                                                                                |                                                                                               
ATAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGACAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCG  >  NC_000913/451785‑451976

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: