Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,722,903 G→C 6.1% D54H (GAC→CAC)  PP_3288 → universal stress protein family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,722,9030GC 6.1% 315.2 / 7.4 98D54H (GAC→CAC) PP_3288universal stress protein family protein
Reads supporting (aligned to +/- strand):  ref base G (46/46);  new base C (2/4);  total (48/50)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

AAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTGGGTGCAGCAGCCCTTGAATGAGATCCTGAT  >  NC_002947/3722757‑3723048
                                                                                                                                                  |                                                                                                                                                   
aaaGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGTTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAATGGCTTTCCACGGCGGGCCTGGTGAATGAc                                                                                                                                                   <  6:220703/149‑1 (MQ=255)
aaaGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGAc                                                                                                                                                   <  3:387785/149‑1 (MQ=255)
     ctgctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTCATTCACCAGgc                                                                                                                                              <  1:356868/149‑1 (MQ=255)
        ctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCt                                                                                                                                           <  7:188549/149‑1 (MQ=255)
        ctgctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTCATTCACCAGgcc                                                                                                                                             <  6:306802/147‑2 (MQ=255)
           ctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTgg                                                                                                                                        >  3:318708/1‑148 (MQ=255)
              aTCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCa                                                                                                                                      <  3:127918/147‑1 (MQ=255)
               tCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTc                                                                                                                                   <  1:236576/149‑1 (MQ=255)
                cGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTCATTCACCAGgccc                                                                                                                                            <  2:368285/140‑3 (MQ=255)
                cGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTCATTCACCAGgccc                                                                                                                                            >  1:368285/1‑138 (MQ=255)
                 gCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCa                                                                                                                                  <  3:366707/148‑1 (MQ=255)
                    tCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATgcg                                                                                                                              >  7:137752/1‑149 (MQ=255)
                       ccccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAATGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGa                                                                                                                           <  4:309877/149‑1 (MQ=255)
                       ccccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATgcgc                                                                                                                             >  3:427274/1‑147 (MQ=255)
                       ccccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATgcgc                                                                                                                             >  2:113419/1‑147 (MQ=255)
                         ccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAgg                                                                                                                         <  7:162422/149‑1 (MQ=255)
                         ccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAgg                                                                                                                         >  7:288589/1‑149 (MQ=255)
                             aTGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTa                                                                                                                     >  2:246803/1‑149 (MQ=255)
                                 gcCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCAt                                                                                                                                 <  6:300283/133‑1 (MQ=255)
                                 gcCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCAt                                                                                                                                 >  5:300283/1‑133 (MQ=255)
                                 gcCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTg                                                                                                                 >  4:155904/1‑149 (MQ=255)
                                 gcCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTg                                                                                                                 <  8:288589/149‑1 (MQ=255)
                                    gAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGCGCAGTCATGCGCGAGGGTTACGTgca                                                                                                               <  1:5324/149‑1 (MQ=255)
                                     cagcccGGTCATAGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcag                                                                                                              >  1:219321/4‑148 (MQ=255)
                                     aaCCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGgc                                                                                                                                              <  6:93755/117‑1 (MQ=255)
                                     aaCCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGgc                                                                                                                                              >  5:93755/1‑117 (MQ=255)
                                      agcccGGTCATAGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcag                                                                                                              >  3:324977/3‑147 (MQ=255)
                                       gcccGGTCATAGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgaagcag                                                                                                           >  8:219466/2‑149 (MQ=255)
                                        cccGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcagcag                                                                                                           >  7:360375/1‑148 (MQ=255)
                                             tCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTCATTCACCAGgcc                                                                                                                                             >  6:430412/1‑109 (MQ=255)
                                             tCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTCATTCACCAGgcc                                                                                                                                             <  5:430412/110‑2 (MQ=255)
                                              cTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGt                                                                                                     <  8:395897/148‑1 (MQ=255)
                                               tATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTg                                                                                                    >  7:327611/1‑148 (MQ=255)
                                                  gACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAAt                                                                                                 <  2:219321/148‑1 (MQ=255)
                                                  gACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATg                                                                                                >  7:90281/1‑149 (MQ=255)
                                                  gACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGACTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATg                                                                                                >  3:54110/1‑149 (MQ=255)
                                                   aCCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcagca                                                                                                            >  3:127365/1‑136 (MQ=255)
                                                   aCCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTgcagca                                                                                                            <  4:127365/136‑1 (MQ=255)
                                                    ccGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGc                                                                                              <  4:318708/149‑1 (MQ=255)
                                                     cGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCt                                                                                             <  1:113566/149‑1 (MQ=255)
                                                     cGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCt                                                                                             >  2:113566/1‑149 (MQ=255)
                                                     cGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACGACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCGCCGTGAATGGCt                                                                                             >  8:67267/1‑149 (MQ=255)
                                                          tgctgcGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTg                                                                                                                                          >  7:282930/1‑100 (MQ=255)
                                                          tgctgcGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTg                                                                                                                                          <  8:282930/100‑1 (MQ=255)
                                                               cgcTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCa                                                                                   <  7:110883/148‑1 (MQ=255)
                                                               cgcTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAg                                                                                  >  5:241968/1‑149 (MQ=255)
                                                               cgcTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAg                                                                                  >  1:9033/1‑149 (MQ=255)
                                                               cgcTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAg                                                                                  >  1:221523/1‑149 (MQ=255)
                                                                cttGGCCTTGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAgg                                                                                 >  4:337682/3‑149 (MQ=255)
                                                                     ccaaggccaaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTg                                                                            >  4:145439/1‑149 (MQ=255)
                                                                      caaggccaaggGCATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGc                                                                           <  2:9033/149‑1 (MQ=255)
                                                                      caaggccaaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAgg                                                                                                                         <  6:85391/104‑1 (MQ=255)
                                                                      caaggccaaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAgg                                                                                                                         >  5:85391/1‑104 (MQ=255)
                                                                      caaggccaaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGc                                                                           <  2:4843/149‑1 (MQ=255)
                                                                      caaggccaaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGc                                                                           <  4:427274/149‑1 (MQ=255)
                                                                       aaggccaaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTAGGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGcc                                                                          >  5:372792/1‑149 (MQ=255)
                                                                       aaggccaaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGcc                                                                          <  8:137752/149‑1 (MQ=255)
                                                                         ggccaaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTTCCACGGCTGGCCTGGTGAATGACCAGGCGCTGTCAGTCATGCGCGTGGGTTACGTGCATCAGCACCGTGAATTGCTTGAAACCCAGGCGTTGccta                                                                        <  4:410867/149‑3 (MQ=255)
                                                                         ggccaaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCATGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGa                                                                        <  7:344910/149‑1 (MQ=255)
                                                                            caaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATgc                                                                     >  3:29322/1‑149 (MQ=255)
                                                                             aaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATgcg                                                                    >  7:107983/1‑149 (MQ=255)
                                                                               gggAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATgcggc                                                                  <  8:53088/149‑1 (MQ=255)
                                                                                ggAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATgcggcg                                                                 <  3:155904/149‑1 (MQ=255)
                                                                                ggAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGACACGGAGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACAGTGAATGGCTGGAAACCCAGGCGTTGCCGATgcggcg                                                                 >  8:201699/1‑149 (MQ=255)
                                                                                 gAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATgcggcg                                                                 >  2:239348/1‑148 (MQ=255)
                                                                                      gcgcTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCTGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATg                                                            <  3:229876/148‑1 (MQ=255)
                                                                                      gcgcTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATgg                                                           >  2:178942/1‑149 (MQ=255)
                                                                                       cgcTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGaa                                                                                        >  3:271518/1‑119 (MQ=255)
                                                                                       cgcTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGaa                                                                                        <  4:271518/119‑1 (MQ=255)
                                                                                       cgcTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGgt                                                          <  3:337682/149‑1 (MQ=255)
                                                                                            cacaTCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCAc                                                     >  8:134130/1‑149 (MQ=255)
                                                                                             acaTCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCaac                                                    >  5:357296/1‑147 (MQ=255)
                                                                                             acaTCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCAcc                                                    <  3:75200/149‑1 (MQ=255)
                                                                                                 cgtcgCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTca                                                >  1:69281/1‑149 (MQ=255)
                                                                                                  gtcgCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTca                                                >  6:87035/1‑148 (MQ=255)
                                                                                                  gtcgCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTca                                                <  2:221523/148‑1 (MQ=255)
                                                                                                    cgCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACaa                                             <  8:360375/149‑1 (MQ=255)
                                                                                                     gCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACaa                                             <  6:241968/148‑1 (MQ=255)
                                                                                                      cTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAgg                                                                                                                         <  2:333886/72‑1 (MQ=255)
                                                                                                      cTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAgg                                                                                                                         >  1:333886/1‑72 (MQ=255)
                                                                                                        ttCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAAcc                                           >  2:52213/1‑147 (MQ=255)
                                                                                                             cTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGTAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAggtg                                     <  1:113419/149‑1 (MQ=255)
                                                                                                                cGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGgtgt                                 <  8:327611/149‑1 (MQ=255)
                                                                                                                cGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGCGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGgtgt                                 <  1:354871/149‑1 (MQ=255)
                                                                                                                  tGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTgg                               <  2:263616/149‑1 (MQ=255)
                                                                                                                   tgAATGGCTTGCCTCGTCGGGCCTGGGTAATGACCAGGCGCTGGATGTCATGCGCGATGGTGACGTGCAGCAGCACCGTGATTGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGgtgt                                 <  4:324977/145‑1 (MQ=255)
                                                                                                                   ggAAGGGCTTGCCACGGTGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTggg                              >  1:261176/1‑149 (MQ=255)
                                                                                                                   ggAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTggg                              >  7:39854/1‑149 (MQ=255)
                                                                                                                   ggAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTggg                              >  3:115478/1‑149 (MQ=255)
                                                                                                                    gAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAg                                        >  7:377441/1‑138 (MQ=255)
                                                                                                                    gAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAg                                        <  8:377441/138‑1 (MQ=255)
                                                                                                                    gAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTggg                              <  3:145439/148‑1 (MQ=255)
                                                                                                                      aGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTGGGTg                            >  8:123130/1‑148 (MQ=255)
                                                                                                                                 cggcggGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTGGGTGCAGCAGCCCtt                 >  3:411659/1‑148 (MQ=255)
                                                                                                                                       gCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTGGGTGCAGCAGCCCTTGAATgag          <  4:115478/149‑1 (MQ=255)
                                                                                                                                        ccTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTGGGTGCAGCAGCCCTTGAATgaga         <  4:54110/149‑1 (MQ=255)
                                                                                                                                            gTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTGGGTGCAGCAGCCCTTGAATGAGATCCt     >  7:323416/1‑149 (MQ=255)
                                                                                                                                               aaTGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTGGGTGCAGCAGCCCTTGAATGAGATCCTGAt  <  6:372792/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGGCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTGGGTGCAGCAGCCCTTGAATGAGATCCTGAT  >  NC_002947/3722757‑3723048

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: