Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 801,498 T→A 7.6% V247D (GTT→GAT)  ybhH → putative PrpF family isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913801,4980TA 7.6% 183.8 / 2.4 60V247D (GTT→GAT) ybhHputative PrpF family isomerase
Reads supporting (aligned to +/- strand):  ref base T (30/25);  new base A (2/3);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTGGAAGGCACCGTCAC  >  NC_000913/801352‑801645
                                                                                                                                                  |                                                                                                                                                   
cgtcATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGtt                                                                                                                                                    <  2:45688/148‑1 (MQ=255)
 gtcATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTAt                                                                                                                                                  <  3:83650/149‑1 (MQ=255)
 gtcATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTAt                                                                                                                                                  >  7:65770/1‑149 (MQ=255)
 gtcATTATTCCGGCTGAATATCT‑GGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATc                                                                                                                                                 >  6:358076/1‑149 (MQ=255)
    attattCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGt                                                                                                                                                     <  4:284355/143‑1 (MQ=255)
    attattCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGt                                                                                                                                                     >  3:284355/1‑143 (MQ=255)
     ttattCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCt                                                                                                                                              <  5:67865/149‑1 (MQ=255)
                gAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAAc                                                                                                                                          >  4:1996/1‑142 (MQ=255)
                gAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAAc                                                                                                                                          <  3:1996/142‑1 (MQ=255)
                 aaTATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCt                                                                                                                                   >  7:104044/1‑148 (MQ=255)
                 aaTATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGACAGTAATATGGTTATCCCTAAACCTGTGCt                                                                                                                                   >  6:162395/1‑148 (MQ=255)
                 aaTATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATACGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCt                                                                                                                                   >  2:312003/1‑148 (MQ=255)
                  atatCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCtt                                                                                                                                  <  8:104044/148‑1 (MQ=255)
                       tGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTtct                                                                                                                            >  5:279970/1‑149 (MQ=255)
                       tGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTtct                                                                                                                            >  8:90315/1‑149 (MQ=255)
                          gAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAATCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCa                                                                                                                         <  6:279970/149‑1 (MQ=255)
                            aaaaCAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAgc                                                                                                                       <  1:312003/149‑1 (MQ=255)
                                 aGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACACGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAg                                                                                                                   >  5:317451/1‑148 (MQ=255)
                                          aTTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCggg                                                                                                          <  5:358076/148‑1 (MQ=255)
                                               cggcggAACAGGATGCCGACAAAGCATTATTAGCCCCCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAAt                                                                                                     >  8:138573/1‑148 (MQ=255)
                                                 gcggAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTaa                                                                                                  <  8:114663/149‑1 (MQ=255)
                                                               cGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTa                                                                                    <  3:73410/149‑1 (MQ=255)
                                                               cGACAAACGATTATTAGCCCGCATTGACTCTATCCGTCTACAAGCGGGTAAAGCAATGGTCTTAGGTGATGTCAGTAATCTGGATATCCCTAACCCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTa                                                                                    <  6:333168/149‑1 (MQ=255)
                                                                   aaaGCATTATTAGCCCGCATTGAATCTATACGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGcc                                                                                >  4:19537/1‑149 (MQ=255)
                                                                   aaaGCATTATTACACAGCATTGAATCTATCCGTCTACAATCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGATATCCCTAAAACTGTGCATATTTCTCCATCGCAGAAACGCGGGGCAATTAATGTCCGTTATTTTATccc                                                                                >  4:42788/1‑149 (MQ=255)
                                                                    aaGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCg                                                                               >  3:150878/1‑149 (MQ=255)
                                                                        attattAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGAGCGTTATTTTATGCCGCAt                                                                            >  1:190182/1‑148 (MQ=255)
                                                                            ttaGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCtt                                                                        <  3:221631/148‑1 (MQ=255)
                                                                               gCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCa                                                                    >  3:319205/1‑149 (MQ=255)
                                                                                 ccacATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATc                                                                  >  5:126399/4‑149 (MQ=255)
                                                                                  cGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGATATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGc                                                                                             >  4:306063/1‑121 (MQ=255)
                                                                                  cGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGATATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGc                                                                                             <  3:306063/121‑1 (MQ=255)
                                                                                     aTTGAATCTATCCGTCTACAAGCGTGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATcgcgc                                                              <  6:370837/149‑1 (MQ=255)
                                                                                          aTCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTgg                                                         <  8:65770/149‑1 (MQ=255)
                                                                                           tCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGc                                                        >  1:184011/1‑149 (MQ=255)
                                                                                              aTCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGa                                                      >  4:193781/1‑148 (MQ=255)
                                                                                                  gTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCTTTCTTGCCATCGCGCTCTGGCGATAAcc                                                 >  4:95466/1‑149 (MQ=255)
                                                                                                  gTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTACTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATATTATGCCGCATTCTTCCCATCTCGCGCTGGCGATAAcc                                                 >  7:383671/1‑149 (MQ=255)
                                                                                                          gCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGtgctat                                         <  2:190182/149‑1 (MQ=255)
                                                                                                           cGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGTGCAATTAATGTGCGTTATTTTATGCCGCATTGTTGCCATCGCGCGCTGGCGATAACCGGtgctatt                                        <  6:317451/149‑1 (MQ=255)
                                                                                                               tAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGtgctattgcta                                    >  5:256219/1‑149 (MQ=255)
                                                                                                                 aaGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTAt                                                                                   >  1:189685/1‑100 (MQ=255)
                                                                                                                 aaGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTAt                                                                                   <  2:189685/100‑1 (MQ=255)
                                                                                                                  aGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTAttt                                 >  2:223445/1‑149 (MQ=255)
                                                                                                                   ggaaTGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTc                                <  3:19537/147‑1 (MQ=255)
                                                                                                                     aaTGGGCTTAGGTGATGTCAGTTATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTTATGTGCGTTTTTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCa                              <  3:42788/149‑1 (MQ=255)
                                                                                                                     aaTGGGCTTAGGTGATGTCAGTAATATGGTTATCACTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCTCTGGCGATAACCGGTGCTATTGCTATTTCCa                              <  1:223445/149‑1 (MQ=255)
                                                                                                                      aTGGGCTTAGTTTTGTTCAGTATTATGGTTATCCCTAACCCTGTGCTTATTTCTCAAGCGCAGAATGGCGGGGAAATAATTGTGC‑TTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGtgctattg                                       <  6:25773/138‑1 (MQ=255)
                                                                                                                      aTGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCag                             >  8:218368/1‑149 (MQ=255)
                                                                                                                        gggCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCag                             >  2:68827/1‑147 (MQ=255)
                                                                                                                         ggCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCagta                           <  7:218368/148‑1 (MQ=255)
                                                                                                                              aGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCgcgcgc                                                            <  8:369643/110‑1 (MQ=255)
                                                                                                                              aGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCgcgcgc                                                            >  7:369643/1‑110 (MQ=255)
                                                                                                                              aGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTtgt                      <  7:138573/148‑1 (MQ=255)
                                                                                                                              aGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTtgt                      <  7:90315/148‑1 (MQ=255)
                                                                                                                                  gATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCAtt                 >  8:98220/1‑149 (MQ=255)
                                                                                                                                  gATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCAtt                 >  7:22099/1‑149 (MQ=255)
                                                                                                                                    tGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTgg               <  5:162395/149‑1 (MQ=255)
                                                                                                                                      tCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGAGTGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTgg               >  6:366643/1‑147 (MQ=255)
                                                                                                                                        aGTAAAAAGGATATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATAAATCTGCCTAATATTAAGCCCCACTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTg                <  8:383671/144‑1 (MQ=255)
                                                                                                                                         gTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTGGAAGGc          <  4:319205/149‑1 (MQ=255)
                                                                                                                                                 gTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTGGAAGGCACCGTCAc  >  4:277591/1‑149 (MQ=255)
                                                                                                                                                 gTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTGGAAGGCACCGTCAc  >  1:32366/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGTCATTATTCCGGCTGAATATCTGGGAAAAACAGGTTATGAATTACCGGCGGAACTGGATGCCGACAAAGCATTATTAGCCCGCATTGAATCTATCCGTCTACAAGCGGGTAAAGCAATGGGCTTAGGTGATGTCAGTAATATGGTTATCCCTAAACCTGTGCTTATTTCTCCAGCGCAGAAAGGCGGGGCAATTAATGTGCGTTATTTTATGCCGCATTCTTGCCATCGCGCGCTGGCGATAACCGGTGCTATTGCTATTTCCAGTAGTTGTGCATTGGAAGGCACCGTCAC  >  NC_000913/801352‑801645

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑