Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,678,456 G→A 6.4% S10S (TCG→TCA ydgH → DUF1471 family periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,678,4560GA 6.4% 203.9 / 2.9 65S10S (TCG→TCAydgHDUF1471 family periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (26/34);  new base A (2/2);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

GGAATGTGATATGAAAAATGACGCAGACAGTTACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTC  >  NC_000913/1678310‑1678598
                                                                                                                                                  |                                                                                                                                              
ggAATGTGATATGAAAAATGACGCAGACAGTTACACCGTTTAAATGTAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCgg                                                                                                                                               >  1:284560/1‑148 (MQ=255)
ggAATGTGATATGAAAAATGACGCAGACAGTTACACCGTTTAAATGTAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGc                                                                                                                                              >  5:190834/1‑149 (MQ=255)
    tgtgATATGAAAAATGACGCAGACAGTTACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTg                                                                                                                                          >  1:254989/1‑149 (MQ=255)
    tgagATATGAAAAATGACGCAGACAGTTACACCGTTTAAATTCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGAGCAAGGCGAAGGATTATTTTTATGAACCTTAAGAACACCCTACTGGCGTCGGCACtc                                                                                                                                          >  1:254993/4‑148 (MQ=255)
        atatGAAAAATGACGCAGACAGTTACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCttt                                                                                                                                      <  2:333000/149‑1 (MQ=255)
        atatGAAAAATGACGCAGACAGTTACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCttt                                                                                                                                      <  8:244190/149‑1 (MQ=255)
         tatGAAAAATGACGCAGACAGTTACACCGGTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTc                                                                                                                                     <  4:221313/149‑1 (MQ=255)
           tGAAAAATGACGCAGACAGTTACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTg                                                                                                                                   <  1:67187/149‑1 (MQ=255)
              aaaaTGACGCAGACAGTTACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAATGATTATTTTTATGAAGCTTAAGAACACCCTCCTTGCGTCCGCACTGCTTTCTGCt                                                                                                                                 <  6:123924/148‑1 (MQ=255)
                      gCAGACAGTTACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATtc                                                                                                                        <  3:193015/149‑1 (MQ=255)
                         gACAGTTACACCGTTTACATGCACTAATCAGCCACGTTTCTCGTTAATAACACTACCAGTACCTGGGTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTGCTGCTATGGCATTCTcc                                                                                                                     <  5:232392/149‑1 (MQ=255)
                         gACAGTTACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTcc                                                                                                                     >  8:111811/1‑149 (MQ=255)
                              ttACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTaa                                                                                                                >  4:36992/1‑149 (MQ=255)
                                 cacCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCAGCACTGCTTTCTGCTATGTCATTCTCCGTTAACg                                                                                                              >  3:61143/1‑148 (MQ=255)
                                 cacCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCAGCACTGCTTTCTGCTATGGCATTCTCCGTTAACg                                                                                                              >  3:57000/1‑148 (MQ=255)
                                  acCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACg                                                                                                              >  4:47034/1‑147 (MQ=255)
                                  acCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACg                                                                                                              >  6:152479/1‑147 (MQ=255)
                                   ccGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACgcag                                                                                                           <  2:107481/149‑1 (MQ=255)
                                    cGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACgcagc                                                                                                          <  3:99646/149‑1 (MQ=255)
                                    cGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCAATGGCATTCTCCGTTAACgcag                                                                                                           >  8:312745/1‑148 (MQ=255)
                                    cGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACAATTACATGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGTCATTCTCCGTTAACgcagc                                                                                                          <  4:47906/149‑1 (MQ=255)
                                      tttAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTTTGGCATTCTCCGTTAACGCAGCaa                                                                                                        <  2:284579/149‑1 (MQ=255)
                                          aaTGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAg                                                                                                     <  2:156327/148‑1 (MQ=255)
                                            tGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAAc                                                                                                  <  8:230087/149‑1 (MQ=255)
                                             gCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACt                                                                                                 >  7:109411/1‑149 (MQ=255)
                                                  aatCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGacac                                                                                            <  8:295980/149‑1 (MQ=255)
                                                   atCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGacac                                                                                            <  3:47034/148‑1 (MQ=255)
                                                             tttCTCGTTAATAACAATACCAGTACCTGGTTTGCGCATGGCGAAGGATTATTTTTATGAAGATTAAGAACACACTCCTGGAGTCGGCAATGCTTTATGATATGGAATTATCCGTTAAAGCAGCAACATAACTGACACCGGAGCAAGCg                                                                                 >  7:326708/1‑149 (MQ=255)
                                                                ctcGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGa                                                                                               >  6:180543/1‑132 (MQ=255)
                                                                ctcGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGa                                                                                               <  5:180543/132‑1 (MQ=255)
                                                                         aaCAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTaaa                                                                     >  4:192191/1‑149 (MQ=255)
                                                                         aaCAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACGGGTTTCTGCTTTGGCATTCTCCGTTAACGCAGCAACAGAAAGGAACCCGCAGCAAGGGGCAGCGGTTaaa                                                                     >  2:89004/1‑149 (MQ=255)
                                                                          aCAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTaaa                                                                     <  6:190834/148‑1 (MQ=255)
                                                                          aCAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTaaa                                                                     <  2:284560/148‑1 (MQ=255)
                                                                              tACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCttt                                                                >  5:255507/1‑149 (MQ=255)
                                                                                ccAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCtttt                                                               <  3:192191/148‑1 (MQ=255)
                                                                                ccAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCtttt                                                               <  1:187509/148‑1 (MQ=255)
                                                                                ccAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGGCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCtttt                                                               <  5:227043/148‑1 (MQ=255)
                                                                                ccAGTACCTGGTTTGCGCAAGGCGAAGGATTAATTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCtttt                                                               <  5:100362/148‑1 (MQ=255)
                                                                                    tACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCg                                                          <  8:66272/149‑1 (MQ=255)
                                                                                    tACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCg                                                          >  1:163706/1‑149 (MQ=255)
                                                                                      ccTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCgtg                                                        <  2:163706/149‑1 (MQ=255)
                                                                                       cTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCttt                                                                >  7:64326/1‑140 (MQ=255)
                                                                                       cTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCttt                                                                <  8:64326/140‑1 (MQ=255)
                                                                                           tttGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTTCTTTCTGCTATGGAATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGtt                                                   <  5:99015/149‑2 (MQ=255)
                                                                                              gcgcAAGGCGAAGGAATATTTTTATGAAGCTTAAGAACACACTCCTGGCGTAGGCACTGCTTTCAGCTATATAATTCTCCCTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAAACTTTTGACCGTGTAGAGGt                                                 >  2:7704/1‑148 (MQ=255)
                                                                                                  aaGGCGTAGGATTATTTTTATGAAGATTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCg                                            <  4:271177/149‑1 (MQ=255)
                                                                                                    ggCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCCGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGt                                                 <  5:272942/142‑1 (MQ=255)
                                                                                                    ggCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCCGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGt                                                 >  6:272942/1‑142 (MQ=255)
                                                                                                       gAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGt                                       >  5:170958/1‑149 (MQ=255)
                                                                                                        aaGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGtt                                      <  5:152479/149‑1 (MQ=255)
                                                                                                         aGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCAGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGttt                                     <  4:61143/149‑1 (MQ=255)
                                                                                                          ggATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCAGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGttt                                     <  4:57000/148‑1 (MQ=255)
                                                                                                              tattTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATg                                >  3:296975/1‑149 (MQ=255)
                                                                                                              tattTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATg                                >  7:3592/1‑149 (MQ=255)
                                                                                                               attTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATg                                <  7:111811/148‑1 (MQ=255)
                                                                                                                tttttATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGc                               >  3:272312/1‑148 (MQ=255)
                                                                                                                     aTGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTgg                         >  1:161837/1‑149 (MQ=255)
                                                                                                                      tGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTgg                         <  6:255507/148‑1 (MQ=255)
                                                                                                                              taGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGt                >  2:520/2‑149 (MQ=255)
                                                                                                                                  acacCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGaaa            <  3:36992/149‑1 (MQ=255)
                                                                                                                                          ctGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGa                                                                                        >  6:43701/1‑65 (MQ=255)
                                                                                                                                          ctGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGa                                                                                        <  5:43701/65‑1 (MQ=255)
                                                                                                                                          ctGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTtc    <  6:238544/149‑1 (MQ=255)
                                                                                                                                             gCGTCGTCACTGCTTTCAGCTATGTCATTCTCCGTTAACGCAGCAACAGAACTGACACCTGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTtctc  <  7:308489/148‑1 (MQ=255)
                                                                                                                                             gCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTtctc  >  6:50224/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GGAATGTGATATGAAAAATGACGCAGACAGTTACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTC  >  NC_000913/1678310‑1678598

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