Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,224,396 T→A 8.1% N51Y (AAC→TAC)  PP_1962 ← site‑specific recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,224,3960TA 8.1% 176.4 / 5.4 62N51Y (AAC→TAC) PP_1962site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base T (29/28);  new base A (2/3);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01

TTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGTTCTCCGGCG  >  NC_002947/2224249‑2224535
                                                                                                                                                   |                                                                                                                                           
tttCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGtt                                                                                                                                             >  3:265257/1‑148 (MQ=255)
    gATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGCATGTAACAtt                                                                                                                                        >  5:210229/1‑148 (MQ=255)
     aTAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCATTTTCAGGTCGTTCACACCTGGTATATCTAGATCAAAGCCAATGTTACATga                                                                                                                                       >  5:109975/1‑148 (MQ=255)
     aTAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCACTTTTCAGGTCGTTCCAACCTGGTATCTCTCGATCAAAGCGAATGTTACATga                                                                                                                                       >  5:143454/1‑148 (MQ=255)
          ccGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTATACAATTGCCCTATTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCttt                                                                                                                                  >  7:189987/1‑149 (MQ=255)
          ccGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCttt                                                                                                                                  >  7:185181/1‑149 (MQ=255)
          ccGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCttt                                                                                                                                  <  7:6427/149‑1 (MQ=255)
          ccGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCttt                                                                                                                                  <  5:277465/149‑1 (MQ=255)
            gATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTttgt                                                                                                                                >  2:121137/1‑149 (MQ=255)
             aTAGTCGTGAGATTTCGCTTTGGGTGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGTGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGCATGTAACATTCGCTttgat                                                                                                                               <  7:209795/149‑3 (MQ=255)
              tAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGCATGTAACATTCGCTttgatc                                                                                                                              <  6:210229/149‑4 (MQ=255)
                gTCGTGAGATTTCTCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTttgtgtt                                                                                                                             <  1:121137/148‑1 (MQ=255)
                gTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTttgtgtt                                                                                                                             <  2:34147/148‑1 (MQ=255)
                 tCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTttgtgttgt                                                                                                                           >  2:275317/1‑149 (MQ=255)
                         tttCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCaa                                                                                                                   <  2:236269/149‑1 (MQ=255)
                         tttCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCaa                                                                                                                   <  1:275317/149‑1 (MQ=255)
                         tttCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGAGCTTTGTTTTGTGTCCGCaa                                                                                                                   >  3:200633/1‑149 (MQ=255)
                         tttCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGACTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGc                                                                                                                     >  5:320678/1‑147 (MQ=255)
                                        tACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTAc                                                                                                                                           >  5:78626/1‑110 (MQ=255)
                                        tACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTAc                                                                                                                                           <  6:78626/110‑1 (MQ=255)
                                        tACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTg                                                                                                     >  6:164778/1‑148 (MQ=255)
                                           gAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGtctc                                                                                                 >  1:122254/1‑149 (MQ=255)
                                           gAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGCTCGTTCCCACCTGGTATCTGTCGATCAAAGCGAATGTTACATGCGCTTTTTGTTGTTTCCTCAAGTTTCCAAATGATGtctc                                                                                                 >  8:122454/1‑149 (MQ=255)
                                                   gAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGc                                                                                                                     >  1:312326/1‑121 (MQ=255)
                                                   gAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGc                                                                                                                     <  2:312326/121‑1 (MQ=255)
                                                         gATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTTCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGActct                                                                                   <  3:172996/149‑1 (MQ=255)
                                                             gAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAAt                                                                                >  2:164994/1‑148 (MQ=255)
                                                             gAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAAt                                                                                <  6:320678/148‑1 (MQ=255)
                                                             gAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAAt                                                                                <  1:164994/148‑1 (MQ=255)
                                                              aaGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATg                                                                               <  8:185181/148‑1 (MQ=255)
                                                              aaGGGAGTAAATCCAGGTAGCAAATCGTAAACTATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTGTCGATCAATGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCATGTTTGCATATGCTGTCTCTGAATTCCGTCTCTAATg                                                                               <  8:189987/148‑1 (MQ=255)
                                                                  gAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACg                                                                          >  6:82081/1‑149 (MQ=255)
                                                                         tGCAGGGAGCAAATTGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAAt                                                                   <  1:5773/149‑1 (MQ=255)
                                                                                   aaaTGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTgg                                                         >  3:304095/1‑149 (MQ=255)
                                                                                     aTGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGa                                                       <  5:164778/149‑1 (MQ=255)
                                                                                          aaaCAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGt                                                   >  8:165298/1‑148 (MQ=255)
                                                                                           aaCAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGt                                                   >  3:136173/1‑147 (MQ=255)
                                                                                              aaTTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATACCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTaagaa                                              <  4:265257/149‑1 (MQ=255)
                                                                                                  gCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGa                                                                                              <  4:288534/97‑1 (MQ=255)
                                                                                                  gCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGa                                                                                              >  3:288534/1‑97 (MQ=255)
                                                                                                   cccTTTTCAGTTCGTGCCCACCTGGTATCTCTCGATCAAAGCGAATGGTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCtt                                          <  8:170971/148‑1 (MQ=255)
                                                                                                    ccTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGa                                        <  7:165298/149‑1 (MQ=255)
                                                                                                      ttttCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGAATCTGAATTCCGACTCTAATGAGAAGTGGAAATCTACGACTAGGAGCGTAAGAAGCTTGAt                                       >  7:235131/1‑148 (MQ=255)
                                                                                                        ttCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATaa                                     <  5:82081/148‑1 (MQ=255)
                                                                                                          cAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCtt                                          <  4:273731/141‑1 (MQ=255)
                                                                                                          cAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCtt                                          >  3:273731/1‑141 (MQ=255)
                                                                                                           aGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTaa                                                 >  2:243703/1‑133 (MQ=255)
                                                                                                           aGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTaa                                                 <  1:243703/133‑1 (MQ=255)
                                                                                                                      aCCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTc                       >  8:185032/1‑148 (MQ=255)
                                                                                                                       ccTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGa                     <  4:304095/149‑1 (MQ=255)
                                                                                                                         tGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAg                    >  2:154504/1‑148 (MQ=255)
                                                                                                                           gTATCTCTCGATCAAATCGAATGTTACATGAGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCTACTCTAATGAGACGGGGAAATCTCCGCGTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGAc                 <  6:109975/149‑1 (MQ=255)
                                                                                                                           gTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGa                  <  1:154504/148‑1 (MQ=255)
                                                                                                                                tctcGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAg            >  8:172551/1‑149 (MQ=255)
                                                                                                                                  tcGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCaa             <  4:136173/146‑1 (MQ=255)
                                                                                                                                  tcGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGt           >  4:196452/1‑148 (MQ=255)
                                                                                                                                      tCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGTTCTcc      >  6:61327/1‑149 (MQ=255)
                                                                                                                                        aaaGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGTTCTCcgg    >  5:146577/1‑149 (MQ=255)
                                                                                                                                          aGCGCATGTAACATTCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGTTCTcccggg  >  1:281281/1‑145 (MQ=255)
                                                                                                                                          aGCGCATGTAACATTCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGTTCTCcggc   <  2:281281/148‑1 (MQ=255)
                                                                                                                                          aGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGTTCTCcggcg  <  3:196452/149‑1 (MQ=255)
                                                                                                                                              aaTGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCaa             >  7:27371/1‑134 (MQ=255)
                                                                                                                                              aaTGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCaa             <  8:27371/134‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                           
TTTCGATAAACCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCTTGATAAGTAGGCCGCTCTTCGAGGACTCAAGTTCTCCGGCG  >  NC_002947/2224249‑2224535

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