Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,224,391 A→C 7.8% I52M (ATT→ATG PP_1962 ← site‑specific recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,224,3910AC 7.8% 195.8 / 5.8 63I52M (ATT→ATGPP_1962site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base A (30/28);  new base C (2/3);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-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)
       aaaCCGATAGTCGTGAGATTTCGCTTGTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCAACCTGGTATCTCTCGATCAAAGCGAAt                                                                                                                                                >  1:171876/1‑139 (MQ=255)
       aaaCCGATAGTCGTGAGATTTCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAAt                                                                                                                                                <  2:171876/138‑1 (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                                                                                                                   <  1:275317/149‑1 (MQ=255)
                         tttCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCaa                                                                                                                   <  2:236269/149‑1 (MQ=255)
                         tttCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGAGCTTTGTTTTGTGTCCGCaa                                                                                                                   >  3:200633/1‑149 (MQ=255)
                         tttCGCTTTGGATGATACGAAGGAACGAATACGATTGAAGGGACTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGc                                                                                                                     >  5:320678/1‑147 (MQ=255)
                                         tACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTAc                                                                                                                                           <  6:78626/110‑1 (MQ=255)
                                         tACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTAc                                                                                                                                           >  5:78626/1‑110 (MQ=255)
                                         tACGAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTg                                                                                                     >  6:164778/1‑148 (MQ=255)
                                            gAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGtctc                                                                                                 >  1:122254/1‑149 (MQ=255)
                                            gAAGGAACGAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGCTCGTTCCCACCTGGTATCTGTCGATCAAAGCGAATGTTACATGCGCTTTTTGTTGTTTCCTCAAGTTTCCAAATGATGtctc                                                                                                 >  8:122454/1‑149 (MQ=255)
                                                    gAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGc                                                                                                                     <  2:312326/121‑1 (MQ=255)
                                                    gAATACGATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGc                                                                                                                     >  1:312326/1‑121 (MQ=255)
                                                          gATTGAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTTCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGActct                                                                                   <  3:172996/149‑1 (MQ=255)
                                                              gAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAAt                                                                                <  6:320678/148‑1 (MQ=255)
                                                              gAAGGGAGTAAATGCAGGGAGCAAATGGTAAACAATTGCCCTTTTCAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAAt                                                                                >  2:164994/1‑148 (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                                          >  3:273731/1‑141 (MQ=255)
                                                                                                           cAGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTAAGAAGCtt                                          <  4:273731/141‑1 (MQ=255)
                                                                                                            aGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTaa                                                 <  1:243703/133‑1 (MQ=255)
                                                                                                            aGGTCGTTCCCACCTGGTATCTCTCGATCAAAGCGAATGTTACATGCGCTTTGTGTTGTGTCCGCAAGTTTCCAAATGCTGTCTCTGAATTCCGACTCTAATGAGACGGGGAAATCTCCGCCTGGGAGCGTaa                                                 >  2:243703/1‑133 (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: