Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,064,598 G→C 13.6% P134A (CCG→GCG)  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,5980GC13.6% 157.6 / 16.1 59P134A (CCG→GCG) PP_1841cytochrome c family protein
Reads supporting (aligned to +/- strand):  ref base G (28/23);  new base C (4/4);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.73e-01

ACAGGGCAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTCGGCGGCAAA  >  NC_002947/2064451‑2064742
                                                                                                                                                   |                                                                                                                                                
aCAGGGCAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTgcggc                                                                                                                                                 <  6:313010/149‑1 (MQ=255)
aCAGGGCAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTgcggc                                                                                                                                                 <  1:330667/149‑1 (MQ=255)
   gggCAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCa                                                                                                                                              <  6:206689/149‑1 (MQ=255)
      cAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTc                                                                                                                                           <  3:146852/149‑1 (MQ=255)
      cAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTc                                                                                                                                           >  1:240137/1‑149 (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                                 <  6:171238/149‑1 (MQ=255)
                                                                                                                tGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTgca                                 <  7:294254/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          >  5:346158/1‑149 (MQ=255)
                                                                                                                                       cGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTcg          >  2:51743/1‑149 (MQ=255)
                                                                                                                                         gtGGACTGCGCCGCAGGCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGcc                         <  6:174496/131‑1 (MQ=255)
                                                                                                                                         gtGGACTGCGCCGCAGGCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGcc                         >  5:174496/2‑132 (MQ=255)
                                                                                                                                         cTGGCCTGCGGCGCAGTCCACACCCTGTAGAAAGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTcggc        <  8:308500/149‑1 (MQ=255)
                                                                                                                                         cTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTcggc        <  2:41787/149‑1 (MQ=255)
                                                                                                                                          tGGACTGCGCCGCAGGCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTgcag                                <  2:262317/124‑1 (MQ=255)
                                                                                                                                          tGGACTGCGCCGCAGGCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTgcag                                >  1:258057/1‑124 (MQ=255)
                                                                                                                                          tGGACTGCGCCGCAGGCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTgcag                                >  1:262317/1‑124 (MQ=255)
                                                                                                                                          tGGACTGCGCCGCAGGCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTgcag                                <  2:258057/124‑1 (MQ=255)
                                                                                                                                            gCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTcggcgg      >  6:318639/1‑148 (MQ=255)
                                                                                                                                            gCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTcggcg       <  7:121900/147‑1 (MQ=255)
                                                                                                                                               tGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTCGGCGGCa    >  1:270190/1‑147 (MQ=255)
                                                                                                                                                 cggcgCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTCGGCGGCaaa  <  2:240137/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ACAGGGCAAAGCCCTGGCGCACCGGGCTGTCGGCCGGCAGCTTCGGGTCGGGCAGCAGCGCAGGAAACCGCTGGTCCACCGGGGCCAGCTTGCGAAGGGTGGCAATCTGGAATGGCCACTGTTCAGGGCGGATGCCGCTGGCCTGCGGCGCAGTCCACACCAGGTAGAACGGCCCGGCACTTGGCTTGCCTTGACCCAGTGGCGGCCACGGCTTGCCTGGGTCTTCCACGGCCAACCAGGCCCGTGCTGGCCCTGGCTGCAGCAATGCCGCAGCGGGCATTTCGGCGGCAAA  >  NC_002947/2064451‑2064742

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