Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,064,592 C→A 13.1% A136S (GCC→TCC)  PP_1841 ← cytochrome c family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,064,5920CA13.1% 173.3 / 14.4 61A136S (GCC→TCC) PP_1841cytochrome c family protein
Reads supporting (aligned to +/- strand):  ref base C (28/25);  new base A (4/4);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

TGGAACAGGGCAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTCGGCGG  >  NC_002947/2064447‑2064738
                                                                                                                                                 |                                                                                                                                                  
tgGAACAGGGCAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTg                                                                                                                                                 >  5:206689/1‑149 (MQ=255)
    aCAGGGCAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTgcggc                                                                                                                                             <  1:330667/149‑1 (MQ=255)
    aCAGGGCAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTgcggc                                                                                                                                             <  6:313010/149‑1 (MQ=255)
       gggCAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCa                                                                                                                                          <  6:206689/149‑1 (MQ=255)
          cAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTc                                                                                                                                       >  1:240137/1‑149 (MQ=255)
          cAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTc                                                                                                                                       <  3:146852/149‑1 (MQ=255)
                ccTGGCGCACGGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGCCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACAcc                                                                                                                                 <  6:22511/149‑1 (MQ=255)
                 cTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCacac                                                                                                                                  >  4:307884/1‑147 (MQ=255)
                    gcgcACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGt                                                                                                                             <  3:291475/149‑1 (MQ=255)
                     cgcACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTa                                                                                                                            <  4:325994/149‑1 (MQ=255)
                      gcACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTAGCGAAGGGTGGCAATCTGGAATGCCCACTGTTAAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACCCCAgg                                                                                                                              >  8:293988/1‑146 (MQ=255)
                            ggCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCGACCGGGGACAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCTCTGGCCTGCGGCGCAGTCAACACCAGGTAGAACGGc                                                                                                                     >  8:107412/1‑149 (MQ=255)
                              cTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGccc                                                                                                                   <  7:159150/149‑1 (MQ=255)
                               tGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGACTGCGCCGCAGGCcagcg                                                                                                                                  <  7:147149/133‑4 (MQ=255)
                               tGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGACTGCGCCGCAGGCcagcg                                                                                                                                  >  8:147149/1‑130 (MQ=255)
                                            cTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGc                                                                                                     >  8:294254/1‑149 (MQ=255)
                                                  gtcgggCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGAc                                                                                               >  5:171238/1‑149 (MQ=255)
                                                        cagcagCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTgg                                                                                                          >  6:9535/1‑132 (MQ=255)
                                                                      aCCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTg                                                                           <  4:51520/149‑1 (MQ=255)
                                                                        cGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGcc                                                                         >  4:218150/1‑149 (MQ=255)
                                                                                          aGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCaa                                                       >  7:16287/1‑149 (MQ=255)
                                                                                            cTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAAc                                                      >  2:80152/1‑148 (MQ=255)
                                                                                            cTTGCGAAGGGTGGCAATATGGAATGGCCACTGTTCAGGGCGGATGCAGCTGGCCTGCGGCGAAGTCCACACCAGGTAGAACGGCCCGGAACTTGGCTTGCCTTGACCCCGTGGCGGCCACGGCTTGACTGGGTCTTCCAAGGCCAAc                                                      >  2:84670/1‑148 (MQ=255)
                                                                                                 gAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTggcggc                                                                                    >  8:280208/1‑113 (MQ=255)
                                                                                                 gAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCtt                                                                 <  5:224066/132‑1 (MQ=255)
                                                                                                 gAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCtt                                                                 >  6:224066/1‑132 (MQ=255)
                                                                                                 gAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAAc                                                      <  3:170432/143‑1 (MQ=255)
                                                                                                 gAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAAc                                                      >  4:170432/1‑143 (MQ=255)
                                                                                                 gAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCacgg                                                 >  1:122906/1‑145 (MQ=255)
                                                                                                 gAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAg                                                   <  2:122906/146‑1 (MQ=255)
                                                                                                      gTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACAAGGTAGAACGGCCAGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGttc                                           >  7:308500/1‑147 (MQ=255)
                                                                                                       tGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCt                                          >  1:41787/1‑149 (MQ=255)
                                                                                                           aaTCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGcc                                      <  8:16287/149‑1 (MQ=255)
                                                                                                                 gAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTgg                                  <  3:307884/147‑1 (MQ=255)
                                                                                                                 gAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCt                                >  3:237551/1‑149 (MQ=255)
                                                                                                                   aTGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTgc                              <  2:259419/149‑1 (MQ=255)
                                                                                                                    tGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTgca                             <  7:294254/149‑1 (MQ=255)
                                                                                                                    tGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTgca                             <  6:171238/149‑1 (MQ=255)
                                                                                                                      gCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTgcag                            <  1:80152/148‑1 (MQ=255)
                                                                                                                      gCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCATGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTgcag                            <  1:84670/148‑1 (MQ=255)
                                                                                                                       ccACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTgca                             <  8:150376/146‑1 (MQ=255)
                                                                                                                       ccACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTgca                             >  7:150376/1‑146 (MQ=255)
                                                                                                                        cACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGCCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTgcagca                          >  5:203405/1‑148 (MQ=255)
                                                                                                                         aCTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGTCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAAt                        <  7:107412/149‑1 (MQ=255)
                                                                                                                             ttCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCg                    >  2:264449/1‑149 (MQ=255)
                                                                                                                                 gggCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGc                >  1:20849/1‑149 (MQ=255)
                                                                                                                                  ggCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCg               >  7:270405/1‑149 (MQ=255)
                                                                                                                                   gCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCgg              >  8:121900/1‑149 (MQ=255)
                                                                                                                                          ccGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTc       <  2:20849/149‑1 (MQ=255)
                                                                                                                                           cGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTcg      >  2:51743/1‑149 (MQ=255)
                                                                                                                                           cGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTcg      >  5:346158/1‑149 (MQ=255)
                                                                                                                                             gtGGACTGCGCCGCAGGCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGcc                     >  5:174496/2‑132 (MQ=255)
                                                                                                                                             gtGGACTGCGCCGCAGGCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGcc                     <  6:174496/131‑1 (MQ=255)
                                                                                                                                             cTGGCCTGCGGCGCAGTCCACACCCTGTAGAAAGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTcggc    <  8:308500/149‑1 (MQ=255)
                                                                                                                                             cTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTcggc    <  2:41787/149‑1 (MQ=255)
                                                                                                                                              tGGACTGCGCCGCAGGCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTgcag                            <  2:262317/124‑1 (MQ=255)
                                                                                                                                              tGGACTGCGCCGCAGGCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTgcag                            <  2:258057/124‑1 (MQ=255)
                                                                                                                                              tGGACTGCGCCGCAGGCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTgcag                            >  1:262317/1‑124 (MQ=255)
                                                                                                                                              tGGACTGCGCCGCAGGCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTgcag                            >  1:258057/1‑124 (MQ=255)
                                                                                                                                                gCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTcggcgg  >  6:318639/1‑148 (MQ=255)
                                                                                                                                                gCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTcggcg   <  7:121900/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TGGAACAGGGCAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTCGGCGG  >  NC_002947/2064447‑2064738

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