Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,062,746 C→G 5.8% P554R (CCG→CGG)  PP_1838 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,062,7460CG 5.8% 267.0 / 4.2 88P554R (CCG→CGG) PP_1838membrane protein
Reads supporting (aligned to +/- strand):  ref base C (42/40);  new base G (3/2);  total (45/43)
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.91e-01

CAGTTCTTCGAGAAGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGG  >  NC_002947/2062604‑2062892
                                                                                                                                              |                                                                                                                                                  
cAGTTCTTCGAGAAGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCAt                                                                                                                                              >  6:323041/1‑149 (MQ=255)
     cttcGAGAAGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAg                                                                                                                                          <  1:141001/148‑1 (MQ=255)
      ttcGAGAAGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAg                                                                                                                                          >  6:14590/1‑147 (MQ=255)
        cGAGAAGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGgag                                                                                                                                       >  5:23087/1‑148 (MQ=255)
         gagaAGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGgag                                                                                                                                       >  8:106314/1‑147 (MQ=255)
           gaAGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGt                                                                                                                                   >  4:188961/1‑149 (MQ=255)
           gaAGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGACATCGTCGGTGACCATGGCTATGGCAACCACCAGCAGGTTACAGAGCAGGACCTGGGCAGATTCAACGTGACGCTGCTGCAGATTGCCCCGTGCATTAAAGAGAAg                                                                                                                                    >  3:341654/1‑148 (MQ=255)
            aaGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGga                                                                                                                                        <  3:32426/143‑1 (MQ=255)
            aaGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGga                                                                                                                                        >  4:32426/1‑143 (MQ=255)
            aaGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGtt                                                                                                                                  <  5:14590/149‑1 (MQ=255)
             aGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTc                                                                                                                                 <  2:257591/149‑1 (MQ=255)
                cgcgTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCTGGGCATTCAGGAGAAGTTcggt                                                                                                                              <  7:55714/149‑1 (MQ=255)
                 gcgTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTcggtg                                                                                                                             >  6:26255/1‑149 (MQ=255)
                    tAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTcggtgcgg                                                                                                                          <  5:26255/149‑1 (MQ=255)
                     aaGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTcggtgcggt                                                                                                                         >  3:42567/1‑149 (MQ=255)
                       ggAGCATTACTTCAAGGACACCCTGTTCGTGATCGTCGGTGTCCATGGCTTTGGCAACCTCCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTGAACGGGCCGCTGGTGCTGATTGCCCCGGGCATTCTGGAGAAGTTCGGTGCGGTGa                                                                                                                       <  5:314601/149‑1 (MQ=255)
                        gAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGaa                                                                                                                      <  5:346126/149‑1 (MQ=255)
                                cTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCggg                                                                                                                                                 >  6:236284/1‑114 (MQ=255)
                                cTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCggg                                                                                                                                                 <  5:236284/114‑1 (MQ=255)
                                 ttCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGt                                                                                                             <  6:23087/149‑1 (MQ=255)
                                 ttCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGt                                                                                                             >  6:161039/1‑149 (MQ=255)
                                  tCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCGCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTg                                                                                                            >  5:327299/1‑149 (MQ=255)
                                  tCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGga                                                                                                                                        <  8:7164/121‑1 (MQ=255)
                                  tCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGga                                                                                                                                        <  8:1560/121‑1 (MQ=255)
                                  tCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTg                                                                                                            >  5:49853/1‑149 (MQ=255)
                                   cAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTgg                                                                                                           >  3:292936/1‑149 (MQ=255)
                                    aaGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTggg                                                                                                          >  2:87090/1‑149 (MQ=255)
                                    aaGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTgg                                                                                                           <  4:323415/148‑1 (MQ=255)
                                      ggACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTa                                                                                                        <  7:66506/149‑1 (MQ=255)
                                         cacCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGt                                                                                                             >  4:224475/1‑141 (MQ=255)
                                         cacCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGt                                                                                                             <  3:224475/141‑1 (MQ=255)
                                          acCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTAcc                                                                                                      >  8:257387/1‑147 (MQ=255)
                                           cccTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAg                                                                                                   <  4:169162/149‑1 (MQ=255)
                                            ccTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAgg                                                                                                  >  1:174672/1‑149 (MQ=255)
                                                ttCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCTAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTc                                                                                                <  5:161039/147‑1 (MQ=255)
                                                ttCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGa                                                                                              <  7:63834/149‑1 (MQ=255)
                                                       tcgtcgGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCAGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTg                                                                                       >  1:269150/1‑149 (MQ=255)
                                                         gtcgGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGc                                                                                      <  4:292936/148‑1 (MQ=255)
                                                         gtcgGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGc                                                                                      <  7:257387/148‑1 (MQ=255)
                                                          tcgGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCg                                                                                    >  5:376222/1‑149 (MQ=255)
                                                             gTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTgg                                                                                                           >  2:286538/1‑123 (MQ=255)
                                                             gTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTgg                                                                                                           <  1:286538/123‑1 (MQ=255)
                                                             gTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGAc                                                                                  >  6:192358/1‑148 (MQ=255)
                                                              tGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGACGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACca                                                                                <  5:292073/149‑1 (MQ=255)
                                                                  cATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACcatcat                                                                            >  2:83000/1‑149 (MQ=255)
                                                                  cATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACcatcat                                                                            <  5:323041/149‑1 (MQ=255)
                                                                  cATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGAGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACcatcat                                                                            >  3:290586/1‑149 (MQ=255)
                                                                    tGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATgg                                                                          <  1:83000/149‑1 (MQ=255)
                                                                            gCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGgccgccag                                                                  >  5:236324/1‑147 (MQ=255)
                                                                             cAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTcg                                                                                                                                <  8:80050/86‑1 (MQ=255)
                                                                             cAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTcg                                                                                                                                >  7:80050/1‑86 (MQ=255)
                                                                             cAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTgg                                                                 <  7:106314/149‑1 (MQ=255)
                                                                              aaCCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGcc                                                                                     >  2:345144/1‑128 (MQ=255)
                                                                              aaCCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGcc                                                                                     <  1:345144/128‑1 (MQ=255)
                                                                              aaCCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTggg                                                                >  8:349779/1‑149 (MQ=255)
                                                                               accaccAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGt                                                               <  2:174672/149‑1 (MQ=255)
                                                                                 caccaGCAGGTTACCGAGCTGGACCTGTTCCGCTTCAACGTGCAGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGACGCCTGGGTgg                                                             >  2:317568/1‑149 (MQ=255)
                                                                                  accaGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGACGCTGCTGCTGATTGCCCGGGGCATTCAGGAGAAGTTCGGTGCTGTGAATCACACCGTGGGTACCCAGGTCGAGATCGTGCCGACCATCATGGGCCGCCTGGGTGGc                                                            <  4:341654/149‑1 (MQ=255)
                                                                                      gcagGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAgg                                                        <  3:188961/149‑1 (MQ=255)
                                                                                         ggTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTc                                                     <  2:215729/149‑1 (MQ=255)
                                                                                          gTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTg                                                              >  5:21080/1‑139 (MQ=255)
                                                                                          gTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTg                                                              <  6:21080/139‑1 (MQ=255)
                                                                                            tACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTc                                                     >  5:283928/1‑146 (MQ=255)
                                                                                            tACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTc                                                     <  6:283928/146‑1 (MQ=255)
                                                                                             aCCGAGCTTGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGCATCACACCGTGGGTACCCAGGTCGACAGCGTGCCGACCATGATGGGCCGCCTGGGTGGCAAGGCACGcca                                                 >  7:249633/1‑149 (MQ=255)
                                                                                              ccGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGccac                                                <  1:87090/149‑1 (MQ=255)
                                                                                               cGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTg                                                              <  5:2203/134‑1 (MQ=255)
                                                                                               cGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTg                                                              >  6:2203/1‑134 (MQ=255)
                                                                                               cGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGccac                                                >  6:176423/1‑148 (MQ=255)
                                                                                                           gggCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGgc                                   >  1:42369/1‑149 (MQ=255)
                                                                                                           gggCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGgc                                   >  5:66430/1‑149 (MQ=255)
                                                                                                               cGCTTCAACGTGCCGCTGCTGCTGATTGCCCCTGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTTGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCACCACCAGTGCTGGGGGCGCGa                               <  4:42567/149‑1 (MQ=255)
                                                                                                                 cTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAgg                                                                                                                                         >  2:12275/1‑41 (MQ=255)
                                                                                                                 cTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAgg                                                                                                                                         <  1:12275/41‑1 (MQ=255)
                                                                                                                 cTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACcatcat                                                                            >  6:12977/1‑102 (MQ=255)
                                                                                                                 cTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACcatcat                                                                            <  5:12977/102‑1 (MQ=255)
                                                                                                                   tCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTg                           >  1:129749/1‑149 (MQ=255)
                                                                                                                     aaCGTGCCGCTTCTGCTTATTGCCCCGTGCATTCAGTAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCCTCATGCGCCGCCTGGGTGGCAAGGCTCGCCACCAGTTCTGGGGGCGCGACCTGc                          <  1:317568/148‑1 (MQ=255)
                                                                                                                      aCGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCt                         >  2:138485/1‑148 (MQ=255)
                                                                                                                                 ctgctgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGa             <  6:66430/149‑1 (MQ=255)
                                                                                                                                   gctgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGAc                              >  7:197352/1‑130 (MQ=255)
                                                                                                                                   gctgATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGAc                              <  8:197352/130‑1 (MQ=255)
                                                                                                                                    ctgAATGCCCGGGGCAATCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAggg          >  5:178458/1‑149 (MQ=255)
                                                                                                                                    ctgAATGCCCGGGGCAATCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAgg           <  5:203922/148‑1 (MQ=255)
                                                                                                                                    ctgAATGCCCGGGGCAATCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAgg           >  6:203922/1‑148 (MQ=255)
                                                                                                                                    ctgAATGCCCGGGGCAAACAGGAGAAGTTCGGTGCGGTCAATCACACCGTCGGTACCCAGCTCCACATCGTGCCGACCATCATCGGCCGCCTGGGTGGCAAGGCTCGCCCCCAGTCCTGGGGGCGCGCCCTGCACAACCTGCCCGAggg          >  6:197205/1‑149 (MQ=255)
                                                                                                                                       aTTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGTCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAggg          <  6:36749/146‑1 (MQ=255)
                                                                                                                                       aTTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGa       >  8:268544/1‑149 (MQ=255)
                                                                                                                                            cccGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAgg  <  2:42369/149‑1 (MQ=255)
                                                                                                                                            cccGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAgg  >  2:290452/1‑149 (MQ=255)
                                                                                                                                            cccGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAgg  >  2:112347/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CAGTTCTTCGAGAAGGCGCGTAAGGAGCCTTACTTCAAGGACACCCTGTTCGTCATCGTCGGTGACCATGGCTTTGGCAACCACCAGCAGGTTACCGAGCTGGACCTGGGCCGCTTCAACGTGCCGCTGCTGCTGATTGCCCCGGGCATTCAGGAGAAGTTCGGTGCGGTGAATCACACCGTGGGTACCCAGGTCGACATCGTGCCGACCATCATGGGCCGCCTGGGTGGCAAGGCTCGCCACCAGTGCTGGGGGCGCGACCTGCTCAACCTGCCCGAGGGTGACCAGG  >  NC_002947/2062604‑2062892

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: ‑