Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,794,412 G→A 6.0% A54T (GCC→ACC)  gabP → gamma‑aminobutyrate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,794,4120GA 6.0% 207.9 / 2.6 66A54T (GCC→ACC) gabPgamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base G (28/34);  new base A (2/2);  total (30/36)
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.55e-01

TCGCAACCACATGAGTTAGGCGGCGGGCTGAAGTCACGCCACGTCACCATGTTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCTGGTTGTACTGGTGGTTTTGGG  >  NC_000913/2794265‑2794559
                                                                                                                                                   |                                                                                                                                                   
tcGCAACCACATGAGTTAGGCGGCGGGCTGAAGTCACGCCACGTCACCATGTTGTCTATTTCCGGTGTTATCGGCGCAAGTCTGTTTGGCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATATGTTCGc                                                                                                                                                    >  1:135436/1‑149 (MQ=255)
    aaCCACATGAGTTAGGCGGCGGGCTGAAGTCACGCCACGTCACCATGTTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCggg                                                                                                                                                <  6:81616/149‑1 (MQ=255)
      ccACATGAGTTAGGCGGCGGGCTGAAGTCACGCCACGTCACCATGTTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGc                                                                                                                                               >  2:95224/1‑148 (MQ=255)
           tGAGTTAGGCGGCGGGCTGAAGTCACGCCACGTCACCATGTTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACtgg                                                                                                                                         <  8:43166/149‑1 (MQ=255)
              gTTAGGCGGCGGGCTGAAGTCACGCCACGTCACCATGTTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACtggtgg                                                                                                                                      <  8:83494/149‑1 (MQ=255)
              gTTAGGCGGCGGGCTGAAGTCACGCCACGTCACCATGTTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACtggtgg                                                                                                                                      >  8:16532/1‑149 (MQ=255)
               ttAGGCGGCGGGCTGAAGTCACGCCACGTCACCATGTTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTggtggt                                                                                                                                     >  4:78452/1‑149 (MQ=255)
                  ggcggcGGGCTGAAGTCACGCCACGTCACCATGTTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTa                                                                                                                                   <  1:136324/148‑1 (MQ=255)
                  ggcggcGGGCTGAAGTCACGCCACGTCACCATGTTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTAt                                                                                                                                  <  3:78452/149‑1 (MQ=255)
                    cggcgtGCTGAAGTCACGCCACGTCACCATGGTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGa                                                                                                                                <  8:308551/149‑1 (MQ=255)
                      gcggGCTGAAGTCACGCCACGTCACCCTGTTGTCTATTGCCGGTGTTATCGGCGCTAGTCTGTTTGTCGGGTCCAGCGGCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGAtt                                                                                                                              <  2:299096/149‑1 (MQ=255)
                           cTGAAGTCACGCCACGTCACCATGTTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCg                                                                                                                         >  4:280825/1‑149 (MQ=255)
                           cTGAAGTCACGCCACGTCACCATGTTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGGATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCg                                                                                                                         >  4:163029/1‑149 (MQ=255)
                            tGAAGTCACGCCACGTCACCATGTTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCg                                                                                                                         >  2:283151/1‑148 (MQ=255)
                              aaGTCACGCCACGTCACCATGTTGTCTATTGCCTGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGa                                                                                                                       <  4:303350/148‑1 (MQ=255)
                                 tCACGCCACGTCACCATGTTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGtt                                                                                                                   <  7:16532/149‑1 (MQ=255)
                                   aCGCCACGTCACCATGTTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCACCGGGCTACTGGTGGTTATGATTATGCGGATGTTgg                                                                                                                 >  4:114803/1‑149 (MQ=255)
                                   aCGCCACGTCACCATGTTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCACCGGGCTACTGGTGGTTATGATTATGCGGATGTTgg                                                                                                                 >  4:110293/1‑149 (MQ=255)
                                    cGCCACGTCACCATGTTGTCTAATGACGGTGTTCTCGGCGCTAGTGTGTTTGTCGGTTCCAGCGACGCCATCGCCGAAGCGGGCCCGGCGGTATTACTCGCCTATCTGGTGGCCGGGCTACTGGTGGTTATCAATATGCGGATGTTGGc                                                                                                                <  2:150788/149‑1 (MQ=255)
                                     gCCACGTCACCATGTTGTCTATTGCCGGTGTTATTGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGa                                                                                                                                >  4:183144/1‑132 (MQ=255)
                                     gCCACGTCACCATGTTGTCTATTGCCGGTGTTATTGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGa                                                                                                                                <  3:183144/132‑1 (MQ=255)
                                      ccACGTCACCATGTTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCgg                                                                                                              <  6:277996/149‑1 (MQ=255)
                                            cACCATGTTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTCTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGCGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATgg                                                                                                        >  1:302960/1‑149 (MQ=255)
                                              ccATGTTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGGCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCg                                                                                                      >  8:86597/1‑149 (MQ=255)
                                                 tgttgtCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTgg                                                                                                                 >  8:304204/1‑135 (MQ=255)
                                                 tgttgtCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTgg                                                                                                                 <  7:304204/135‑1 (MQ=255)
                                                      tCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCACCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCAc                                                                                              <  3:110293/149‑1 (MQ=255)
                                                      tCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCACCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCAc                                                                                              <  3:114803/149‑1 (MQ=255)
                                                      gcTATTGCCGGGGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCAc                                                                                              <  4:198248/148‑1 (MQ=255)
                                                        tATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGaaa                                                                                                           <  2:230361/134‑1 (MQ=255)
                                                        tATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGc                                                                                            >  3:69962/1‑149 (MQ=255)
                                                        tATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGc                                                                                            >  6:117646/1‑149 (MQ=255)
                                                                   ttATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATg                                                                                                         >  1:168642/1‑125 (MQ=255)
                                                                   ttATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATg                                                                                                         <  2:168642/125‑1 (MQ=255)
                                                                     aTCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTc                                                                               >  4:23166/1‑149 (MQ=255)
                                                                      tCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCg                                                                              >  4:114345/1‑149 (MQ=255)
                                                                         gcgcAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGAGGTTATCATTATGAGGATGTTGGCGGAAATCGCCGTTGCCACACCCGATACCGGTTCGttt                                                                           >  4:104955/1‑149 (MQ=255)
                                                                           gcAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTc                                                                         >  7:323563/1‑149 (MQ=255)
                                                                                  tggttGTCGGTTCAAGCGTCGCCATCGCCGAAGCGGGCCAGGCGGTATTACTTGCCTATCTGTTCGCCGGGCTACTGGTTGTGATGATTATGCGGATGTTGGCGGAAATGGCGGGTGCCACGCCAGATACCGGTTCGTTTTCCACCTa                                                                   <  1:245635/145‑1 (MQ=255)
                                                                                    tttGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATgg                                                                                                        <  5:6942/109‑1 (MQ=255)
                                                                                    tttGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATgg                                                                                                        >  6:6942/1‑109 (MQ=255)
                                                                                       gTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGa                                                             >  6:51803/1‑149 (MQ=255)
                                                                                        tCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGTTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGAt                                                            <  7:71187/149‑1 (MQ=255)
                                                                                           gTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATaaa                                                         <  5:117646/149‑1 (MQ=255)
                                                                                               cAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGcc                                                      <  7:316641/148‑1 (MQ=255)
                                                                                                aGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCGACTGGTGGTTATGATTATGCGGATGTTGCCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCAt                                                    <  5:179876/149‑1 (MQ=255)
                                                                                                  cggcGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTGCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGGTGCCTCGCCCGATACCGGTTCGTTTTCCTCCTATGCCGATAAAGCCCTTg                                                  <  4:241537/146‑1 (MQ=255)
                                                                                                        cATCGCCGAAGTGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCt                                            <  5:128042/149‑1 (MQ=255)
                                                                                                        cATCGCCGAAGTGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCt                                            >  6:128042/1‑149 (MQ=255)
                                                                                                            gCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTgggc                                        >  5:10134/1‑149 (MQ=255)
                                                                                                              cGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTgggcgg                                      <  4:69962/149‑1 (MQ=255)
                                                                                                              cGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTgggcgg                                      >  2:64100/1‑149 (MQ=255)
                                                                                                              cGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCGACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTgggcgg                                      <  4:206519/149‑1 (MQ=255)
                                                                                                                   cGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCta                                  >  8:5214/1‑148 (MQ=255)
                                                                                                                        ccGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCAt                            >  4:134507/1‑149 (MQ=255)
                                                                                                                             ggTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTg                                           >  8:78172/1‑129 (MQ=255)
                                                                                                                             ggTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTg                                           <  7:78172/129‑1 (MQ=255)
                                                                                                                                 ttACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCTGGtt                   >  1:275092/1‑149 (MQ=255)
                                                                                                                                    cTGGCCTATCTGTTCGGCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGTCTATACCATCGGCTGGTTGTa                <  3:280825/149‑1 (MQ=255)
                                                                                                                                    cTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCTGGTTGTa                <  3:163029/149‑1 (MQ=255)
                                                                                                                                        ccTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCTGGTTGTACtgg            <  1:95224/149‑1 (MQ=255)
                                                                                                                                              tGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCTGGTTGTACTGGTGGttt      >  7:300812/1‑149 (MQ=255)
                                                                                                                                                ttCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGTAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATATAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCTGGTTGTACTGGTGGtt       <  3:104955/146‑1 (MQ=255)
                                                                                                                                                ttCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCTGGTTGTACTGGTGGtttt     <  2:275092/148‑1 (MQ=255)
                                                                                                                                                ttCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCTGGTTGTACTGGTGGtt       <  6:10134/146‑1 (MQ=255)
                                                                                                                                                  cGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCAATGGTCGCTGGGCGGGCTATACCATCGGCTGGTTGTACTGGTGGTTTTgg   <  3:23166/148‑1 (MQ=255)
                                                                                                                                                   gCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCTGGTTGTACTGGTGGTTTTggg  >  6:320356/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCGCAACCACATGAGTTAGGCGGCGGGCTGAAGTCACGCCACGTCACCATGTTGTCTATTGCCGGTGTTATCGGCGCAAGTCTGTTTGTCGGTTCCAGCGTCGCCATCGCCGAAGCGGGCCCGGCGGTATTACTGGCCTATCTGTTCGCCGGGCTACTGGTGGTTATGATTATGCGGATGTTGGCGGAAATGGCGGTTGCCACGCCCGATACCGGTTCGTTTTCCACCTATGCCGATAAAGCCATTGGTCGCTGGGCGGGCTATACCATCGGCTGGTTGTACTGGTGGTTTTGGG  >  NC_000913/2794265‑2794559

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