Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,785,059 A→T 8.1% E85V (GAA→GTA) ‡ PP_3348 → GGDEF domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,785,0590AT 8.1% 213.2 / 6.9 74E85V (GAA→GTA) ‡PP_3348GGDEF domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (29/39);  new base T (4/2);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 3.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

CGTCATGCTGTGCTGGTTGGGCTGGGAGCGTTTCGACCCGCAAAGGGTCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCGGCATGTCTGTCAGCGCGGTGTC  >  NC_002947/3784916‑3785206
                                                                                                                                               |                                                                                                                                                   
cgTGATGCTGTGCTGGTTGGGCTGGGAGCGTTTCGACCCGCAAAGGGGCATGGTCTGGATGCCGTTATTGGGCGTGTCGACCGGCGTGCGTACGGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCATCCAGCGCGCGGCGGAgcgcg                                                                                                                                                >  6:329414/1‑148 (MQ=255)
  tCATGCTGTGCTGGTTGGGCTGGGAGCGTTTCGACCCGCAAAGGGTCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCTCCCGCCAACCAGCGCACGGCCGAACGCTgg                                                                                                                                              <  1:33347/149‑1 (MQ=255)
  tCATGCTGTGCTGGTTGGGCTGGGAGCGTTTCGACCCGCAAAGGGTCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTgg                                                                                                                                              <  5:224477/149‑1 (MQ=255)
  tCATGCTGTGCTGGTTGGGCTGGGAGCGTTTCGACCCGCAAAGGGTCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTgg                                                                                                                                              <  5:10443/149‑1 (MQ=255)
  tCATGCTGTGCTGGTTGGGCTGGGAGCGTTTCGACCCGCAAAGGGTCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTgg                                                                                                                                              <  6:287655/149‑1 (MQ=255)
  tCACGCAGTGCGGGCTCGGCTGCGAGCGTTTCGACCCGCAAACGGTCATGGTCGGGATGACGTAATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCggg                                                                                                                                              <  8:6696/149‑1 (MQ=255)
    aTGCTGTGCTGGTTGGGCGGGGAGCGTTTCGACCCGCAAAGGGTCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGa                                                                                                                                            <  3:289893/149‑1 (MQ=255)
     tgctgtgctgGTTGGGCTGGGAGCGTTTCGACCCGCAAAGGGTCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGTTCGGCCGTGCGCTggttg                                                                                                                                           <  3:144907/149‑4 (MQ=255)
     tgctgtgctgGTTGGGCTGGGAGCGTTTCGACCCGCAAAGGGTCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAg                                                                                                                                           <  4:348493/149‑1 (MQ=255)
       ctgtgctgGTTGGGCTGGGAGCGTTTCGACCCGCAAAGGGTCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGTTCGGCCGTGCGCTggttggc                                                                                                                                         <  3:256533/149‑6 (MQ=255)
       ctgtgctgGTTGGGCTGGGAGCGTTTCGACCCGCAAAGGGTCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGt                                                                                                                                          <  8:370532/148‑1 (MQ=255)
          tgctgGTTGGGCTGGGAGCGTTTCGACCCGCAAAGGGTCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCg                                                                                                                                      <  5:135103/149‑1 (MQ=255)
                ttGGGCTGGGAGCGTTTCGACCCGCAAAGGGTCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTg                                                                                                                                <  1:236043/149‑1 (MQ=255)
                  gggCTGGGAGCGTTTCGACCCGCAAAGGGTCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTgg                                                                                                                               <  3:10939/148‑1 (MQ=255)
                             gTTTCGAGCCGCAAAGGGTCATGGTCTGGATGACGTTATTGGGGGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGGCTATCTGCGTGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATg                                                                                                                   <  2:67166/149‑1 (MQ=255)
                             gTTTCGACCCGCAAAGGGTCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATg                                                                                                                   <  6:376938/149‑1 (MQ=255)
                             gTTTCGACCCGCAAAGGGTCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATg                                                                                                                   >  4:41302/1‑149 (MQ=255)
                              tttCGACCCGCAAAGGGTCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGc                                                                                                                  >  5:139885/1‑149 (MQ=255)
                                      cGCAAAGGGTCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTcggcg                                                                                                          <  7:366038/149‑1 (MQ=255)
                                        cAAAGGGTCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCggcggc                                                                                                        >  8:307891/1‑149 (MQ=255)
                                               tCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATctgg                                                                                                 <  5:333349/149‑1 (MQ=255)
                                               tCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATctgg                                                                                                 <  7:213962/149‑1 (MQ=255)
                                                 aTGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATctgggg                                                                                               >  1:248223/1‑149 (MQ=255)
                                                 aTGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATctggg                                                                                                <  1:170302/148‑1 (MQ=255)
                                                  tGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATctggggc                                                                                              >  3:204543/1‑149 (MQ=255)
                                                               gTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTg                                                                                 <  2:248223/149‑1 (MQ=255)
                                                               gTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCAAGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGCAGCTCTCGCGGGCTATCTGCGGCTGGGGCGCCTTTg                                                                                 >  5:129204/1‑149 (MQ=255)
                                                                ttATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCttt                                                                                  <  5:384148/147‑1 (MQ=255)
                                                                ttATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCttt                                                                                  >  6:384148/1‑147 (MQ=255)
                                                                ttATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGc                                                                                <  3:41302/149‑1 (MQ=255)
                                                                 tATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGTTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGcc                                                                               >  4:148735/1‑149 (MQ=255)
                                                                  aTTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCg                                                                              >  2:78989/1‑149 (MQ=255)
                                                                   ttGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCg                                                                              <  6:139885/148‑1 (MQ=255)
                                                                    tGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGt                                                                             >  8:62053/1‑148 (MQ=255)
                                                                     gggCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATctggggctggg                                                                                          >  1:93610/1‑134 (MQ=255)
                                                                     gggCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATctggggctggg                                                                                          <  2:93610/134‑1 (MQ=255)
                                                                      ggCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTaa                                                                           <  2:312230/148‑1 (MQ=255)
                                                                         gggtCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGcc                                                                       <  1:203541/147‑1 (MQ=255)
                                                                              gACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCg                                                                     >  5:284405/1‑146 (MQ=255)
                                                                                ccGGAGTGCGTACCGCCATGTTCTTAGCCTTTCTGGGAGCACCTGTTAACCAGTGCTCTGCCGAAGGCTGGGAGTGGGGTTACTGGGTGTGGCTGATGCTGTCTGCGGCTATCTGGGGCTGGGGCTCCTTTGCCGTAATGCCCGCAga                                                                 <  4:204543/148‑1 (MQ=255)
                                                                                    cgtgcgtACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGacgacg                                                            <  7:62053/149‑1 (MQ=255)
                                                                                          tACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGAGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTgg                                                      >  8:238864/1‑149 (MQ=255)
                                                                                             cGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGc                                                   <  3:120307/149‑1 (MQ=255)
                                                                                              gCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTAGCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCa                                                  <  7:307891/149‑1 (MQ=255)
                                                                                                cATGTTCTTCGCCTATCTGCGAGCACCCCCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGa                                                >  1:219462/1‑149 (MQ=255)
                                                                                                  tGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACg                                              >  5:115186/1‑149 (MQ=255)
                                                                                                      cttcGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGa                                          >  2:261612/1‑149 (MQ=255)
                                                                                                      cttcGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGa                                          >  4:271959/1‑149 (MQ=255)
                                                                                                         cGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATc                                        <  6:284405/148‑1 (MQ=255)
                                                                                                           ccTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGGGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATg                                     <  5:329414/149‑1 (MQ=255)
                                                                                                             tATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCt                                   >  3:131002/1‑149 (MQ=255)
                                                                                                                tgccGAACACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGtt                                >  6:293348/4‑149 (MQ=255)
                                                                                                                  gCGAGCACCCGCCAACCAGCGGACGGCCGAACGCTGGGTGGGGCGTTACTGGGGCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTg                              <  3:295303/149‑1 (MQ=255)
                                                                                                                  gCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGCCGCGCTGGGGCAGACGCTGCTCATGCTGTTTg                              >  8:163669/1‑149 (MQ=255)
                                                                                                                  gCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTg                              >  5:105713/1‑149 (MQ=255)
                                                                                                                    gAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGcc                            <  6:115186/149‑1 (MQ=255)
                                                                                                                      gCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGt                          >  2:177072/1‑149 (MQ=255)
                                                                                                                          ccGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCgg                       >  2:2371/1‑148 (MQ=255)
                                                                                                                          ccGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACACTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCgg                       >  4:303206/1‑148 (MQ=255)
                                                                                                                           cGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCGGCa                     >  1:3548/1‑149 (MQ=255)
                                                                                                                            gCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCGGCAt                    <  1:261612/149‑1 (MQ=255)
                                                                                                                            gCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCGGCAt                    <  3:271959/149‑1 (MQ=255)
                                                                                                                              ctcccAGCGTTCGGCCGTGCGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCGGCAtgt                  >  8:334270/4‑149 (MQ=255)
                                                                                                                               tcccAGCGTTCGGCCGTGCGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCGGCAtgtc                 >  5:27543/3‑149 (MQ=255)
                                                                                                                               tcccAGCGTTCGGCCGTGCGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCGGCAtgtc                 >  2:60899/3‑149 (MQ=255)
                                                                                                                                aCCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCGGCAtgtct                <  3:148735/149‑1 (MQ=255)
                                                                                                                                 ccAGCGTTCGGCCGTGCGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCGGCAtgtctg               >  4:341375/1‑149 (MQ=255)
                                                                                                                                  cAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCGGCAtgtctgt              >  5:118758/1‑149 (MQ=255)
                                                                                                                                      gcACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCGGCATGTCTGTCAgc          <  2:219462/149‑1 (MQ=255)
                                                                                                                                          ggCCGAACGCTGGGAGTGGCGTTACTGGGTCACACTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCGGCATGTCTGTCAGCGCgg      <  3:303206/149‑1 (MQ=255)
                                                                                                                                           gCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCGGCATGTCTGCCAGCGCGgt     <  1:177072/149‑1 (MQ=255)
                                                                                                                                            cctaaCGGTGGGAGTGGCGTTACTGGTTCACGGTGATGCTGTCGGCGGCTATGTGGTGCTGTGGCGCCTTTGCCGTAATGCCGGCAGACGAGGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCGGAATGTCTGTCAGCGCGgtg    <  6:105713/146‑1 (MQ=255)
                                                                                                                                             cGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCGGCATGTCTGTCAGCGCGgtgt   >  8:357170/1‑149 (MQ=255)
                                                                                                                                              gAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCGGCATGTCTGTCAGCGCGgtgt   >  2:31432/1‑148 (MQ=255)
                                                                                                                                              gAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCGGCATGTCTGTCAGCGCGGTGTc  <  7:357170/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CGTCATGCTGTGCTGGTTGGGCTGGGAGCGTTTCGACCCGCAAAGGGTCATGGTCTGGATGACGTTATTGGGCGTGTCGACCGGCGTGCGTACCGCCATGTTCTTCGCCTATCTGCGAGCACCCGCCAACCAGCGCACGGCCGAACGCTGGGAGTGGCGTTACTGGGTCACGCTGATGCTGTCGGCGGCTATCTGGGGCTGGGGCGCCTTTGCCGTAATGCCCGCAGACGACGCGCTGGGGCAGACGCTGATCATGCTGTTTGCCGTCGGCATGTCTGTCAGCGCGGTGTC  >  NC_002947/3784916‑3785206

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