Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,159,105 G→C 9.8% A361P (GCA→CCA)  opuA → glycine betaine ABC transporter ATP‑binding protein/permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,159,1050GC 9.8% 352.5 / 24.5 114A361P (GCA→CCA) opuAglycine betaine ABC transporter ATP‑binding protein/permease
Reads supporting (aligned to +/- strand):  ref base G (52/50);  new base C (7/4);  total (59/55)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

Only 100 of 120 total aligned reads displayed.

TTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATGGCGCGCCCTGACGCCAGAGGCTGCCGGG  >  NC_002947/3158956‑3159254
                                                                                                                                                     |                                                                                                                                                     
ttCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTgg                                                                                                                                                       <  1:137538/150‑1 (MQ=255)
    ccAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGt                                                                                                                                                   <  1:1408505/150‑1 (MQ=255)
     cAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGtg                                                                                                                                                  >  2:1682281/1‑150 (MQ=255)
      aaaGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGAGAGGGGGCGGGGAGGCTATCGGCCACCGCACTGGCAGtcg                                                                                                                                                 >  2:416668/1‑148 (MQ=255)
            tgtTGGGCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTgg                                                                                                                                           >  1:1016775/1‑150 (MQ=255)
                  gCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCCGCCCTGGCCATTATTGCTGGCGCGGGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATc                                                                                                                                     <  1:900055/150‑1 (MQ=255)
                       gCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGccca                                                                                                                                <  2:1334697/150‑4 (MQ=255)
                       gCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTg                                                                                                                                <  1:27642/150‑1 (MQ=255)
                       gCAGCGCTCGCTGCCGCCGGACGAGGGGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTg                                                                                                                                >  1:1507803/1‑150 (MQ=255)
                           cgctcgctGCCGCCGGACGAGGTGGTTTGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGtt                                                                                                                            <  1:416668/150‑1 (MQ=255)
                            gctcgctGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTtg                                                                                                                           >  1:1751925/1‑150 (MQ=255)
                              tcgctGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTtgtg                                                                                                                         >  1:255152/1‑150 (MQ=255)
                                 ctGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCtt                                                                                                                      <  1:1333782/150‑1 (MQ=255)
                                  tGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCt                                                                                                                                  <  1:179316/137‑1 (MQ=255)
                                  tGCCGCCGGACGAGGTGGGTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATGGGCCACAGCACTGGCAGTGCTGCTGCTCAATCGCTTGGGTTGTGCTTg                                                                                                                     >  1:1048426/1‑150 (MQ=255)
                                   gccgccGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTgg                                                                                                                    <  1:452064/150‑1 (MQ=255)
                                       ccgGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGa                                                                                                                <  2:1582277/150‑1 (MQ=255)
                                        cgGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAg                                                                                                               <  1:1951062/150‑1 (MQ=255)
                                         ggACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGtt                                                                                                                            <  1:1874805/136‑1 (MQ=255)
                                         ggACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGtt                                                                                                                            >  2:1874805/1‑136 (MQ=255)
                                         ggACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGc                                                                                                              <  1:1272142/150‑1 (MQ=255)
                                             gAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGc                                                                                                          >  1:467175/1‑150 (MQ=255)
                                             gAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGATGCTGGTCAAACGCTTGGGTTGGGCTGGGCGGAGCCtac                                                                                                          >  2:139883/1‑148 (MQ=255)
                                                       gttCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGtgctg                                                                                                                                               <  2:1778627/103‑1 (MQ=255)
                                                       gttCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGtgctg                                                                                                                                               >  1:1778627/1‑103 (MQ=255)
                                                        ttCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGtt                                                                                                                            <  2:2131332/121‑1 (MQ=255)
                                                        ttCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGtt                                                                                                                            >  1:2131332/1‑121 (MQ=255)
                                                        ttCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGcacggccacggc                                                                                               <  2:579725/150‑1 (MQ=255)
                                                        ttCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGCCAGtgctgtg                                                                                                                                             <  1:1679653/104‑3 (MQ=255)
                                                        ttCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGCCAGtgctgt                                                                                                                                              <  1:1070557/103‑2 (MQ=255)
                                                        ttCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGCCAGtgctgt                                                                                                                                              >  2:1070557/1‑102 (MQ=255)
                                                           tggctggcGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGt                                                                                            >  2:777207/1‑150 (MQ=255)
                                                             gctgggGCTGCGGGGGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGAACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGa                                                                                          <  2:1048426/150‑1 (MQ=255)
                                                              ctggcGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGgccgtggccgtgca                                                                                         <  1:477971/150‑3 (MQ=255)
                                                                     tgctgGTGGAGGTGGTGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCgg                                                                                  >  2:2403286/1‑150 (MQ=255)
                                                                         ggtggAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGgtgt                                                                              <  2:1507803/150‑1 (MQ=255)
                                                                             gAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGtgctg                                                                                                                                               >  1:706299/1‑81 (MQ=255)
                                                                             gAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGtgctg                                                                                                                                               <  2:706299/81‑1 (MQ=255)
                                                                                gTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGgt                                                                                <  1:1825400/141‑1 (MQ=255)
                                                                                gTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGgt                                                                                >  2:1825400/1‑141 (MQ=255)
                                                                                 tGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCAAGGCCGTGACTCAAGGGGTGTCGCttaag                                                                      >  1:945583/1‑147 (MQ=255)
                                                                                   gcggcCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTg                                                                                                                                >  2:2057396/1‑90 (MQ=255)
                                                                                     ggcCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTggg                                                                                                                   >  2:857944/1‑101 (MQ=255)
                                                                                     ggcCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTggg                                                                                                                   <  1:857944/101‑1 (MQ=255)
                                                                                     ggcCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTAccc                                                                    >  2:359917/1‑148 (MQ=255)
                                                                                       cccTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTgg                                                                                                                    >  2:1712081/1‑98 (MQ=255)
                                                                                       cccTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTgg                                                                                                                    <  1:1712081/98‑1 (MQ=255)
                                                                                         cTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCt                                                                          <  2:1433182/138‑1 (MQ=255)
                                                                                          tGGCCATGCTTGCTGGCGCGTGGGGCGAAGGGGCGGGGTGGCTATCGGCCACAGCACTGAAAACGGCATACGAGATTCTCTTGGGTTGTTCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGgcgc                                                             <  1:2223340/150‑1 (MQ=255)
                                                                                          tGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGgcgc                                                             <  2:1424670/150‑1 (MQ=255)
                                                                                            gCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCgctg                                                           >  1:2049349/1‑150 (MQ=255)
                                                                                             ccATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGt                                                                                                                  <  2:1180111/94‑1 (MQ=255)
                                                                                             ccATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGt                                                                                                                  >  1:1180111/1‑94 (MQ=255)
                                                                                             ccATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCgctgc                                                          >  1:2158144/1‑150 (MQ=255)
                                                                                              cATTCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCgctgct                                                         >  2:1830697/1‑150 (MQ=255)
                                                                                              cATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCgctgct                                                         >  1:98593/1‑150 (MQ=255)
                                                                                                  cttgctGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTg                                                     >  2:1425059/1‑150 (MQ=255)
                                                                                                   ttgctGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTgg                                                    >  2:1763898/1‑150 (MQ=255)
                                                                                                      ctGGCGGGTGGGGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGa                                                 <  2:61638/150‑1 (MQ=255)
                                                                                                      ctGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGt                                                   >  2:752135/1‑148 (MQ=255)
                                                                                                       tggcgcgtggcgcgAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCgctg                                                           <  2:1815012/139‑1 (MQ=255)
                                                                                                       tggcgcgtggcgcgAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAg                                                >  1:1794736/1‑150 (MQ=255)
                                                                                                       tggcgcgtggcgcgAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAg                                                >  1:1972140/1‑150 (MQ=255)
                                                                                                          cgcgtggcgcgAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCa                                             <  1:1996176/150‑1 (MQ=255)
                                                                                                          cgcgtggcgcgAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCa                                             <  1:1425059/150‑1 (MQ=255)
                                                                                                           gcgtggcgcgAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAg                                            >  1:1915426/1‑150 (MQ=255)
                                                                                                             gtggcgcgAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCg                                                                                   <  2:1313729/109‑1 (MQ=255)
                                                                                                                gcgcgAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACt                                       >  1:99400/1‑150 (MQ=255)
                                                                                                                   cgAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTCAc                                    <  1:2403286/150‑1 (MQ=255)
                                                                                                                     aGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCg                                  >  1:930433/1‑150 (MQ=255)
                                                                                                                        gggCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGa                                     >  2:2407876/1‑144 (MQ=255)
                                                                                                                        gggCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGa                                     <  1:2407876/144‑1 (MQ=255)
                                                                                                                        gggCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTAt                               <  2:2340845/150‑1 (MQ=255)
                                                                                                                         ggCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATg                              <  1:39348/150‑1 (MQ=255)
                                                                                                                          gCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATgg                             <  2:2158144/150‑1 (MQ=255)
                                                                                                                               gTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATGgcgcgc                        >  1:2363267/1‑150 (MQ=255)
                                                                                                                                tGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATGGCGCGcc                       >  2:929286/1‑150 (MQ=255)
                                                                                                                                tGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATGGCGCGcc                       <  2:1016775/150‑1 (MQ=255)
                                                                                                                                    tATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTAcc                                                                     >  2:1860426/1‑100 (MQ=255)
                                                                                                                                    tATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTAcc                                                                     <  1:1860426/100‑1 (MQ=255)
                                                                                                                                    tATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGCGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTAccc                                                                    <  1:736721/101‑1 (MQ=255)
                                                                                                                                     aTCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATGGCGCGCCCTGAc                  <  1:1316071/150‑1 (MQ=255)
                                                                                                                                         gCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGgt                                                                                >  2:667321/1‑84 (MQ=255)
                                                                                                                                         gCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGgt                                                                                <  1:667321/84‑1 (MQ=255)
                                                                                                                                         gCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATGGCGCGCCCTGACGCCa              >  1:586473/1‑150 (MQ=255)
                                                                                                                                         ccagcaGCACTGCCAGTGCTGTTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGgcgc                                                             <  1:1672733/99‑1 (MQ=255)
                                                                                                                                         ccagcaGCACTGCCAGTGCTGTTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGgcgc                                                             >  2:1672733/5‑103 (MQ=255)
                                                                                                                                         ccagcaGCACTGCCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAgg                                           <  2:1504787/117‑1 (MQ=255)
                                                                                                                                         ccagcaGCACTGCCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAgg                                           >  1:1504787/5‑121 (MQ=255)
                                                                                                                                          ccACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGcacggcca                                                                                                   >  1:2107707/1‑64 (MQ=255)
                                                                                                                                          ccACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGcacggcca                                                                                                   <  2:2107707/64‑1 (MQ=255)
                                                                                                                                          cagcaGCACTGCCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGGGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATGGCGCGCCCTGACGCCag             >  1:1524875/4‑150 (MQ=255)
                                                                                                                                            acaGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATGGCGCGCCCTGACGCCagag           >  2:890551/1‑150 (MQ=255)
                                                                                                                                            acaGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATGGCGCGCCCTGACGCCagag           <  2:255152/150‑1 (MQ=255)
                                                                                                                                             caGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATGGCGCGCCCTGACGCCAGAgg          <  2:1751925/150‑1 (MQ=255)
                                                                                                                                                    ggCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATGGCGCGCCCTGACGCCAGAGGCTGCCgg   <  1:1830697/150‑1 (MQ=255)
                                                                                                                                                     gCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTgg                                                                <  2:2168119/88‑1 (MQ=255)
                                                                                                                                                     gCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTgg                                                                >  1:2168119/1‑88 (MQ=255)
                                                                                                                                                     gCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATGGCGCGCCCTGACGCCAGAGGCTGCCggg  >  1:2008950/1‑150 (MQ=255)
                                                                                                                                                     gCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATGGCGCGCCCTGACGCCAGAGGCTGCCggg  >  1:573358/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCGCCAAAGCGTGTTGGGCCCGGCAGCGCTCGCTGCCGCCGGACGAGGTGGTTGGTTCTGGCTGGCGCTGCTGGTGGAGGTGGCGGCCCTGGCCATGCTTGCTGGCGCGTGGCGCGAGGGGGCGGGGTGGCTATCGGCCACAGCACTGGCAGTGCTGCTGGTCAATCGCTTGGGTTGTGCTTGGGTGAGCCTGCACGGCCACGGCCGTGACTCACCGGGTGTCGCTTTACCCTTGGCGCTGCTACTGGTGAGTCAGGCACTGACCGTATGGCGCGCCCTGACGCCAGAGGCTGCCGGG  >  NC_002947/3158956‑3159254

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