Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,043,276 C→A 6.4% S361R (AGC→AGA colS → integral membrane sensor signal transduction histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,043,2760CA 6.4% 200.5 / 3.1 63S361R (AGC→AGAcolSintegral membrane sensor signal transduction histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (31/27);  new base A (2/2);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

GATGGCCGCGTGAGCGCTAGCCCTGTGCTGTACAACGCCACCTTCCTGCAGTCGGTGATGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCTGGTGCAGCGGATCTGCGACGACCAAGGTTGGCGTGTCACCCTGACA  >  NC_002947/1043127‑1043420
                                                                                                                                                     |                                                                                                                                                
gATGGCCGCGTGAGCGCTAGCCCTGTGCTGTACAACGCCACCTTCCTGCAGTCGGTGATGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAgc                                                                                                                                                  <  1:660098/150‑1 (MQ=255)
      cgcgTGAGCGCTAGCCCTGTGCTGTACAACGCCACCTTCCTGCAGTCGGTGATGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTg                                                                                                                                            <  2:1292974/150‑1 (MQ=255)
       gcgTGAGCGCTAGCCCTGTGCTGTACAACGCCACCTTCCTGCAGTCGGTGATGGGCAACCTGCTGCTCAATGCGGCGAACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTgg                                                                                                                                           <  1:100595/150‑1 (MQ=255)
       gcgTGAGCGCTAGCCCTGTGCTGTACAACGCCACCTTCCTGCAGTCGGTGATGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTgg                                                                                                                                           <  2:996306/150‑1 (MQ=255)
                  aGCCCTGTGCTGTACAACGCCACCTTCCTGCAGTCGGTGATGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGaa                                                                                                                                <  1:959458/150‑1 (MQ=255)
                     ccTGTGCTGTACAACGCCACCTTCCTGCAGTCGGTGATGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAagcg                                                                                                                             <  2:646996/150‑3 (MQ=255)
                     ccTGTGCTGTACAACGCCACCTTCCTGCAGTCGGTGATGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAagag                                                                                                                             >  1:1329448/1‑150 (MQ=255)
                         tgCTGTACAACGCCACCTTCCTGCAGTCGGTGATGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGagcagc                                                                                                                         >  1:783689/1‑150 (MQ=255)
                             gTACAACGCCACCTTCCTGCAGTCGGTGATGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGt                                                                                                                                             <  1:1089635/126‑1 (MQ=255)
                             gTACAACGCCACCTTCCTGCAGTCGGTGATGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGt                                                                                                                                             >  2:1089635/1‑126 (MQ=255)
                              tACAACGCCACCTTCCTGCAGTCGGTGATGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAagag                                                                                                                             >  2:159211/1‑141 (MQ=255)
                              tACAACGCCACCTTCCTGCAGTCGGTGATGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAagag                                                                                                                             <  1:159211/141‑1 (MQ=255)
                              tACAACGCCACCTTCCTGCAGTCGGTGATGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAg                                                                                                                    >  1:166961/1‑150 (MQ=255)
                              tACAACGCCACCTTCCTGCAGTCGGTGATGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAg                                                                                                                    >  1:83958/1‑150 (MQ=255)
                                          ttCCTGCAGTCGGTGATGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAgg                                                                                                        <  2:166961/150‑1 (MQ=255)
                                                       tGATGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAgg                                                                                                                   <  1:669680/126‑1 (MQ=255)
                                                       tGATGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAgg                                                                                                                   >  2:669680/1‑126 (MQ=255)
                                                         aTGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCg                                                                                                                                  >  2:1285898/1‑109 (MQ=255)
                                                         aTGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCg                                                                                                                                  <  1:1285898/109‑1 (MQ=255)
                                                            ggCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCCCCGGATAGGACAGCGGCGTGGGTATCCCCGAAGCGCAGCGTGAGGCCATGTTCAGGGCGTTCGTGCGGGGCGAt                                                                                      >  1:672206/1‑150 (MQ=255)
                                                             gCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCggg                                                                                           <  1:604541/144‑1 (MQ=255)
                                                             gCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCggg                                                                                           >  2:604541/1‑144 (MQ=255)
                                                              cAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGa                                                                                    >  1:297229/1‑150 (MQ=255)
                                                               aaCCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGaa                                                                                   >  2:458449/1‑150 (MQ=255)
                                                                aCCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAAc                                                                                  <  1:501275/150‑1 (MQ=255)
                                                                 ccTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGaa                                                                                                                                <  2:847640/103‑1 (MQ=255)
                                                                 ccTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGaa                                                                                                                                >  1:847640/1‑103 (MQ=255)
                                                                 ccTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCAGCGTGGGTATCCGCCAAGAGCAGCGAGAGGCCATGTTCAGGCCCTTCGTGCGAGGCGACGAACg                                                                                 >  2:429270/1‑150 (MQ=255)
                                                                 ccTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGAGGCGTGGGTa                                                                                                                                        >  2:1328894/1‑95 (MQ=255)
                                                                 ccTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGAGGCGTGGGTa                                                                                                                                        <  1:1328894/95‑1 (MQ=255)
                                                                 ccTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGAGGCGTGGGTa                                                                                                                                        <  1:1329376/95‑1 (MQ=255)
                                                                 ccTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGAGGCGTGGGTa                                                                                                                                        >  2:1329376/1‑95 (MQ=255)
                                                                   tgctgcGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGc                                                                               >  2:99982/1‑150 (MQ=255)
                                                                     ctgcGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGt                                                                             <  1:99982/150‑1 (MQ=255)
                                                                           aaTGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAg                                                                       >  1:632814/1‑150 (MQ=255)
                                                                                        acacCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTgg                                                          <  2:783689/150‑1 (MQ=255)
                                                                                         cacCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTTGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCtt                                                                 <  1:319243/142‑1 (MQ=255)
                                                                                         cacCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCtt                                                                 >  2:319243/1‑142 (MQ=255)
                                                                                         cacCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTggg                                                         <  1:1088897/150‑1 (MQ=255)
                                                                                          acCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGAt                                                                                      >  1:1058588/1‑120 (MQ=255)
                                                                                          acCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGAt                                                                                      <  2:1058588/120‑1 (MQ=255)
                                                                                           ccGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGc                                                       <  1:321887/150‑1 (MQ=255)
                                                                                                aGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCg                                                  <  2:83958/150‑1 (MQ=255)
                                                                                                  tGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCt                                                >  2:320742/1‑150 (MQ=255)
                                                                                                  tGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCt                                                >  1:740865/1‑150 (MQ=255)
                                                                                                     gTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGTGTGGGTCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGTGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCTGGt                                             >  2:200755/1‑150 (MQ=255)
                                                                                                     gTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAgg                                                                                                        <  1:47120/91‑1 (MQ=255)
                                                                                                     gTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAgg                                                                                                        >  2:47120/1‑91 (MQ=255)
                                                                                                     gTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCTGGt                                             <  2:1329448/150‑1 (MQ=255)
                                                                                                          tCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGa                                                                        <  2:492735/118‑1 (MQ=255)
                                                                                                          tCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGa                                                                        >  1:492735/1‑118 (MQ=255)
                                                                                                          tCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCTGGTGCAGc                                        >  2:746547/1‑150 (MQ=255)
                                                                                                                tGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCTGGTGCAGCGGATCt                                  >  1:849690/1‑150 (MQ=255)
                                                                                                                    acTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCTGGTGCAGCGGATCTGcga                              >  1:276038/2‑150 (MQ=255)
                                                                                                                      tGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAggg                                                                     >  1:999993/1‑109 (MQ=255)
                                                                                                                      tGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAggg                                                                     <  2:999993/109‑1 (MQ=255)
                                                                                                                            ccAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCTGGTGCAGCGGATCTGCGACGACCAAg                      <  2:632814/150‑1 (MQ=255)
                                                                                                                              aaCGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCTGGTGCAGCGGATCTGCGACGACCAAGGt                    >  1:1399841/1‑150 (MQ=255)
                                                                                                                              aaCGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCTGGTGCAGCGGATCTGCGACGACCAAGGt                    >  1:1248989/1‑150 (MQ=255)
                                                                                                                              aaCGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTCGCGACGGGCATGGGCTGGGGCTTTCGCTGGTGCAGCGGATCTGCGACGACCAAGGt                    >  2:749695/1‑150 (MQ=255)
                                                                                                                                   tttCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCTGGTGCAGCGGATCTGCGGCGACCAAGGTTGGCg               >  2:992526/1‑150 (MQ=255)
                                                                                                                                     tCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCTGGTGCAGCGGATCTGCGACGACCAAGGTTGGCgtg             <  2:672206/150‑1 (MQ=255)
                                                                                                                                       aGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCTGGTGCAGCGGATCTGCGACGACCAAGGTTGGCGTGTc           <  1:458449/150‑1 (MQ=255)
                                                                                                                                                gACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCTGGTGCAGCGGATCTGCGACGACCAAGGTTGGCGTGTCACCCTGACa  >  1:787501/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GATGGCCGCGTGAGCGCTAGCCCTGTGCTGTACAACGCCACCTTCCTGCAGTCGGTGATGGGCAACCTGCTGCGCAATGCGGCGCACTACACCGACAGTGGGTATATCCGCCTGAGCCTGGAAGCCAACGGTTTCAGCGTAGAGGACAGCGGCGTGGGTATCCCCGAAGAGCAGCGTGAGGCCATGTTCAGGCCGTTCGTGCGGGGCGATGAACGGCGTGGCGAGGGGCTTGGGCTGGGGCTTTCGCTGGTGCAGCGGATCTGCGACGACCAAGGTTGGCGTGTCACCCTGACA  >  NC_002947/1043127‑1043420

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: