Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,240,567 C→G 5.5% H43Q (CAC→CAG)  PP_4618 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,240,5670CG 5.5% 236.7 / 3.8 73H43Q (CAC→CAG) PP_4618hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/31);  new base G (2/2);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

TGTGGAGATACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAGGTAC  >  NC_002947/5240425‑5240715
                                                                                                                                              |                                                                                                                                                    
tgtgGAGATACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTg                                                                                                                                                >  3:86231/1‑149 (MQ=255)
     agaTACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAgcg                                                                                                                                           >  2:204567/1‑149 (MQ=255)
         aCCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAgcgcagc                                                                                                                                       <  7:70612/149‑4 (MQ=255)
         aCCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAgcgc                                                                                                                                          >  2:61179/1‑146 (MQ=255)
          ccccATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAgcgcagc                                                                                                                                       <  7:75774/148‑4 (MQ=255)
             cATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCTCGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGAt                                                                                                                                     >  7:201533/1‑147 (MQ=255)
               tGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCGTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAgc                                                                                                                                 <  2:108695/149‑1 (MQ=255)
                 aaGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAgcg                                                                                                                                <  3:71745/148‑1 (MQ=255)
                 aaGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAgcg                                                                                                                                <  3:307508/148‑1 (MQ=255)
                      aCGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGAc                                                                                                                                                   >  7:277727/1‑124 (MQ=255)
                      aCGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGAc                                                                                                                                                   <  8:277727/124‑1 (MQ=255)
                         tCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGc                                                                                                                       <  8:201533/149‑1 (MQ=255)
                          cAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGcc                                                                                                                      >  6:199633/1‑149 (MQ=255)
                             gcctgcctggcctgcctgGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGAt                                                                                                                   <  4:86231/149‑1 (MQ=255)
                               ctgccgggcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCATCCTGCCGATGc                                                                                                                 >  7:346210/1‑149 (MQ=255)
                                 gcctggcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCtt                                                                                                               >  6:373013/1‑149 (MQ=255)
                                  cctggcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTc                                                                                                              >  2:140720/1‑149 (MQ=255)
                                   ctggcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGcgcg                                                                                                                                         >  2:189819/1‑121 (MQ=255)
                                   ctggcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGcgcg                                                                                                                                         <  1:189819/121‑1 (MQ=255)
                                      gcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCa                                                                                                            <  6:384830/147‑1 (MQ=255)
                                      gcctgcctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCAc                                                                                                           <  5:199633/148‑1 (MQ=255)
                                           cctggAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGc                                                                                                      <  5:92382/148‑1 (MQ=255)
                                                                 ctgctgGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCAGGTCGTGCAgcgca                                                                                                                                         >  4:71092/1‑90 (MQ=255)
                                                                 ctgctgGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCAGGTCGTGCAgcgca                                                                                                                                         <  3:71092/91‑2 (MQ=255)
                                                                 ctgctgGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTgcg                                                                               >  2:288114/1‑149 (MQ=255)
                                                                 ctgctgGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTgcg                                                                               >  4:152663/1‑149 (MQ=255)
                                                                  tgctgGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTgcgc                                                                              >  1:82566/1‑149 (MQ=255)
                                                                  tgctgGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTgcgc                                                                              >  3:69193/1‑149 (MQ=255)
                                                                          gCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAAGCTACCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGaa                                                                      >  1:311063/1‑149 (MQ=255)
                                                                           ccGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAAt                                                                     <  1:140720/149‑1 (MQ=255)
                                                                                 tgtgGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTg                                                               >  6:152458/1‑149 (MQ=255)
                                                                                   tgGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTgg                                                              <  1:61179/148‑1 (MQ=255)
                                                                                   tgGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTg                                                               <  4:69193/147‑1 (MQ=255)
                                                                                     gATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCTAACTGGCCCAGCCACTGCTGCGCCAGCTGAATgcgc                                                                 >  8:15703/1‑143 (MQ=255)
                                                                                     gATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCTAACTGGCCCAGCCACTGCTGCGCCAGCTGAATgcgc                                                                 <  7:15703/143‑1 (MQ=255)
                                                                                        tgctgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCttt                                                         >  8:344961/1‑148 (MQ=255)
                                                                                        tgctgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCttt                                                         >  8:131900/1‑148 (MQ=255)
                                                                                          ctgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcag                                                      <  5:373013/149‑1 (MQ=255)
                                                                                          ctgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcag                                                      <  7:131900/149‑1 (MQ=255)
                                                                                          ctgcCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcag                                                      >  7:126100/1‑149 (MQ=255)
                                                                                             ccGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcagcag                                                   <  5:152458/149‑1 (MQ=255)
                                                                                             ccGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcagc                                                     <  1:288114/147‑1 (MQ=255)
                                                                                                aGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGAt                                                >  6:237789/1‑149 (MQ=255)
                                                                                                  cACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGa                                              <  7:344961/149‑1 (MQ=255)
                                                                                                    cGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTggg                                                             >  4:287131/1‑132 (MQ=255)
                                                                                                    cGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTggg                                                             <  3:287131/132‑1 (MQ=255)
                                                                                                    cGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAat                                            <  3:152663/149‑1 (MQ=255)
                                                                                                     gACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACtgctgc                                                                                >  6:117387/1‑112 (MQ=255)
                                                                                                     gACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACtgctgc                                                                                <  5:117387/112‑1 (MQ=255)
                                                                                                     gACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAata                                           >  6:230675/1‑149 (MQ=255)
                                                                                                     gACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAat                                            >  5:234693/1‑148 (MQ=255)
                                                                                                        ccTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGATGGGCTTTCAGCAGGATGAATATCa                                        >  8:245074/1‑149 (MQ=255)
                                                                                                                  aGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATg                                               >  1:317278/1‑132 (MQ=255)
                                                                                                                  aGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATg                                               <  2:317278/132‑1 (MQ=255)
                                                                                                                  aGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTgcgc                              >  8:268259/1‑149 (MQ=255)
                                                                                                                      cGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcagca                                                    >  2:259145/1‑123 (MQ=255)
                                                                                                                      cGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTcagca                                                    <  1:259145/123‑1 (MQ=255)
                                                                                                                      cGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGcccc                           >  7:109836/1‑148 (MQ=255)
                                                                                                                       gCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCActgctg                                                                                 >  2:172703/1‑93 (MQ=255)
                                                                                                                       gCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCActgctg                                                                                 <  1:172703/93‑1 (MQ=255)
                                                                                                                        cTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCAt                        <  5:237789/149‑1 (MQ=255)
                                                                                                                         tGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCAccc                                     >  2:360553/1‑135 (MQ=255)
                                                                                                                         tGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCAccc                                     <  1:360553/135‑1 (MQ=255)
                                                                                                                         tGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATg                       >  4:329076/1‑149 (MQ=255)
                                                                                                                             aGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCa                                        >  2:337288/1‑128 (MQ=255)
                                                                                                                             aGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCa                                        <  1:337288/128‑1 (MQ=255)
                                                                                                                                cTTGCTGCGCTGCACAACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCgat                >  4:89426/1‑149 (MQ=255)
                                                                                                                                   gctgcgctgcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCgatgat             >  3:143410/1‑149 (MQ=255)
                                                                                                                                    ctgcgctgcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGatgatg            <  6:234693/149‑1 (MQ=255)
                                                                                                                                     tgcgctgcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGatgatg            >  2:161810/1‑148 (MQ=255)
                                                                                                                                       cgctgcAGGTCGTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGa                                                 <  4:187896/109‑1 (MQ=255)
                                                                                                                                       cgctgcAGGTCGTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGa                                                 >  3:187896/1‑109 (MQ=255)
                                                                                                                                       cgctgcACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGa          >  3:153071/1‑148 (MQ=255)
                                                                                                                                            cacaACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAgg     <  3:89426/148‑1 (MQ=255)
                                                                                                                                              cGACCTGCAGCGCGAGATAAGCGCCATCCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAGGTAc  <  8:346210/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TGTGGAGATACCCCATGAAGCCACGTCAAGCCTGCCTGGCCTGCCTGGAACGCGAGCCTGTCGACCTGCTGGAAGCCGCACTGTGGATTGCTGCCGAGCACGACCCTGGTGTCGAGCCCGCTGCCAGCCTTGCTGCGCTGCACGACCTGCAGCGCGAGATAAGCGCCAACCTGCCGATGCTTCCACTATGCGAACTGGCCCAGCCACTGCTGCGCCAGCTGAATGCGCTGGGCTTTCAGCAGGATGAATATCACCCCTTGCGCCCCCATGCGGCGATGATGGACAAGGTAC  >  NC_002947/5240425‑5240715

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑