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.5% 754.6 / 2.6 229intergenic (‑271/+189)cyoA/ampGcytochrome o ubiquinol oxidase subunit II/muropeptide transporter
Reads supporting (aligned to +/- strand):  ref base A (5/3);  new base G (120/98);  total (127/102)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 229 total aligned reads displayed.

TAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGACAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCG  >  NC_000913/451786‑451976
                                                                                               |                                                                                               
taACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCaa                                                                                              <  1:343893/99‑1 (MQ=255)
  aCAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACt                                                                                            >  1:1507270/1‑99 (MQ=255)
  aCAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACt                                                                                            >  1:218473/1‑99 (MQ=255)
  aCAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAACGTGAAATAAAACAATTATTTAACAAATAATTATCATTATAGGGGCAAca                                                                                            >  2:264714/1‑98 (MQ=255)
   cAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTg                                                                                           >  1:2916835/1‑99 (MQ=255)
    atctggggAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTgg                                                                                          <  1:5386289/95‑1 (MQ=255)
    aaaaGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTgg                                                                                          <  2:1582476/99‑1 (MQ=255)
     aaaGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTggg                                                                                         <  2:4611201/99‑1 (MQ=255)
      aaGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGTAACTGGGt                                                                                        <  1:3957958/99‑1 (MQ=255)
      aaGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGt                                                                                        <  2:800530/99‑1 (MQ=255)
        ggggAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGtga                                                                                      <  2:2985533/99‑1 (MQ=255)
         gggAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGtgat                                                                                     <  2:852170/99‑1 (MQ=255)
         gggAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGagat                                                                                     >  2:5269415/1‑99 (MQ=255)
          ggAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGtgatg                                                                                    >  1:2656829/1‑99 (MQ=255)
           gAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGtgatga                                                                                   >  1:3335028/1‑99 (MQ=255)
           gAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGtgatga                                                                                   >  2:3222777/1‑99 (MQ=255)
             gTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATgaga                                                                                 >  1:193979/1‑99 (MQ=255)
              tGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGaa                                                                                >  1:2866915/1‑99 (MQ=255)
               gCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGaaa                                                                               >  1:4545049/1‑99 (MQ=255)
               gCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGaaa                                                                               <  1:3848769/99‑1 (MQ=255)
                cTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGaaaa                                                                              >  1:3472740/1‑99 (MQ=255)
                cTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGACAACTGGGTGATGAGaaaa                                                                              >  1:4257130/1‑99 (MQ=255)
                 tGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAAt                                                                             <  1:3190209/99‑1 (MQ=255)
                  gAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATa                                                                            <  2:1754484/99‑1 (MQ=255)
                  gAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATa                                                                            <  1:1727463/99‑1 (MQ=255)
                   aaGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATaa                                                                           >  2:2444013/1‑99 (MQ=255)
                     ggAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAg                                                                         >  2:3101827/1‑99 (MQ=255)
                     ggAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAg                                                                         >  2:4219646/1‑99 (MQ=255)
                      gAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGa                                                                        <  2:2519372/99‑1 (MQ=255)
                      gAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGa                                                                        <  2:5184545/99‑1 (MQ=255)
                        gTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGACAACTGGGTGATGAGAAAATAAAGACt                                                                      <  1:3306234/99‑1 (MQ=255)
                         tCTGGGCGGGCAATTGGTATAACCAATGTGAAATAGAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACta                                                                     >  2:3029412/1‑98 (MQ=255)
                           tGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATGTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAAGAAAGACTCAt                                                                   >  2:7058/1‑99 (MQ=255)
                             ggcgggcAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTa                                                                 >  1:4501898/1‑99 (MQ=255)
                                gggcAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCt                                                              <  2:4330864/99‑1 (MQ=255)
                                  gcAATTTGTATAACCAATGTGAAATATAACAATTATTTAGCAAATTATTATCATTAGAGTGGCAACTGGGAGATGTGAAAATAAAGACTCATTAGCTGa                                                            <  1:5265555/99‑1 (MQ=255)
                                  gcAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGa                                                            >  1:4134014/1‑99 (MQ=255)
                                  gcAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTAGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGa                                                            >  2:112291/1‑99 (MQ=255)
                                    aaTTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAAt                                                          >  2:2750529/1‑99 (MQ=255)
                                     aTTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATc                                                         <  1:4662375/99‑1 (MQ=255)
                                     aTTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATc                                                         <  2:5496949/99‑1 (MQ=255)
                                     aTTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATc                                                         >  2:509293/1‑99 (MQ=255)
                                      ttGGTATAACCAATGTGAAATAAAACAATTATTCTACAAATAATTATCATTAGAGTGTCAACGGAGCAATGAGAAAATAAAGACTCATTAGCTGAATCg                                                        >  2:5413651/1‑99 (MQ=255)
                                        ggTATACCAATGCGAAATAAAACAATCATTCAACAAAGCATTATCATTAGAGTGGCAACTGGGAGATGAGACAATAAAGAATCATTAGCTGAATCttgt                                                     >  2:147878/1‑99 (MQ=255)
                                         gTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCgtgt                                                     <  1:1773914/99‑1 (MQ=255)
                                         gTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCgtgt                                                     <  1:1578007/99‑1 (MQ=255)
                                         gTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCgtgt                                                     >  2:505693/1‑99 (MQ=255)
                                         gTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCgtgt                                                     <  1:1022340/99‑1 (MQ=255)
                                         gTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTCGCAACTGGGTGATGAGAACATAAAGACTCCTTAGCTGAATCgtgt                                                     >  2:5507963/1‑99 (MQ=255)
                                         ctataACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCACCTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCgtgt                                                     <  2:3325114/98‑1 (MQ=255)
                                          tataACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTa                                                    >  1:704545/1‑99 (MQ=255)
                                           ataACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTaa                                                   <  1:3598934/99‑1 (MQ=255)
                                            taACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAg                                                  >  2:3867941/1‑99 (MQ=255)
                                            taACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAg                                                  <  1:1871877/99‑1 (MQ=255)
                                             aaCCAAGGTGAAATAGAACAATTATTTAACAAATAATTATCAGTAGAGTGGCAACTGGTTGATGAGAAAATAAAGACGCATTAGCTGAATCGAGTAAGa                                                 >  2:1891246/1‑99 (MQ=255)
                                               ccAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGaaa                                               >  1:2199080/1‑99 (MQ=255)
                                                 aaTGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGaaaaa                                             <  1:3630526/99‑1 (MQ=255)
                                                  aTGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAAt                                            <  2:1292467/99‑1 (MQ=255)
                                                  aTGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAAt                                            >  2:4762763/1‑99 (MQ=255)
                                                      gAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGaa                                        <  1:3325428/99‑1 (MQ=255)
                                                      gAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGaa                                        >  1:619414/1‑99 (MQ=255)
                                                      gAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGaa                                        >  1:4740149/1‑99 (MQ=255)
                                                       aaaTAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAAt                                       <  1:3209492/99‑1 (MQ=255)
                                                       aaaTAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAAt                                       <  2:170151/99‑1 (MQ=255)
                                                         aTAAAACAATTATTTAACACATAATTATCATTAGAGTGGCAACTGGATGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAAttt                                     >  2:2292775/1‑99 (MQ=255)
                                                         aTAAAACAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAAttt                                     >  1:1070213/1‑99 (MQ=255)
                                                           aaaaCAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTtat                                   >  2:2340818/1‑99 (MQ=255)
                                                              aCAATTATTTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATaaa                                >  1:2006550/1‑99 (MQ=255)
                                                              aCAATTATTAAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATaaa                                <  1:2253361/99‑1 (MQ=255)
                                                                  ttattTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGgc                            >  2:3155408/1‑99 (MQ=255)
                                                                  ttattTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGgc                            >  2:4208309/1‑99 (MQ=255)
                                                                    attTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGgcgc                          >  1:11757/1‑99 (MQ=255)
                                                                    attTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGgcgc                          <  1:1250102/99‑1 (MQ=255)
                                                                    attTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGgcgc                          >  2:1073355/1‑99 (MQ=255)
                                                                    attTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGgcgc                          <  2:3498748/99‑1 (MQ=255)
                                                                    attTAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGgcgc                          <  2:4030146/99‑1 (MQ=255)
                                                                     tttAACAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAATGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATGTATAAATGGCGcc                         >  1:5088438/1‑99 (MQ=255)
                                                                        aaCAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCaat                      >  1:2298765/1‑99 (MQ=255)
                                                                         aCAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCaata                     >  2:1449434/1‑99 (MQ=255)
                                                                         aCAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCaata                     >  2:1197744/1‑99 (MQ=255)
                                                                         aCAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCaata                     >  2:343324/1‑99 (MQ=255)
                                                                          cAAATAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCaataa                    >  1:488350/1‑99 (MQ=255)
                                                                           aaaTAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCaataat                   >  1:483101/1‑99 (MQ=255)
                                                                           aaaTAATTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCaataat                   >  2:50618/1‑99 (MQ=255)
                                                                            aataatTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAAtt                  >  1:2750841/1‑99 (MQ=255)
                                                                            aataatTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGAATCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAAtt                  <  2:4969824/99‑1 (MQ=255)
                                                                             ataatTATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAAttt                 <  1:5639214/99‑1 (MQ=255)
                                                                                 ttATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCaa             <  1:553953/99‑1 (MQ=255)
                                                                                  tATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCaaa            <  1:1339822/99‑1 (MQ=255)
                                                                                  tATCATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCAGTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCaaa            <  2:566099/99‑1 (MQ=255)
                                                                                  tATCATTAGAGTGGAAAAAGGGTAATGAGAAAATAAAGACTCATTAGCTGAAAAAAAAAAAAAAAATAGAATTTATAAATGGCGCCAATAATTTCCaaa            <  2:1343067/99‑1 (MQ=255)
                                                                                  tATCATTAGAGTGACAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCaaa            >  1:2518804/1‑99 (MQ=255)
                                                                                     gaTTAGAGTGACAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGc         >  2:5113047/2‑99 (MQ=255)
                                                                                     cATTAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGc         >  2:379671/1‑99 (MQ=255)
                                                                                       ttAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTa       >  2:4305612/1‑99 (MQ=255)
                                                                                        tAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTATCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTaa      <  1:2683564/99‑1 (MQ=255)
                                                                                        tAGAGTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTaa      >  2:2159716/1‑99 (MQ=255)
                                                                                          gagTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGc    >  1:1533607/1‑99 (MQ=255)
                                                                                          gagTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGc    >  2:3654461/1‑99 (MQ=255)
                                                                                            gTGGCAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCg  >  2:308358/1‑99 (MQ=255)
                                                                                               |                                                                                               
TAACAAAAGGGGAGTGCTGAAGGAGTCTGGGCGGGCAATTGGTATAACCAATGTGAAATAAAACAATTATTTAACAAATAATTATCATTAGAGTGACAACTGGGTGATGAGAAAATAAAGACTCATTAGCTGAATCGTGTAAGAAAAATAGAATTTATAAATGGCGCCAATAATTTCCAAACGCTAAGCCG  >  NC_000913/451786‑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: