Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,269,968 G→T 10.9% intergenic (+366/+248) yhaC → / ← rnpB pentapetide repeats‑related protein/RNase P, M1 RNA component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,269,9680GT10.9% 166.7 / 8.9 69intergenic (+366/+248)yhaC/rnpBpentapetide repeats‑related protein/RNase P, M1 RNA component
Reads supporting (aligned to +/- strand):  ref base G (28/33);  new base T (3/5);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCA  >  NC_000913/3269825‑3270103
                                                                                                                                               |                                                                                                                                       
tAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCg                                                                                                                                    >  2:203436/1‑149 (MQ=255)
  aatcgcACTCAAAAATCTTAAACCCGCCGTATCGGTTTTTTACGTAAATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCgg                                                                                                                                   <  8:111103/145‑1 (MQ=11)
  aaGCGCACTCAACAATCTGAAACGCTGAGGATCGGTTTATTACGTTCATCAGGTAAAAGTGACCGATTAGCCGCTTTCTTTTGTGTTTAGTGTGGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAAGCCGCCGAAGCgg                                                                                                                                   <  8:438895/148‑1 (MQ=18)
  aaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCggg                                                                                                                                  <  2:925/149‑1 (MQ=255)
  aaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                   <  1:203436/148‑1 (MQ=255)
  aaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCgg                                                                                                                                   <  3:341265/148‑1 (MQ=255)
  aaGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGAAAAAACCCGCCGCAGCggg                                                                                                                                  >  7:111103/1‑149 (MQ=255)
   aGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGt                                                                                                                                 <  3:338287/149‑1 (MQ=255)
      gcACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCTTCGGCGGGtttt                                                                                                                              <  2:268362/149‑1 (MQ=255)
         cTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGAACTATAAGCCTCTTTATTTTGTGTATAGTGTCGTGGACAGTCAATCATCTATCTGCCCCTCCAAAATTAAAAACCCGCTTCGGCGGGTTTTTAc                                                                                                                           <  6:141896/149‑1 (MQ=18)
          tCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACg                                                                                                                          >  4:350811/1‑149 (MQ=255)
            aacaaTCTGAAACCCACCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGt                                                                                                                         >  3:346894/1‑148 (MQ=255)
            aacaaGCTGAAACCCGCCGGAACGGTTTATTCCGTACATCAGGTAAACCTGACCGATCAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAACAACCCGCCGAAGCGGGTTTTTACg                                                                                                                          <  1:469606/147‑1 (MQ=17)
               aaTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTaaaa                                                                                                                     >  5:266892/1‑149 (MQ=255)
                aTCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAAc                                                                                                                    <  8:432854/149‑1 (MQ=255)
                     aatcccGCCGGTCCGGGTTATTCCGTACATCAGGTAAAACTGACCTATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCGCTCCAAATGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                                                                                                 <  6:67228/144‑1 (MQ=14)
                     aaaCCCGCCGGAACGGTTTATTACGTACATCAGTTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                                <  3:350811/148‑1 (MQ=255)
                     aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAgg                                                                                                                 <  1:165294/147‑1 (MQ=255)
                     aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                                <  7:353316/148‑1 (MQ=255)
                     aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                                <  6:266892/148‑1 (MQ=255)
                     aaaCCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                                <  1:425164/148‑1 (MQ=255)
                          cgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGCCAGGCATTCATCTTTCTGCCCCGCCAAAACTAAAAACCCGCCGAAGCGGGGTTTTACGTAAAACAGGTGaga                                                                                                            >  8:272488/1‑145 (MQ=17)
                          cgccgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGt                                                                                                                <  7:272488/143‑1 (MQ=255)
                             cgGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                       <  4:346894/149‑1 (MQ=255)
                             cgGAACGGTTTATTACGTACATAAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTTCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                       <  4:451574/149‑1 (MQ=255)
                              ggAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAcc                                                                                                      <  8:213869/149‑1 (MQ=255)
                               gAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTAAAGTGTCGTGGGCAGTCATAAAACATTCAGCCCCGACAAAAGTAAAAAACAGCCGTAGCGGGTTTTTAAGTAAAACAGGTGCAACGGAcc                                                                                                      >  3:240933/1‑148 (MQ=17)
                                 aCGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGAt                                                                                                   <  4:284022/149‑1 (MQ=255)
                                                cATCAGTTAATACTGACCGATAAGCCGCTTGCTTTTGCTTTTAGTGTCGTGGACTCTGTTTCATATTTAGGCCCCTCAAAAAGTATAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTTAAACTGACCGATAAGCCGctttcttt                                                                                     <  8:360204/148‑1 (MQ=11)
                                                cATCAAGTAAAACTTACCGATAAGAAGCATTCATTTTGGTATAGTGTCGGGGACAGTCATACATCTTTCTGCCCCTCCAAAAGTATAAACCCGCCGAATCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                    <  7:203992/149‑1 (MQ=14)
                                                  tCATGTAACACTTTCCGTTCAGCCGCTTGCTGGTGGTTAGAGTGTCGGCGACTGTCTTTCCTCTTTCTGGCCCTCCAAAATTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                  <  7:448540/149‑1 (MQ=11)
                                                   ctgttAAAACTGAGCGATTAGGAGCTTGCATTTGTCGTTAGTCTCGTGGAAACTGTTTCTTATATCTGCCCCTCCAAAAGTAAAAACCCGCCGTAGCGGTTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                  <  3:202594/144‑1 (MQ=2)
                                                                           ctttctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                    <  7:57285/122‑1 (MQ=35)
                                                                           ctttctttTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                    >  8:57285/1‑122 (MQ=35)
                                                                                      gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACa                                                                 >  5:352536/1‑130 (MQ=38)
                                                                                      gTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACa                                                                 <  6:352536/130‑1 (MQ=38)
                                                                                                        tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                                                                                     <  6:481387/92‑1 (MQ=37)
                                                                                                        tcattcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                  >  8:280117/1‑95 (MQ=37)
                                                                                                        tcattcatCTTTCTGCCCCTCCAAAAGGAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                  <  6:262327/95‑1 (MQ=18)
                                                                                                        tcattcatCTTTCTGCCCCTCCAAAACTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                  >  5:262327/1‑95 (MQ=37)
                                                                                                        tcattcatCTTTCTGCCCCCCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGctttcttt                                                                                     >  5:481387/1‑92 (MQ=37)
                                                                                                          attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGAc                                                                                                       >  1:434128/1‑72 (MQ=255)
                                                                                                          attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACccgc                          >  7:326306/1‑149 (MQ=37)
                                                                                                          attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtc                                                              >  4:87504/1‑113 (MQ=17)
                                                                                                          attcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCAGAAGCGGGGTTTTACGTAAAACACGTGAAACTGACCGATAAGACGATTTCTTTTGGGTATAGTGTCGTGGACAGtc                                                              >  5:310399/1‑113 (MQ=14)
                                                                                                           ttcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcat                                                            >  2:316375/1‑114 (MQ=17)
                                                                                                            tcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                   <  8:124224/90‑1 (MQ=38)
                                                                                                            tcatCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTg                                                                                   >  7:124224/1‑90 (MQ=38)
                                                                                                             catCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGttttt                                                                                                                             >  7:186824/1‑47 (MQ=38)
                                                                                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCtttt                                                                                    >  5:40944/1‑87 (MQ=38)
                                                                                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGtcattc                                                          >  5:22731/1‑113 (MQ=17)
                                                                                                              atCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACAGATAAGCCGCTTTCtttt                                                                                    <  6:40944/87‑1 (MQ=255)
                                                                                                                      gCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGttt              >  3:221522/1‑149 (MQ=255)
                                                                                                                       ccccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                      >  6:10621/1‑124 (MQ=38)
                                                                                                                       ccccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                      <  5:10621/124‑1 (MQ=38)
                                                                                                                        cccTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGttttt            <  4:221522/149‑1 (MQ=255)
                                                                                                                           tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                  >  4:276302/1‑76 (MQ=39)
                                                                                                                           tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGGGTCGTGGACAGTCATTCATCTTTc                                                   >  3:433755/1‑107 (MQ=25)
                                                                                                                           tCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGCGTCGTGGACAGTCATTCATCTTTc                                                   <  4:433755/107‑1 (MQ=17)
                                                                                                                            ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGccc                                              >  7:254844/1‑111 (MQ=39)
                                                                                                                            ccAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGccc                                              <  8:254844/111‑1 (MQ=39)
                                                                                                                             cAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGtttttgcttt       >  8:362484/1‑144 (MQ=11)
                                                                                                                             cAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCTTCGGCGGGtttttgcttt       >  4:109775/1‑144 (MQ=11)
                                                                                                                                   tAAAAACCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                      <  6:461861/112‑1 (MQ=11)
                                                                                                                                   tAAAAACCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCaaaa                                      >  5:461861/1‑112 (MQ=11)
                                                                                                                                    aaaaaCCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCa  <  3:263167/147‑1 (MQ=255)
                                                                                                                                    aaaaaCCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCa  >  4:263167/1‑147 (MQ=255)
                                                                                                                                     aaaaCCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaaaa                                <  7:489647/116‑1 (MQ=11)
                                                                                                                                     aaaaCCCGCTTCGGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCaaaa                                >  8:489647/1‑116 (MQ=11)
                                                                                                                                            gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                  >  8:107736/1‑59 (MQ=35)
                                                                                                                                            gccgAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTgg                                                                                  <  7:107736/59‑1 (MQ=35)
                                                                                                                                               |                                                                                                                                       
TAAAGCGCACTCAACAATCTGAAACCCGCCGGAACGGTTTATTACGTACATCAGGTAAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGTAAAAACCCGCCGAAGCGGGTTTTTACGTAAAACAGGTGAAACTGACCGATAAGCCGCTTTCTTTTGGGTATAGTGTCGTGGACAGTCATTCATCTTTCTGCCCCTCCAAAAGCAAAAACCCGCCGAAGCGGGTTTTTACGTAAACCA  >  NC_000913/3269825‑3270103

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: