Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,722,901 A→T 6.1% N53I (AAT→ATT)  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,9010AT 6.1% 289.1 / 6.6 98N53I (AAT→ATT) PP_3288universal stress protein family protein
Reads supporting (aligned to +/- strand):  ref base A (46/46);  new base T (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 = 6.81e-01

CCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTGGGTGCAGCAGCCCTTGAATGAGATCCTGAT  >  NC_002947/3722753‑3723048
                                                                                                                                                    |                                                                                                                                                     
ccGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGaa                                                                                                                                                       <  7:258313/149‑1 (MQ=255)
 cGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAAt                                                                                                                                                      <  5:107525/149‑1 (MQ=255)
    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)
               ctgATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTg‑g                                                                                                                                        >  3:318708/1‑148 (MQ=255)
                  aTCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCa                                                                                                                                      <  3:127918/147‑1 (MQ=255)
                   tCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTc                                                                                                                                   <  1:236576/149‑1 (MQ=255)
                    cGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTCATTCACCAGgccc                                                                                                                                            >  1:368285/1‑138 (MQ=255)
                    cGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTCATTCACCAGgccc                                                                                                                                            <  2:368285/140‑3 (MQ=255)
                     gCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCa                                                                                                                                  <  3:366707/148‑1 (MQ=255)
                        tCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATgcg                                                                                                                              >  7:137752/1‑149 (MQ=255)
                           ccccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAATGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGa                                                                                                                           <  4:309877/149‑1 (MQ=255)
                           ccccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATgcgc                                                                                                                             >  3:427274/1‑147 (MQ=255)
                           ccccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATgcgc                                                                                                                             >  2:113419/1‑147 (MQ=255)
                             ccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAgg                                                                                                                         >  7:288589/1‑149 (MQ=255)
                             ccTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAgg                                                                                                                         <  7:162422/149‑1 (MQ=255)
                                 aTGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTa                                                                                                                     >  2:246803/1‑149 (MQ=255)
                                     gcCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCAt                                                                                                                                 <  6:300283/133‑1 (MQ=255)
                                     gcCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCAt                                                                                                                                 >  5:300283/1‑133 (MQ=255)
                                     gcCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTg                                                                                                                 >  4:155904/1‑149 (MQ=255)
                                     gcCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTg                                                                                                                 <  8:288589/149‑1 (MQ=255)
                                        gAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTGCGCAGTCATGCGCGAGGGTTACGTgca                                                                                                               <  1:5324/149‑1 (MQ=255)
                                         cagcccGGTCATAGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTgcag                                                                                                              >  1:219321/4‑148 (MQ=255)
                                         aaCCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGgc                                                                                                                                              <  6:93755/117‑1 (MQ=255)
                                         aaCCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGgc                                                                                                                                              >  5:93755/1‑117 (MQ=255)
                                          agcccGGTCATAGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTgcag                                                                                                              >  3:324977/3‑147 (MQ=255)
                                           gcccGGTCATAGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTgaagcag                                                                                                           >  8:219466/2‑149 (MQ=255)
                                            cccGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTgcagcag                                                                                                           >  7:360375/1‑148 (MQ=255)
                                                 tCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTCATTCACCAGgcc                                                                                                                                             >  6:430412/1‑109 (MQ=255)
                                                 tCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTCATTCACCAGgcc                                                                                                                                             <  5:430412/110‑2 (MQ=255)
                                                  cTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGt                                                                                                     <  8:395897/148‑1 (MQ=255)
                                                   tATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTg                                                                                                    >  7:327611/1‑148 (MQ=255)
                                                      gACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAAt                                                                                                 <  2:219321/148‑1 (MQ=255)
                                                      gACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATg                                                                                                >  7:90281/1‑149 (MQ=255)
                                                      gACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGACTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATg                                                                                                >  3:54110/1‑149 (MQ=255)
                                                       aCCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTgcagca                                                                                                            <  4:127365/136‑1 (MQ=255)
                                                       aCCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTgcagca                                                                                                            >  3:127365/1‑136 (MQ=255)
                                                        ccGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGc                                                                                              <  4:318708/149‑1 (MQ=255)
                                                         cGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCt                                                                                             >  2:113566/1‑149 (MQ=255)
                                                         cGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCt                                                                                             <  1:113566/149‑1 (MQ=255)
                                                         cGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACGACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCGCCGTGAATGGCt                                                                                             >  8:67267/1‑149 (MQ=255)
                                                              tgctgcGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTg                                                                                                                                          <  8:282930/100‑1 (MQ=255)
                                                              tgctgcGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTg                                                                                                                                          >  7:282930/1‑100 (MQ=255)
                                                                   cgcTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCa                                                                                   <  7:110883/148‑1 (MQ=255)
                                                                   cgcTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAg                                                                                  >  1:221523/1‑149 (MQ=255)
                                                                   cgcTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAg                                                                                  >  5:241968/1‑149 (MQ=255)
                                                                   cgcTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAg                                                                                  >  1:9033/1‑149 (MQ=255)
                                                                    cttGGCCTTGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAgg                                                                                 >  4:337682/3‑149 (MQ=255)
                                                                         ccaaggccaaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTg                                                                            >  4:145439/1‑149 (MQ=255)
                                                                          caaggccaaggGCATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGc                                                                           <  2:9033/149‑1 (MQ=255)
                                                                          caaggccaaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAgg                                                                                                                         <  6:85391/104‑1 (MQ=255)
                                                                          caaggccaaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAgg                                                                                                                         >  5:85391/1‑104 (MQ=255)
                                                                          caaggccaaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGc                                                                           <  4:427274/149‑1 (MQ=255)
                                                                          caaggccaaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGc                                                                           <  2:4843/149‑1 (MQ=255)
                                                                           aaggccaaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTAGGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGcc                                                                          >  5:372792/1‑149 (MQ=255)
                                                                           aaggccaaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGcc                                                                          <  8:137752/149‑1 (MQ=255)
                                                                             ggccaaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTTCCACGGCTGGCCTGGTGAATGACCAGGCGCTG‑TCAGTCATGCGCGTGGGTTACGTGCATCAGCACCGTGAATTGCTTG‑AAACCCAGGCGTTGccta                                                                        <  4:410867/149‑3 (MQ=255)
                                                                             ggccaaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCATGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGa                                                                        <  7:344910/149‑1 (MQ=255)
                                                                                caaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATgc                                                                     >  3:29322/1‑149 (MQ=255)
                                                                                 aaggGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATgcg                                                                    >  7:107983/1‑149 (MQ=255)
                                                                                   gggAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATgcggc                                                                  <  8:53088/149‑1 (MQ=255)
                                                                                    ggAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATgcggcg                                                                 <  3:155904/149‑1 (MQ=255)
                                                                                    ggAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGACACGGAGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACAGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATgcggcg                                                                 >  8:201699/1‑149 (MQ=255)
                                                                                     gAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATgcggcg                                                                 >  2:239348/1‑148 (MQ=255)
                                                                                          gcgcTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCTGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATg                                                            <  3:229876/148‑1 (MQ=255)
                                                                                          gcgcTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATgg                                                           >  2:178942/1‑149 (MQ=255)
                                                                                           cgcTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑aa                                                                                        >  3:271518/1‑119 (MQ=255)
                                                                                           cgcTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑aa                                                                                        <  4:271518/119‑1 (MQ=255)
                                                                                           cgcTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGgt                                                          <  3:337682/149‑1 (MQ=255)
                                                                                                cacaTCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCAc                                                     >  8:134130/1‑149 (MQ=255)
                                                                                                 acaTCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCaac                                                    >  5:357296/1‑147 (MQ=255)
                                                                                                 acaTCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCAcc                                                    <  3:75200/149‑1 (MQ=255)
                                                                                                     cgtcgCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTca                                                >  1:69281/1‑149 (MQ=255)
                                                                                                      gtcgCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTca                                                <  2:221523/148‑1 (MQ=255)
                                                                                                      gtcgCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTca                                                >  6:87035/1‑148 (MQ=255)
                                                                                                        cgCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACaa                                             <  8:360375/149‑1 (MQ=255)
                                                                                                         gCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACaa                                             <  6:241968/148‑1 (MQ=255)
                                                                                                          cTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAgg                                                                                                                         <  2:333886/72‑1 (MQ=255)
                                                                                                          cTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAgg                                                                                                                         >  1:333886/1‑72 (MQ=255)
                                                                                                            ttCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAAcc                                           >  2:52213/1‑147 (MQ=255)
                                                                                                                 cTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGGTAAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAggtg                                     <  1:113419/149‑1 (MQ=255)
                                                                                                                    cGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGgtgt                                 <  8:327611/149‑1 (MQ=255)
                                                                                                                    cGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGCGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGgtgt                                 <  1:354871/149‑1 (MQ=255)
                                                                                                                      tGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTgg                               <  2:263616/149‑1 (MQ=255)
                                                                                                                       tgAATGGCTTGCCTCGTCGGGCCTGGGTAATGACCAGGCGCTG‑GATGTCATGCGCGATGGTGACGTGCAGCAGCACCGTGATTGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGgtgt                                 <  4:324977/145‑1 (MQ=255)
                                                                                                                       ggAAGGGCTTGCCACGGTGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTggg                              >  1:261176/1‑149 (MQ=255)
                                                                                                                       ggAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTggg                              >  7:39854/1‑149 (MQ=255)
                                                                                                                       ggAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTggg                              >  3:115478/1‑149 (MQ=255)
                                                                                                                        gAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAg                                        >  7:377441/1‑138 (MQ=255)
                                                                                                                        gAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAg                                        <  8:377441/138‑1 (MQ=255)
                                                                                                                        gAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTggg                              <  3:145439/148‑1 (MQ=255)
                                                                                                                          aGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTGGGTg                            >  8:123130/1‑148 (MQ=255)
                                                                                                                                     cggcggGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTGGGTGCAGCAGCCCtt                 >  3:411659/1‑148 (MQ=255)
                                                                                                                                           gCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTGGGTGCAGCAGCCCTTGAATgag          <  4:115478/149‑1 (MQ=255)
                                                                                                                                            ccTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTGGGTGCAGCAGCCCTTGAATgaga         <  4:54110/149‑1 (MQ=255)
                                                                                                                                                gTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTGGGTGCAGCAGCCCTTGAATGAGATCCt     >  7:323416/1‑149 (MQ=255)
                                                                                                                                                   aaTGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTGGGTGCAGCAGCCCTTGAATGAGATCCTGAt  <  6:372792/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGCAAAGACTGCTGCTGATCGCCTCCCCCCTGATGCGCCGAACCCCGGTCTATGACCGGGCTGCTGCGCTGGCCAAGGCCAAGGGAATGGCGCTGCACATCGTCGCTTTCGACTACGTGGAAGGGCTTGCCACGGCGGGCCTGGTGAATGACCAGGCGCTG‑GCAGTCATGCGCGAGGGTTACGTGCAGCAGCACCGTGAATGGCTGG‑AAACCCAGGCGTTGCCGATGCGGCGCAATGGTGTCACCGTCACAACCGAGGTGGTGTGGGTGCAGCAGCCCTTGAATGAGATCCTGAT  >  NC_002947/3722753‑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: ‑