Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,499,480 2 bp→AC intergenic (+57/‑106) trpS → / → zapE tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,499,4800CA96.7% 176.5 / NA 61intergenic (+57/‑107)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (24/35);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0029471,499,4810AC98.4% 180.0 / NA 61intergenic (+58/‑106)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (25/35);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGT‑GT‑T‑TCGGGCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAG  >  NC_002947/1499337‑1499614
                                                                                                                                                ||                                                                                                                                         
gCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACACGCGCCGCGGGGCATGGCACCGGc                                                                                                                                       >  8:176755/1‑149 (MQ=39)
   gACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCttt                                                                                                                                    <  6:349894/149‑1 (MQ=255)
      gggCCGCAGTTCGACAACGCTGAAGCCCGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGcc                                                                                                                                 >  6:165754/1‑149 (MQ=255)
        gCCGCAGTTCGACAACGCTGAAGCCCGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCgg                                                                                                                               <  5:365779/149‑1 (MQ=255)
        gCCGCAGTTCGACAACGCTGAAGCCCGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCgg                                                                                                                               >  1:230491/1‑149 (MQ=255)
         ccGCAGTTCGACAACGCTGAAGCCCGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgt                                                                                                                              <  2:288972/149‑1 (MQ=255)
          cGCAGTTCGACAACGCTGAAGCCCGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACACCAGGACTGAGGCTGCCAAACGGTCGGAGCGGGCGTGACCGCGAAACACGCGCCGCCGTGCAAGGCACCGCCATTGCCGgt‑g                                                                                                                            >  2:244646/1‑149 (MQ=255)
              gTTCGACAACGCTGAAGCCCGTGATGGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACACGCGCCCCGGTGCATGGCCACGGCGTTGCCGgt‑gt                                                                                                                           >  1:3036/1‑147 (MQ=39)
               ttCGACAACGCTGAAGCCCGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCTGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑Gt‑tcgc                                                                                                                      <  1:81219/149‑4 (MQ=255)
                  gACAACGCTGAAGCCCGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑Gt‑tcgcggg                                                                                                                   >  5:279932/1‑143 (MQ=255)
                     aaCGCTGAAGCCCGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                                <  2:138359/135‑1 (MQ=255)
                     aaCGCTGAAGCCCGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                                >  1:138359/1‑135 (MQ=255)
                       cGCTGAAGCCCGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCAcg‑c                                                                                                             <  1:150689/149‑1 (MQ=37)
                         cTGAAGCCCGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCAcg‑ccc                                                                                                           <  6:363195/149‑2 (MQ=37)
                          tGAAGCCCGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGGCCGCGAAACAGGCGCCCCCGGGCATGGCACCCGGTTTGGCGGT‑GT‑TAGCGGGCCcg‑cccg                                                                                                          >  1:136248/1‑147 (MQ=255)
                              gCCCGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACCGGCGCCGCGGTGCATGGCCCCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTcc                                                                                                       >  6:66359/1‑149 (MQ=25)
                              gCCCGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTcc                                                                                                       <  7:245238/149‑1 (MQ=25)
                              gCCCGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGGGGGAGCGGGCGGGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTcc                                                                                                       <  6:278801/149‑1 (MQ=25)
                                 cGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCca                                                                                                     >  6:87907/1‑148 (MQ=25)
                                 cGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCca                                                                                                     >  5:87562/1‑148 (MQ=25)
                                     aT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGAc                                                                                                <  5:292494/149‑1 (MQ=21)
                                        gCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCCCAGACt                                                                                               >  8:4904/1‑148 (MQ=21)
                                           gCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAAGGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGAtt                                                                                           <  6:165496/149‑1 (MQ=21)
                                             cATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGc                                                                                         <  5:165754/149‑1 (MQ=21)
                                              aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgt                                                                                                                              <  1:162961/113‑1 (MQ=255)
                                              aTTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgt                                                                                                                              >  2:162961/1‑113 (MQ=255)
                                                   gagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACAGGCGCCCCGGTGCCTGGCACCGGCGTTGGCGGT‑GT‑TAGCGGGCACGCCCGCGCCCCCAGACAGATTGCCAtt                                                                                     >  2:220603/1‑147 (MQ=255)
                                                       cttaaCCAAGCCTTGCCACCCCCGCCGGACTGAGGCTGCCAAACGGGGGGAGCGGGGGGGCCCGCGAAACAGCCGCCGCGGTGCAGGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                               <  3:316582/146‑1 (MQ=18)
                                                       cTGAACGAAGCCTTGGCACCACTGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGACCTGGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGTGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                               <  1:244646/149‑1 (MQ=21)
                                                          aaCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                             <  5:298614/148‑1 (MQ=21)
                                                               aGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                       <  5:87907/149‑1 (MQ=21)
                                                                  cTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTa                                                                    <  6:87562/149‑1 (MQ=21)
                                                                       caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgt‑gt                                                                                                                           >  8:331716/1‑90 (MQ=255)
                                                                       caccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgt‑gt                                                                                                                           <  7:331716/90‑1 (MQ=255)
                                                                        accacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGa                                                              >  1:48128/1‑149 (MQ=21)
                                                                         ccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCAcg‑cc                                                                                                            <  5:126586/101‑1 (MQ=18)
                                                                           acaGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACgg                                                           >  5:89360/1‑149 (MQ=21)
                                                                                       cggCTGCCAAACGGTGGGAGCGGGCGGGCCCGCGAAACAGGCGCCGCGGTGCAGGGCACCGGCTTTGCCGGT‑GTGC‑GCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgt                                               <  3:336998/148‑1 (MQ=25)
                                                                                       aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGTATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttg                                                <  5:20116/148‑1 (MQ=21)
                                                                                       aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                 <  6:89360/147‑1 (MQ=21)
                                                                                       aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                                 <  5:167507/147‑1 (MQ=21)
                                                                                        ggCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                              <  2:48128/149‑1 (MQ=21)
                                                                                           tGCCAAACGGTGGGAGCGGGCGTGCCCGCCAAAAAGGCGCCGCGGTGCATGGCACCCGCTTTGTAGGT‑GT‑TAGCGGGCACGCCCCCTCCCACAGACAGATTGCCATTAAGCGGGCCGATCGCGACAGTGACGGCCCGGTTTtttttggc                                           >  8:305093/1‑147 (MQ=255)
                                                                                             ccAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGGT‑C‑GCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGcc                                           <  7:176755/147‑1 (MQ=21)
                                                                                             ccAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTc                                         <  1:220603/149‑1 (MQ=21)
                                                                                             ccAAACGGGGGGCGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTGGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTc                                         <  3:311597/149‑1 (MQ=21)
                                                                                                   ggTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAAc                                   <  1:224437/149‑1 (MQ=21)
                                                                                                       ggAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt                               <  7:111619/149‑1 (MQ=25)
                                                                                                          gCGGGCGTGCCCGCTAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                            <  2:230491/149‑1 (MQ=25)
                                                                                                             ggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGAc                         >  8:251146/1‑149 (MQ=25)
                                                                                                                  gCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                    >  3:150709/1‑149 (MQ=37)
                                                                                                                   cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                    <  6:279932/148‑1 (MQ=37)
                                                                                                                     cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                    >  6:287195/1‑146 (MQ=37)
                                                                                                                     cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg                 <  2:3036/149‑1 (MQ=37)
                                                                                                                         aaaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACg             >  4:112128/1‑149 (MQ=37)
                                                                                                                          aaCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCCCGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGc            >  5:232363/1‑149 (MQ=37)
                                                                                                                           aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGc            >  3:151891/1‑148 (MQ=37)
                                                                                                                                cgccgcGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCaa      <  2:136248/149‑1 (MQ=37)
                                                                                                                                 gccgcGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg     >  3:328815/1‑149 (MQ=37)
                                                                                                                                  ccgcGGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGc    <  7:4904/149‑1 (MQ=37)
                                                                                                                                    gcgGTGCATGGCACCGGCTTTGCCGGT‑GT‑TCGCGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCat  >  1:355100/1‑148 (MQ=37)
                                                                                                                                                ||                                                                                                                                         
GCCGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGAT‑GCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGT‑GT‑T‑TCGGGCACG‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAGCAG  >  NC_002947/1499337‑1499614

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: ‑