Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,691,171 T→A 10.9% noncoding (170/184 nt) glmY ← sRNA activator of glmS mRNA, glmZ processing antagonist

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,691,1710TA10.9% 197.0 / 15.2 64noncoding (170/184 nt)glmYsRNA activator of glmS mRNA, glmZ processing antagonist
Reads supporting (aligned to +/- strand):  ref base T (30/27);  new base A (4/3);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.88e-01

AGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGG  >  NC_000913/2691028‑2691312
                                                                                                                                                |                                                                                                                                             
agcagAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTAccc                                                                                                                                          >  4:625029/1‑149 (MQ=255)
 gcagAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCg                                                                                                                                         >  2:686062/1‑149 (MQ=255)
    gAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCt                                                                                                                                      <  1:239181/149‑1 (MQ=255)
     aaTCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCGGCtt                                                                                                                                     <  4:252702/149‑1 (MQ=255)
     aaTCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCtt                                                                                                                                     <  4:186326/149‑1 (MQ=255)
      aTCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATT‑CCCGGCttt                                                                                                                                    <  2:163717/149‑1 (MQ=255)
       tCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGAGAAAAACAGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCGTGATGCTAACCGAGGGGAAGTTCAGGTACAACAAAGCCGGG‑AATCACCCGGCTTTg                                                                                                                                   >  4:427379/1‑149 (MQ=255)
            aaaaaTGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATg                                                                                                                              >  4:531887/1‑149 (MQ=255)
                aTGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAAt                                                                                                                          >  4:347619/1‑149 (MQ=255)
                      gCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCg                                                                                                                    >  4:406422/1‑149 (MQ=255)
                         ttACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCGGCTTTGtt                                                                                                                                 >  1:229952/1‑133 (MQ=255)
                         ttACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCGGCTTTGtt                                                                                                                                 <  2:229952/133‑1 (MQ=255)
                         ttACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTg                                                                                                                 >  3:620736/1‑149 (MQ=255)
                          tACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGc                                                                                                                <  3:190599/149‑1 (MQ=255)
                           aCCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTAt                                                                                                                               >  3:355868/1‑133 (MQ=255)
                           aCCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTAt                                                                                                                               <  4:355868/133‑1 (MQ=255)
                               gTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGc                                                                                                                <  3:641510/144‑1 (MQ=255)
                               gTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGc                                                                                                                >  4:641510/1‑144 (MQ=255)
                                   tCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGa                                                                                                       >  2:300713/1‑149 (MQ=255)
                                    cGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGAc                                                                                                      <  3:531887/149‑1 (MQ=255)
                                     gTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACg                                                                                                     >  4:660842/1‑149 (MQ=255)
                                      tAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCGGCtt                                                                                                                                     >  4:561592/1‑116 (MQ=255)
                                      tAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑TAATTCCCGGCtt                                                                                                                                     <  3:561592/116‑1 (MQ=255)
                                      tAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTAt                                                                                                                               <  4:324126/122‑1 (MQ=255)
                                      tAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTAt                                                                                                                               >  3:324126/1‑122 (MQ=255)
                                       aaTGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGtt                                                                                                   >  3:340360/1‑149 (MQ=255)
                                          gAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCg                                                                                                >  3:358888/1‑149 (MQ=255)
                                                    aaaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                      <  3:625029/149‑1 (MQ=255)
                                                      aCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGa                                                                                    <  1:300713/149‑1 (MQ=255)
                                                             aGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCt                                                                             >  1:596445/1‑149 (MQ=255)
                                                               cGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTAc                                                                           >  3:495571/1‑149 (MQ=255)
                                                                     aaGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGat                                                                     <  4:340360/149‑1 (MQ=255)
                                                                     aaGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGat                                                                     >  1:252859/1‑149 (MQ=255)
                                                                     aaGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGat                                                                     <  3:646803/149‑1 (MQ=255)
                                                                            tACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGaa                                                              <  2:596445/149‑1 (MQ=255)
                                                                            tACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGaa                                                              <  3:498916/149‑1 (MQ=255)
                                                                               tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGtt                                                           >  4:243206/1‑149 (MQ=255)
                                                                               tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGtt                                                           >  1:185867/1‑149 (MQ=255)
                                                                               tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGtt                                                           >  1:570693/1‑149 (MQ=255)
                                                                               tctcGTCAAACGCGCATAGCCTGATGCTAACTGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCCGTTCTGTTATGGAATTAGGCGCTGCCTAACTCGAAGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGtt                                                           <  3:427379/149‑1 (MQ=255)
                                                                                   gTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATAGAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGAc                                                       >  3:395589/1‑149 (MQ=255)
                                                                                      aGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggc                                                    <  2:630622/149‑1 (MQ=255)
                                                                                       gACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggca                                                   >  3:98973/1‑149 (MQ=255)
                                                                                          gcgAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTa                                                                            >  4:238508/1‑121 (MQ=255)
                                                                                          gcgAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTa                                                                            <  3:238508/121‑1 (MQ=255)
                                                                                                  ccTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCttgtt                                        <  3:347619/149‑1 (MQ=255)
                                                                                                  ccTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCttgtt                                        <  3:243206/149‑1 (MQ=255)
                                                                                                  ccTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCttgtt                                        <  1:686062/149‑1 (MQ=255)
                                                                                                  ccTGATGATAACCGAGGGGAAGTTCAGATACAACAAAGCCAGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCttgtt                                        <  3:73779/149‑1 (MQ=255)
                                                                                                    tGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTtgttgt                                      <  4:620736/149‑1 (MQ=255)
                                                                                                    tGATGCTAACCGAGGGGAAGTTCAGAGACAACAAAGCCAGG‑AATTACCCGGCTTTGTTATGGAATACGGCGGTGCCTAACTCGACGTTTCGCCCGATGCTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTtgttgt                                      >  3:429676/1‑149 (MQ=255)
                                                                                                      aTGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGc                                    <  1:515527/149‑1 (MQ=255)
                                                                                                         cTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTc                                 >  3:209008/1‑149 (MQ=255)
                                                                                                           aaCCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCAt                               <  3:225901/149‑1 (MQ=255)
                                                                                                              cGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCg                            <  2:570693/149‑1 (MQ=255)
                                                                                                                                  acaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTccc        >  3:575175/1‑149 (MQ=255)
                                                                                                                                  acaaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTccc        >  1:162503/1‑149 (MQ=255)
                                                                                                                                   caaAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCg       >  3:681906/1‑149 (MQ=255)
                                                                                                                                    aaaGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGa      >  4:156995/1‑149 (MQ=255)
                                                                                                                                     aaGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGaa     <  4:358888/149‑1 (MQ=255)
                                                                                                                                     aaGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGaa     >  4:330306/1‑149 (MQ=255)
                                                                                                                                       gCCGGG‑TAATTCCCGGCTTTGTTGTGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAgg   >  3:1082/1‑149 (MQ=255)
                                                                                                                                       gCCGGG‑TAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAgg   >  1:615355/1‑149 (MQ=255)
                                                                                                                                        ccGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAggg  <  4:681906/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                             
AGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGG‑AATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGG  >  NC_000913/2691028‑2691312

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑