Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,143,487 C→G 5.1% S89R (AGC→AGG)  PP_1899 → ExbD/TolR family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,143,4870CG 5.1% 256.7 / 3.3 78S89R (AGC→AGG) PP_1899ExbD/TolR family protein
Reads supporting (aligned to +/- strand):  ref base C (34/40);  new base G (2/2);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

CGCGTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAG  >  NC_002947/2143347‑2143634
                                                                                                                                            |                                                                                                                                                    
cgcgTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGcc                                                                                                                                                <  1:191507/147‑1 (MQ=255)
 gcgTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCAt                                                                                                                                              <  5:109143/148‑1 (MQ=255)
 gcgTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATa                                                                                                                                             <  1:60436/149‑1 (MQ=255)
      gAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAgc                                                                                                                                         <  5:26125/148‑1 (MQ=255)
       aGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAgcgc                                                                                                                                       <  1:156854/149‑1 (MQ=255)
         cTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAgcgc                                                                                                                                       >  2:216625/1‑147 (MQ=255)
         cTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGa                                                                                                                                     <  3:197357/149‑1 (MQ=255)
          tGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCTCGCTCAGGGTGGcc                                                                                                                                                >  8:147726/1‑137 (MQ=255)
          tGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCTCGCTCAGGGTGGcc                                                                                                                                                <  7:147726/137‑1 (MQ=255)
           gCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGt                                                                                                                                   >  3:104836/1‑149 (MQ=255)
            cccGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTc                                                                                                                                  >  8:356946/1‑149 (MQ=255)
             ccGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCg                                                                                                                                 >  6:93633/1‑149 (MQ=255)
               gAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCggg                                                                                                                               <  2:90227/149‑1 (MQ=255)
                 aGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGggcg                                                                                                                             <  4:242554/149‑1 (MQ=255)
                   ccaccaGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGggcggc                                                                                                                           <  7:6450/149‑1 (MQ=255)
                      ccaGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑Gac                                                                                                                       <  6:381585/149‑1 (MQ=255)
                      ccaGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑Ga                                                                                                                        >  6:328486/1‑148 (MQ=255)
                      ccaGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑Ga                                                                                                                        >  6:334207/1‑148 (MQ=255)
                       caGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCggg                                                                                                                               >  5:161740/1‑141 (MQ=255)
                       caGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCggg                                                                                                                               <  6:161740/141‑1 (MQ=255)
                         gcgcCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑Gacaa                                                                                                                     >  5:333221/1‑148 (MQ=255)
                           gcCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGGGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑Gacaac                                                                                                                    >  1:336324/1‑147 (MQ=255)
                              gAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGtt                                                                                                               <  8:26224/149‑1 (MQ=255)
                                ctgtgcgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GAcaacaa                                                                                                                  >  4:206454/5‑144 (MQ=255)
                                   ggcgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGc                                                                                                          >  2:315860/1‑149 (MQ=255)
                                   ggcgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGc                                                                                                          >  4:130395/1‑149 (MQ=255)
                                   ggcgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGc                                                                                                          <  4:104836/149‑1 (MQ=255)
                                    gcgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCt                                                                                                         >  6:216614/1‑149 (MQ=255)
                                     cgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTg                                                                                                        <  3:80277/149‑1 (MQ=255)
                                     cgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTg                                                                                                        <  5:93633/149‑1 (MQ=255)
                                     cgccggcCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTg                                                                                                        <  1:315860/149‑1 (MQ=255)
                                      gccggccgACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCGTGAGCGAGGCCATTGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTgg                                                                                                       <  2:178944/149‑1 (MQ=255)
                                           ccgACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATc                                                                                                  <  1:216625/149‑1 (MQ=255)
                                            cgACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCa                                                                                                 <  6:398347/149‑1 (MQ=255)
                                                ccGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCAGGCCATCAgcgc                                                                                             <  7:356946/149‑1 (MQ=255)
                                                   ggCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGa                                                                                          <  2:336324/149‑1 (MQ=255)
                                                       aGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGc                                                                                      <  3:130395/149‑1 (MQ=255)
                                                       aGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGc                                                                                      <  5:216614/149‑1 (MQ=255)
                                                          tggtggAAATCACCTTCAGCGCCGCAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAg                                                                                   <  7:76099/149‑1 (MQ=255)
                                                              ggAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGAccc                                                                               >  1:357301/1‑149 (MQ=255)
                                                               gAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAgcgc                                                                                                                                       <  1:195909/93‑1 (MQ=255)
                                                               gAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAgcgc                                                                                                                                       >  2:195909/1‑93 (MQ=255)
                                                                 aaTCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGc                                                                            <  6:333221/149‑1 (MQ=255)
                                                                      ccATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAgcg                                                                                              >  8:219072/1‑126 (MQ=255)
                                                                      ccATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAgcg                                                                                              <  7:219072/126‑1 (MQ=255)
                                                                      ccATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAg                                                                       >  2:213117/1‑149 (MQ=255)
                                                                       cATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGc                                                                                                                                                 <  6:228272/75‑1 (MQ=255)
                                                                       cATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGc                                                                                                                                                 >  5:228272/1‑75 (MQ=255)
                                                                       cATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAgg                                                                      >  2:334457/1‑149 (MQ=255)
                                                                        aTCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGCAGACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAgg                                                                      <  1:101671/149‑1 (MQ=255)
                                                                        aTCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGc                                                                     >  4:186232/1‑149 (MQ=255)
                                                                               ccGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGt                                                               >  2:315051/1‑148 (MQ=255)
                                                                               ccGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTc                                                              <  1:213117/149‑1 (MQ=255)
                                                                                       gtgtATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAAt                                                      >  6:162473/1‑149 (MQ=255)
                                                                                                gTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCAGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGAtgctgc                                             >  2:129030/1‑149 (MQ=255)
                                                                                                  gAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCg                                           >  6:187935/1‑149 (MQ=255)
                                                                                                  gAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCg                                           >  6:194844/1‑149 (MQ=255)
                                                                                                      aacCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGc                                         >  7:162558/1‑147 (MQ=255)
                                                                                                            ctgctgCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCgg                                 >  5:15706/1‑149 (MQ=255)
                                                                                                               ctgcCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGttt                              >  7:316022/1‑149 (MQ=255)
                                                                                                                  cccAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCa                            <  6:15706/148‑1 (MQ=255)
                                                                                                                  cccAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCa                            <  5:194844/148‑1 (MQ=255)
                                                                                                                   ccAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGc                                                                                                      <  3:97459/73‑1 (MQ=255)
                                                                                                                   ccAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGc                                                                                                      >  4:97459/1‑73 (MQ=255)
                                                                                                                   ccAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGc                          >  7:34785/1‑149 (MQ=255)
                                                                                                                    cAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGa                                                   <  7:351045/123‑1 (MQ=255)
                                                                                                                    cAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGa                                                   >  8:351045/1‑123 (MQ=255)
                                                                                                                    cAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGc                          >  5:93459/1‑148 (MQ=255)
                                                                                                                     aaaaGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCa                        <  6:93459/149‑1 (MQ=255)
                                                                                                                           tctatGGCCTCGCTCAGGGTGGCCATAGAGCGCGAGTCGGGCGGC‑GGCAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCa                            >  5:114946/5‑139 (MQ=255)
                                                                                                                           tctatGGCCTCGCTCAGGGTGGCCATAGAGCGCGAGTCGGGCGGC‑GGCAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCa                            <  6:114946/135‑1 (MQ=255)
                                                                                                                            aCCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTgcgc                 <  1:334457/149‑1 (MQ=255)
                                                                                                                            aCCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTgcgc                 <  5:187935/149‑1 (MQ=255)
                                                                                                                                 gCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGa             <  1:315051/148‑1 (MQ=255)
                                                                                                                                      ccTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCg       >  3:17842/1‑149 (MQ=255)
                                                                                                                                      ccTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCg       >  1:383841/1‑149 (MQ=255)
                                                                                                                                           gCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAg  <  3:186232/149‑1 (MQ=255)
                                                                                                                                           gCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAg  <  2:383841/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
CGCGTCGAGCTGCCCGAAGCCACCAGCGCCGAGCCGGCGCCGGCCGACCCGGGCAAGCTGGTGGAAATCACCATCAGCGCCGAAGGCGTGTATTCGGTGAACAACCACCTGCTGCCCAAAAGCGACCTGGCCACCCTGAGCGAGGCCATAGAGCGCGAGTCGGGCGGC‑GACAACAAGTTGCCGCTGGCCATCAGCGCCGACGGCAAGACCCCGCATCAGGCTGTGGTCACCGCAATGGATGCTGCCGGCAAGCTCGGTTTCAGCCAGTTGCGCATGACCACCGTCGAG  >  NC_002947/2143347‑2143634

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