Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,383,235 A→G 6.6% R14R (AGA→AGG)  PP_5547 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,383,2350AG 6.6% 174.1 / 2.1 61R14R (AGA→AGG) PP_5547hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/30);  new base G (2/2);  total (29/32)
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.79e-01

GTGCGCCGGCGACGGAGAAGGTCATTAGTGCGATAGCTTCGAAAAGCCATCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGTGCGCTTGCCTGGCTGGGGAAACCTCACGCCCAAGAGGCAGTCGCCAAGCT  >  NC_002947/3383087‑3383378
                                                                                                                                                    |                                                                                                                                               
gTGCGCCGGCGACGGAGAAGGTCATTAGTGCGATAGCTTCGAAAAGCCATCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCgaga                                                                                                                                                 <  5:350668/149‑1 (MQ=255)
             ggAGAAGGTCATTAGTGCGATAGCTTCGAAAAGCCATCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAg                                                                                                                                    >  1:140983/1‑149 (MQ=255)
             ggAGAAGGTCATTAGTGCGATAGCTTCGAAAAGCCATCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAg                                                                                                                                    <  3:95393/149‑1 (MQ=255)
             ggAGAAGGTCATTAGTGCGATAGCTTCGAAAAGCCATCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAg                                                                                                                                    <  2:203148/149‑1 (MQ=255)
             ggAGAAGGTCATTAGTGCGATAGCTTCGAAAAGCCATCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAg                                                                                                                                    >  7:128553/1‑149 (MQ=255)
                 aaGGTCATTAGTGCGATAGCTTCGAAAAGCCATCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCaa                                                                                                                                <  5:9050/149‑1 (MQ=255)
                 aaGGTCATTAGTGCGATAGCTTCGAAAAGCCATCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCaa                                                                                                                                <  5:3736/149‑1 (MQ=255)
                 aaGGTCATTAGTGCGATAGCTTCGAAAAGCCATCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCaa                                                                                                                                <  2:140983/149‑1 (MQ=255)
                    gacATTAGTGCGATATCTTCGAAAAGCCATCAGGTAATCGCTTTAATCGGTGCCCGCGCAGCGAACCTACTAGCGCTAACGCAATAGGTGTCGATTAGCGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCAGCACTTAGGCAAgg                                                                                                                              >  6:337514/3‑148 (MQ=255)
                       aTTAGTGCGATAGCTTCGAAAAGCCATCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTg                                                                                                                                           >  1:121878/1‑132 (MQ=255)
                       aTTAGTGCGATAGCTTCGAAAAGCCATCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTg                                                                                                                                           <  2:121878/132‑1 (MQ=255)
                              cGATAGCTTCGAAAAGCCATCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGtaat                                                                                                                    <  3:91907/148‑1 (MQ=255)
                                 tAGCTTCGAAAAGCCATCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGGACCCTGCAGTTAGGCAAGGCTCAGGTAATAATc                                                                                                                >  1:246113/1‑149 (MQ=255)
                                 tAGCTTCGAAAAGCCATCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGGACCCTGCAGTTAGGCAAGGCTCAGGTAATAATc                                                                                                                >  1:240816/1‑149 (MQ=255)
                                   gCTTCGAAAAGCCATCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGa                                                                                                              <  5:314529/149‑1 (MQ=255)
                                     ttCGAAAAGCCATCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAg                                                                                                             >  8:7995/1‑148 (MQ=255)
                                        gAAAAGCCATCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCAt                                                                                                          <  5:7293/148‑1 (MQ=255)
                                         aaaaGCCATCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATgc                                                                                                        >  1:209323/1‑149 (MQ=255)
                                                aTCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAg                                                                                                 >  1:17599/1‑149 (MQ=255)
                                                 tCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGt                                                                                                <  7:240123/149‑1 (MQ=255)
                                                   aGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTgg                                                                                              >  5:188263/1‑149 (MQ=255)
                                                     gTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGAGGc                                                                                             >  3:295080/1‑148 (MQ=255)
                                                        aTCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCaaa                                                                                          <  5:346910/148‑1 (MQ=255)
                                                        aTCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAAt                                                                                         <  5:93256/149‑1 (MQ=255)
                                                        aTCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAAt                                                                                         >  5:5245/1‑149 (MQ=255)
                                                          cGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCg                                                                                       >  1:192536/1‑149 (MQ=255)
                                                           gCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCg                                                                                       <  6:188263/148‑1 (MQ=255)
                                                                   cGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAg                                                                              <  2:209323/149‑1 (MQ=255)
                                                                    ggTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGAGGCAAATCGCTCTCGGAGc                                                                             <  4:295080/149‑1 (MQ=255)
                                                                        cccGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGAc                                                                          >  1:298215/1‑148 (MQ=255)
                                                                        cccGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGAc                                                                          <  1:137770/148‑1 (MQ=255)
                                                                           gcgcAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAAACCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCAtt                                                                      >  7:93582/1‑149 (MQ=255)
                                                                                 cGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGc                                                                >  8:99915/1‑149 (MQ=255)
                                                                                    aCCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCa                                                                                                                                         <  2:292844/73‑1 (MQ=255)
                                                                                    aCCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCa                                                                                                                                         >  1:292844/1‑73 (MQ=255)
                                                                                     ccTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGAc                                                             >  7:147501/1‑148 (MQ=255)
                                                                                            gggCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGct                                                                                     >  4:169617/1‑117 (MQ=255)
                                                                                            gggCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGct                                                                                     <  3:169617/117‑1 (MQ=255)
                                                                                             ggCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGt                                                    <  7:34365/149‑1 (MQ=255)
                                                                                              gCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGTg                                                   <  5:337514/149‑1 (MQ=255)
                                                                                                 aaCGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGa                                                                                                              >  8:102530/1‑87 (MQ=255)
                                                                                                 aaCGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGa                                                                                                              <  7:102530/87‑1 (MQ=255)
                                                                                                    gCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTg                                                                  >  2:294483/1‑128 (MQ=255)
                                                                                                    gCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTg                                                                  <  1:294483/128‑1 (MQ=255)
                                                                                                      aaTAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGTGCGCTTGcc                                           <  1:366506/149‑1 (MQ=255)
                                                                                                              tCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGGACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGTGCGCTTGCCTGGCTggg                                   <  2:246113/149‑1 (MQ=255)
                                                                                                              tCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGGACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGTGCGCTTGCCTGGCTggg                                   <  2:240816/149‑1 (MQ=255)
                                                                                                              tCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGt                                                    >  1:168982/1‑132 (MQ=255)
                                                                                                              tCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGt                                                    <  2:168982/132‑1 (MQ=255)
                                                                                                              tCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGTGCGCTTGCctggctgg                                    <  2:192536/148‑1 (MQ=255)
                                                                                                                gATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGTGCGCtt                                              >  7:204414/1‑136 (MQ=255)
                                                                                                                gATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGTGCGCtt                                              <  8:204414/136‑1 (MQ=255)
                                                                                                                               tcttACCGGCTGACGGACGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATTCGCCTCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGAGACGCTGTACGTGCGCTTGCCTGGCTGGGGAAACCTCACGCCCaa                   <  8:318350/146‑1 (MQ=255)
                                                                                                                               tGTTACCGGCGGACGGCAGAGAACCCTGCAGTGAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGTGCGCTTGCCTGGCTGGGGAAACCTCACGCCCaa                   <  8:323273/148‑1 (MQ=255)
                                                                                                                                gTTACCGGCGGACGGCCGAGCACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGGGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGTGCGCTTGCCTGGCTGGGGAAACCTCCCGCCCaa                   >  7:323273/1‑147 (MQ=255)
                                                                                                                                gTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGATGGAGACGCTGTACGTGAGCTTGCCTGGCTGGGGAAACCTAACGCCCaa                   >  7:318350/1‑147 (MQ=255)
                                                                                                                                 ttACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGTGCGCTTGCCTGGCTGGGGAAACCTCACGCCCAagag                >  6:129480/1‑149 (MQ=255)
                                                                                                                                    ccGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGTGCGCTTGCCTGGCTGGGGAAACCTCACGcc                      >  3:64398/1‑140 (MQ=255)
                                                                                                                                    ccGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGTGCGCTTGCCTGGCTGGGGAAACCTCACGcc                      <  4:64398/140‑1 (MQ=255)
                                                                                                                                       gcggACGGCCGAGAACCCTGCAGTCAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGTGCGCTTGCCTGGCTGGGGAAACCTCACGCCCAAGAGGCAGTc          >  6:101087/1‑149 (MQ=255)
                                                                                                                                             ggCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGTGCGCTTGCCTGGCTGGGGAAACCTCACGCCCAAGAGGCAGTCGCCAAg    >  4:322624/1‑149 (MQ=255)
                                                                                                                                               ccGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGTGCGCTTGCCTGGCTGGGGAAACCTCACGCCCAAGAGGCAGTCGCCAAGCt  <  2:298215/149‑1 (MQ=255)
                                                                                                                                                  agaACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGTGCGCTTGCctggctg                                     <  6:339446/111‑1 (MQ=255)
                                                                                                                                                  agaACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGTGCGCTTGCctggctg                                     >  5:339446/1‑111 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GTGCGCCGGCGACGGAGAAGGTCATTAGTGCGATAGCTTCGAAAAGCCATCAGGTAATCGCTTTAAGCGGTGCCCGCGCAGCGAACCTTCTAGGGCTAACGCAATAGGTGTCGATTAGGGAAGTATTTGTTACCGGCGGACGGCCGAGAACCCTGCAGTTAGGCAAGGCTCAGGTAATAATCGAGCATGCGCCGCAGTGGCAAATCGCTCTCGGAGCGACCATTGCTGGCGACGCTGTACGTGCGCTTGCCTGGCTGGGGAAACCTCACGCCCAAGAGGCAGTCGCCAAGCT  >  NC_002947/3383087‑3383378

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