Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 26,844 A→T 11.9% intergenic (+40/+145) PP_0021 → / ← PP_0022 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294726,8440AT11.9% 147.6 / 9.0 57intergenic (+40/+145)PP_0021/PP_0022hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/22);  new base T (3/4);  total (31/26)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

CGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGGTGTCC  >  NC_002947/26699‑26989
                                                                                                                                                  |                                                                                                                                                  
cGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTg                                                                                                                                                 >  1:310401/1‑149 (MQ=255)
  aaTCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGt                                                                                                                                               >  8:354237/1‑149 (MQ=255)
      ccTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTa                                                                                                                                            <  7:142426/148‑1 (MQ=255)
       cTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACa                                                                                                                                          <  7:354237/149‑1 (MQ=255)
           tCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTg                                                                                                                                       >  1:128275/1‑148 (MQ=255)
           tCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTg                                                                                                                                       >  5:63617/1‑148 (MQ=255)
             cGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAAACCACTGTGGTTACAGTGaa                                                                                                                                     <  4:151069/148‑1 (MQ=255)
             cGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGaa                                                                                                                                     <  5:151418/148‑1 (MQ=255)
                     ggTTGATCATCCTTATCGCAGAATGGATAGAGTGAGGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAAACCACTGTGGTTACAGTGAAAAAAAGTgg                                                                                                                           <  7:193307/149‑1 (MQ=37)
                       ttGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGtt                                                                                                                         <  7:199957/149‑1 (MQ=255)
                             atcCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACa                                                                                                                     >  4:252633/1‑147 (MQ=255)
                              tcCTTATCGCAGAATGGATTGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACg                                                                                                                  >  4:68908/1‑149 (MQ=255)
                               ccTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGt                                                                                                                                                  >  8:143322/1‑117 (MQ=255)
                               ccTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGt                                                                                                                                                  <  7:143322/117‑1 (MQ=255)
                               ccTTATCGCAGAATGGATAGAGTAAGGTTGCCCGGGGGGGACTGGGGGGTGGGGGTGGGGGGCTGAGGGCGGGGGGGGGCTGTTGTTAACCACATTTTTTAACTGTAAACAAAGTGGGTTAAAGTGAAAAAAAg                                                                                                                              >  8:102400/1‑134 (MQ=255)
                                   aTCGCAGAATGGATAGAGTGAAGCTGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGc                                                                                                              <  4:193735/148‑1 (MQ=255)
                                      gCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTAc                                                                                                           <  1:340789/148‑1 (MQ=255)
                                      gCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTAc                                                                                                           >  2:340789/1‑148 (MQ=255)
                                      gCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTAc                                                                                                           >  3:19678/1‑148 (MQ=255)
                                       cAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTAc                                                                                                           >  1:321377/1‑147 (MQ=255)
                                           aTGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACTGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGacac                                                                                                      <  2:329410/148‑1 (MQ=255)
                                           aTGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGacac                                                                                                      <  2:310864/148‑1 (MQ=255)
                                            tGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACAcgc                                                                                                    >  3:108396/1‑149 (MQ=255)
                                            tGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACAcgc                                                                                                    >  2:314507/1‑149 (MQ=255)
                                               aTAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAAACCACTGTGGTTACAGTGaaaaaaa                                                                                                                               <  7:175704/119‑1 (MQ=25)
                                               aTAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAAACCACTGTGGTTACAGTGaaaaaaa                                                                                                                               >  8:175704/1‑119 (MQ=25)
                                               aTAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGcc                                                                                                 >  4:215220/1‑149 (MQ=255)
                                                  gagTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTg                                                                                              <  1:86983/149‑1 (MQ=255)
                                                  gagTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTg                                                                                              >  8:148674/1‑149 (MQ=255)
                                                   agTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTgg                                                                                             <  3:215220/149‑1 (MQ=255)
                                                           tGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTccc                                                                                     <  3:26158/149‑1 (MQ=255)
                                                               ccGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAgg                                                                                 >  5:224897/1‑149 (MQ=255)
                                                                cGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGc                                                                                <  2:128275/149‑1 (MQ=255)
                                                                  ggaggaACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCa                                                                              <  2:321377/149‑1 (MQ=255)
                                                                       aaCTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAgcgcg                                                                         >  6:111393/1‑149 (MQ=255)
                                                                                  ctgctgGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGa                                                              >  5:347800/1‑149 (MQ=255)
                                                                                    gctgGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGAtt                                                            >  7:325823/1‑149 (MQ=255)
                                                                                     ctgGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCAtt                                                                <  7:148674/144‑1 (MQ=255)
                                                                                     ctgGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGa                                                              <  6:63617/146‑1 (MQ=255)
                                                                                     ctgGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGAt                                                             <  4:108396/147‑1 (MQ=255)
                                                                                     ctgGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTc                                                           <  5:32983/149‑1 (MQ=255)
                                                                                              aCTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTc                                                   <  6:224897/148‑1 (MQ=255)
                                                                                               cTCAGGCATGAGGAGGCCTGTTGTAAAGCACTTTTTTTCACTGTATCCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCACAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCg                                                  <  8:325823/148‑1 (MQ=255)
                                                                                                tCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCggg                                                >  1:13204/1‑149 (MQ=255)
                                                                                                   ggCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGt                                             >  7:32771/1‑149 (MQ=255)
                                                                                                    gCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTc                                            <  4:311240/149‑1 (MQ=255)
                                                                                                        gaggagGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAac                                        >  6:13442/1‑149 (MQ=255)
                                                                                                         aggaggCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGaa                                         <  3:252633/147‑1 (MQ=255)
                                                                                                                  gttgtAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTAtgctg                              >  4:314753/1‑149 (MQ=255)
                                                                                                                   ttgtAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGTTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTAtgctgc                             >  2:248019/1‑149 (MQ=255)
                                                                                                                       aaaCCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAac                                        <  5:26612/134‑1 (MQ=255)
                                                                                                                       aaaCCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAac                                        >  6:26612/1‑134 (MQ=255)
                                                                                                                           cACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGcc                     >  2:79237/1‑149 (MQ=255)
                                                                                                                                      aCTGTAAACCACTGTGGTTACAGTGAAGAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCAcgcc           >  1:16607/1‑149 (MQ=38)
                                                                                                                                      aCTGTAAACCACTGTGGTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCAcgccg          >  5:142267/1‑149 (MQ=255)
                                                                                                                                      aCTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCAcgccg          >  1:350398/1‑149 (MQ=255)
                                                                                                                                      aCTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCAcgccg          >  1:149144/1‑149 (MQ=255)
                                                                                                                                               cacaGTGGTTTACAGTGAAAAAAATTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGGTGTcc  <  2:310401/149‑1 (MQ=255)
                                                                                                                                               cacaGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGGTGTcc  <  5:111393/149‑1 (MQ=255)
                                                                                                                                               cacaGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGGTGTcc  <  1:314507/149‑1 (MQ=255)
                                                                                                                                               cacaGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGGTGTcc  <  1:148456/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGAATCCCTTATCCGGATGCCGGTTGATCATCCTTATCGCAGAATGGATAGAGTGAAGATGGCCCGGGAGGAACTGGGCGATCTGCTGGAGCTGACTCAGGCCTGAGGAGGCCTGTTGTAAACCACTTTTTTTCACTGTAACCACAGTGGTTTACAGTGAAAAAAAGTGGTTTACAACGTAGCTACGACACGCGCCCTGGCCACTCCCCAGGCCAGCGCGCTCCTCATTGATTCGTTCAGTCGGGCGTCGAACACCTATGCTGCTCGGGGCCTTGCCACGCCGCGGTGTCC  >  NC_002947/26699‑26989

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: