Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,489,780 A→C 6.1% M88L (ATG→CTG) ‡ PP_4825 → MarC family membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,489,7800AC 6.1% 249.8 / 6.9 82M88L (ATG→CTG) ‡PP_4825MarC family membrane protein
Reads supporting (aligned to +/- strand):  ref base A (38/37);  new base C (3/2);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

GCAACTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGAAA  >  NC_002947/5489632‑5489924
                                                                                                                                                       |                                                                                                                                                
gcaaCTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCa                                                                                                                                                  <  1:145153/149‑1 (MQ=255)
gcaaCTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCa                                                                                                                                                  >  5:41905/1‑149 (MQ=255)
gcaaCTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCa                                                                                                                                                  >  7:277311/1‑149 (MQ=255)
 caaCTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCAt                                                                                                                                                 <  3:325808/149‑1 (MQ=255)
    cTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGc                                                                                                                                              >  1:4796/1‑149 (MQ=255)
     tGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGcc                                                                                                                                             >  2:55919/1‑149 (MQ=255)
     tGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGcc                                                                                                                                             >  2:280964/1‑149 (MQ=255)
       gCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCa                                                                                                                                           <  6:190360/149‑1 (MQ=255)
       gCCTCGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGcc                                                                                                                                             >  7:142535/1‑147 (MQ=255)
         cTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAgg                                                                                                                                         <  5:60974/149‑1 (MQ=255)
         cTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAgg                                                                                                                                         >  8:313821/1‑149 (MQ=255)
          tGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCATCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAggg                                                                                                                                        >  7:306733/1‑149 (MQ=255)
            gAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCCCTCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAgg                                                                                                                                         <  1:110791/148‑1 (MQ=255)
            gAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCATGCCCAggg                                                                                                                                        <  4:167398/147‑1 (MQ=255)
            gAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAggg                                                                                                                                        <  8:255867/147‑1 (MQ=255)
            gAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAggg                                                                                                                                        <  5:3312/147‑1 (MQ=255)
            gAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGGTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGc                                                                                                                                       <  4:239056/148‑1 (MQ=255)
              aGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAg                                                                                                                                    >  8:352699/1‑149 (MQ=255)
                 tCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCcg                                                                                                                                 >  7:339932/1‑149 (MQ=255)
                  cgcgTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCcgg                                                                                                                                >  1:135891/1‑149 (MQ=255)
                    cgTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCcggcg                                                                                                                              <  8:142535/149‑1 (MQ=255)
                                 ggTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGa                                                                                                                   <  5:360703/147‑1 (MQ=255)
                                  gTCGCCAGCGTGTTGGTGAACGTTTTGGGGGGGGCTTTCTTCGTATTTTTCGGCATCC‑‑CCGCGGTTGCTTTCCGCATGGTTGCTGGCAGCGTGGTGTTCATTTCGCCCGTGGGAAGGGCCCTGGGCAAGCCGGCGGTACAGGCAGACa                                                                                                                 <  2:135891/148‑1 (MQ=255)
                                     gCCTGCGTGTTGCTGTACCTGTTCGGGCGGGCTATTTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCTGGGCAAGCCGGCGGTACAGGCAGACAACGt                                                                                                             <  6:138639/149‑1 (MQ=255)
                                       cAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGtt                                                                                                            >  5:299862/1‑148 (MQ=255)
                                       cAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTc                                                                                                           <  7:69713/149‑1 (MQ=255)
                                           gtgtTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTcagca                                                                                                       <  3:366550/149‑1 (MQ=255)
                                            tgtTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAggg                                                                                                                                        >  8:339782/1‑115 (MQ=255)
                                            tgtTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAggg                                                                                                                                        <  7:339782/115‑1 (MQ=255)
                                            tgtTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTcagcag                                                                                                      <  7:313821/149‑1 (MQ=255)
                                                 cTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTc                                                                                                           <  4:348420/139‑1 (MQ=255)
                                                 cTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTc                                                                                                           >  3:348420/1‑139 (MQ=255)
                                                    tACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGGCTTCCGCATCGGTGGTGGCAGCGTGCTGTTCATTTCGGCCCTGCGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGAtgtg                                                                                                <  1:280964/148‑1 (MQ=255)
                                                     aCCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGAc                                                                                              >  4:300265/1‑148 (MQ=255)
                                                       cTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCAt                                                                                           >  1:313302/1‑149 (MQ=255)
                                                            cgggcgggCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTAc                                                                                                                          <  8:94286/113‑1 (MQ=255)
                                                            cgggcgggCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTAc                                                                                                                          >  7:94286/1‑113 (MQ=255)
                                                            cgggcgggCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGGGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTa                                                                                       <  3:177827/148‑1 (MQ=255)
                                                             gggcgggcGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTAcc                                                                                     >  3:236137/1‑149 (MQ=255)
                                                              ggcgggcgATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCg                                                                                    >  4:1309/1‑149 (MQ=255)
                                                              ggcgggcgATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCg                                                                                    >  4:169465/1‑149 (MQ=255)
                                                                 gggcgATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTg                                                                                 >  4:31260/1‑149 (MQ=255)
                                                                   gcgATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGAc                                                                               <  7:118658/149‑1 (MQ=255)
                                                                   gcgATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGAc                                                                               <  1:55919/149‑1 (MQ=255)
                                                                     gATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGcc                                                                                                                                  <  1:33081/96‑1 (MQ=255)
                                                                     gATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGcc                                                                                                                                  >  2:33081/1‑96 (MQ=255)
                                                                             gAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCATGCCCAggg                                                                                                                                        >  2:323773/1‑82 (MQ=255)
                                                                             gAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCATGCCCAggg                                                                                                                                        <  1:323773/82‑1 (MQ=255)
                                                                             gAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCt                                                                      >  1:340600/1‑148 (MQ=255)
                                                                             gAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTg                                                                     >  7:260448/1‑149 (MQ=255)
                                                                                         tCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGcccc                                                         >  5:187932/1‑149 (MQ=255)
                                                                                           a‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCg                                                        <  2:4796/148‑1 (MQ=255)
                                                                                                   gATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCAt                                                 <  8:339932/149‑1 (MQ=255)
                                                                                                     tGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCg                                               <  8:260448/149‑1 (MQ=255)
                                                                                                        ggaaggCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTg                                            >  5:43983/6‑149 (MQ=255)
                                                                                                         ttCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGc                                           >  2:218543/1‑149 (MQ=255)
                                                                                                             gCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCctg                                       >  1:324739/1‑149 (MQ=255)
                                                                                                             gCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCctg                                       >  1:237979/1‑149 (MQ=255)
                                                                                                             ccATGGGGGCGGGCTGCGTGCTTTCCTTTTGGGCCTTGGGCTTGGCCCTGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCGATGGTGccggg                                       <  2:340600/148‑4 (MQ=255)
                                                                                                                 cGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGTTGCCCTGCTgg                                   >  3:337771/1‑149 (MQ=255)
                                                                                                                  ggTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTgg                                   <  4:236137/148‑1 (MQ=255)
                                                                                                                       tGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATgg                             <  7:352699/149‑1 (MQ=255)
                                                                                                                         gCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATgggg                           >  3:31966/1‑149 (MQ=255)
                                                                                                                          cAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATggggg                          <  6:43983/149‑1 (MQ=255)
                                                                                                                           aGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGt                         >  6:217404/1‑149 (MQ=255)
                                                                                                                           aGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGt                         <  3:169465/149‑1 (MQ=255)
                                                                                                                             ctggCTGTTCATTTCGGCCCTGGGATGGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGa                       <  6:41905/146‑1 (MQ=255)
                                                                                                                                    ttCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTcagc                                                                                                        >  1:153637/1‑61 (MQ=255)
                                                                                                                                    ttCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTcagc                                                                                                        <  2:153637/61‑1 (MQ=255)
                                                                                                                                     tCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCg               >  3:103904/1‑149 (MQ=255)
                                                                                                                                       aTTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCa             >  2:329310/1‑149 (MQ=255)
                                                                                                                                       aTTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCa             >  3:314555/1‑149 (MQ=255)
                                                                                                                                       aTTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCa             <  4:103904/149‑1 (MQ=255)
                                                                                                                                        tttCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCAc            >  7:20746/1‑149 (MQ=255)
                                                                                                                                        tttCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCAc            <  5:217404/149‑1 (MQ=255)
                                                                                                                                         gcttgCCCTGGGCCATGCCCAGGGCCAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGCTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACt           >  4:29699/5‑149 (MQ=255)
                                                                                                                                          tCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACt           <  7:212763/148‑1 (MQ=255)
                                                                                                                                            cgtccTGGGCAGGGCCGAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGAAGGATGTTTCCTTCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTgg         <  8:306733/145‑1 (MQ=255)
                                                                                                                                             gCCCTGGGCCATGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGa       >  2:131551/1‑149 (MQ=255)
                                                                                                                                             gCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGa       >  2:321932/1‑149 (MQ=255)
                                                                                                                                              cccTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGAt      >  7:76230/1‑149 (MQ=255)
                                                                                                                                               ccTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATg                              <  2:346936/124‑1 (MQ=255)
                                                                                                                                               ccTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATg                              >  1:346936/1‑124 (MQ=255)
                                                                                                                                                  gggCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGaaa  <  6:299862/149‑1 (MQ=255)
                                                                                                                                                    ggCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGaa   <  6:187932/147‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                
GCAACTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGAAA  >  NC_002947/5489632‑5489924

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