Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,489,782 G→T 6.2% M88I (ATG→ATT) ‡ 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,7820GT 6.2% 263.8 / 6.2 81M88I (ATG→ATT) ‡PP_4825MarC family membrane protein
Reads supporting (aligned to +/- strand):  ref base G (38/38);  new base T (3/2);  total (41/40)
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.92e-01

CTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGAAA  >  NC_002947/5489636‑5489924
                                                                                                                                                     |                                                                                                                                              
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                                                                                                                                         >  8:313821/1‑149 (MQ=255)
     cTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAgg                                                                                                                                         <  5:60974/149‑1 (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                                                                                                                                        <  7:339782/115‑1 (MQ=255)
                                        tgtTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAggg                                                                                                                                        >  8:339782/1‑115 (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                                                                                                                          >  7:94286/1‑113 (MQ=255)
                                                        cgggcgggCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTAc                                                                                                                          <  8:94286/113‑1 (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                                                                                                                                        <  1:323773/82‑1 (MQ=255)
                                                                         gAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCCATGCCCAggg                                                                                                                                        >  2:323773/1‑82 (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:237979/1‑149 (MQ=255)
                                                                                                         gCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCctg                                       >  1:324739/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                         <  3:169465/149‑1 (MQ=255)
                                                                                                                       aGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGt                         >  6:217404/1‑149 (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             <  4:103904/149‑1 (MQ=255)
                                                                                                                                   aTTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCa             >  3:314555/1‑149 (MQ=255)
                                                                                                                                   aTTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCa             >  2:329310/1‑149 (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                              >  1:346936/1‑124 (MQ=255)
                                                                                                                                           ccTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATg                              <  2:346936/124‑1 (MQ=255)
                                                                                                                                              gggCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGaaa  <  6:299862/149‑1 (MQ=255)
                                                                                                                                                ggCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGaa   <  6:187932/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
CTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCA‑‑CCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTCGGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGAAA  >  NC_002947/5489636‑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: