Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,816,594 T→A 8.5% I21N (ATT→AAT)  PP_2472 → MerR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,816,5940TA 8.5% 197.5 / 8.3 71I21N (ATT→AAT) PP_2472MerR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base A (4/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

CAAGACAGGGGAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCGT  >  NC_002947/2816447‑2816739
                                                                                                                                                   |                                                                                                                                                 
cAAGACAGGGGAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCAtt                                                                                                                                                  >  2:148896/1‑149 (MQ=255)
 aaGACAGGGGAAGAGATCCCGATCACCGCTCGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAACGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTg                                                                                                                                                 <  6:236793/149‑1 (MQ=255)
        gggAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGt                                                                                                                                          >  6:330692/1‑149 (MQ=255)
         ggAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCAATGGTGAAGta                                                                                                                                         <  3:24710/149‑2 (MQ=255)
                aTCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGt                                                                                                                                  >  2:364862/1‑149 (MQ=255)
                      aTCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTtgtgtg                                                                                                                            <  8:28731/149‑1 (MQ=255)
                         aCCGCACTGCGCGTCGTCACATCTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGTAACCAAGCCTTAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGcc                                                                                                                          <  5:16638/148‑1 (MQ=255)
                              aCGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAg                                                                                                                    <  7:89496/149‑1 (MQ=255)
                                ggCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGcc                                                                                                                  <  4:9174/149‑1 (MQ=255)
                                 gcgcGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCca                                                                                                                 >  4:167301/1‑149 (MQ=255)
                                   gcgTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCcaca                                                                                                               >  4:42938/1‑149 (MQ=255)
                                          cACCTTTCGGCCCCGGCAGAAGTTCAAGGCCCGGGTTGCGGCCTATGCGGGAACCAAGCCCTAACGCCGAACCGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGtgctg                                                                                                        <  2:300940/149‑1 (MQ=255)
                                                tCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTa                                                                                                   <  3:42938/148‑1 (MQ=255)
                                                tCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTAt                                                                                                  >  3:335113/1‑149 (MQ=255)
                                                  gTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTAtt                                                                                                 <  6:135993/148‑1 (MQ=255)
                                                  gTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTACCGTAAAGCCACACGTGCTGCGGTAtt                                                                                                 <  1:16305/148‑1 (MQ=255)
                                                   tCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTgg                                                                                               <  7:187660/149‑1 (MQ=255)
                                                   tCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTgg                                                                                               <  6:226467/149‑1 (MQ=255)
                                                   tCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTgg                                                                                               >  8:192243/1‑149 (MQ=255)
                                                     cAGGGCAGACGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACTAACTGCCCCCCAAACCGGGCAACCGGTATTTCACCAGTGGTGAAGTGATCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTAtttcga                                                                                             <  5:268706/149‑5 (MQ=255)
                                                     cAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTTGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTggg                                                                                              <  3:209604/148‑1 (MQ=255)
                                                      aGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAg                                                                                            >  2:192053/1‑149 (MQ=255)
                                                          cAGAAGTTGAAGGTCCGGTTTGAGGCGTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAgg                                                                                        <  8:268160/149‑1 (MQ=255)
                                                             aaGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTTAGCGAGTTGTGTGCCGTAAAGCAACACGTGCTGCGGTATTGGGAGCAGGAGt                                                                                     >  3:314047/1‑149 (MQ=255)
                                                                  gAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGATGCGGTATTGGGAGCAGGAGTTCGAc                                                                                >  7:121587/1‑149 (MQ=255)
                                                                   aaGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGTCGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGAcc                                                                               <  5:330692/149‑1 (MQ=255)
                                                                      gCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGc                                                                            >  4:233370/1‑149 (MQ=255)
                                                                        ccGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCt                                                                           >  1:182371/1‑148 (MQ=255)
                                                                               gAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCtg                                                                   >  8:238564/1‑149 (MQ=255)
                                                                               cgggCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCtg                                                                   <  1:279320/147‑1 (MQ=255)
                                                                                  gCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCtgtg                                                                 >  4:72206/1‑148 (MQ=255)
                                                                                   ccTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGaa                                                               <  3:101334/149‑1 (MQ=255)
                                                                                    cTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTTCTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGAcc                                                                               <  8:352900/132‑1 (MQ=255)
                                                                                    cTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGAcc                                                                               >  7:352900/1‑132 (MQ=255)
                                                                                    cTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAg                                                              <  1:148896/149‑1 (MQ=255)
                                                                                     tATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAgc                                                             >  3:2302/1‑149 (MQ=255)
                                                                                     tATGCTGGAACAAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAgc                                                             >  1:127790/1‑149 (MQ=255)
                                                                                          tGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCt                                                                    >  5:321898/1‑137 (MQ=255)
                                                                                          tGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCt                                                                    <  6:321898/137‑1 (MQ=255)
                                                                                           ggAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCggg                                                       >  8:8021/1‑149 (MQ=255)
                                                                                           ggAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCggg                                                       >  1:68563/1‑149 (MQ=255)
                                                                                            gAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCggg                                                       <  2:182371/148‑1 (MQ=255)
                                                                                               ccAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACAAGCTGAACCCTGTGAAGCGGCGGGGGaa                                                   >  3:227236/1‑149 (MQ=255)
                                                                                                cAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAg                                                                                      <  4:150051/113‑1 (MQ=255)
                                                                                                cAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAg                                                                                      >  3:150051/1‑113 (MQ=255)
                                                                                                cAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAAc                                                  <  3:321950/149‑1 (MQ=255)
                                                                                                     cATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACcggcg                                             <  1:364862/149‑1 (MQ=255)
                                                                                                         acgacgAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTAt                                         >  2:157529/1‑149 (MQ=255)
                                                                                                            acgaACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTAc                                      <  8:223277/149‑1 (MQ=255)
                                                                                                                      ccccATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTAc                                      >  5:47343/1‑139 (MQ=255)
                                                                                                                      ccccATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTAc                                      <  6:47343/139‑1 (MQ=255)
                                                                                                                      ccccATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAg                             >  8:17947/1‑148 (MQ=255)
                                                                                                                       cccATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGa                           >  1:361830/1‑149 (MQ=255)
                                                                                                                        ccATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGAt                          >  1:155048/1‑149 (MQ=255)
                                                                                                                        ccATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGAt                          >  7:136264/1‑149 (MQ=255)
                                                                                                                           tACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGAtgtg                       <  3:233370/149‑1 (MQ=255)
                                                                                                                           tACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGAtgtg                       <  4:335113/149‑1 (MQ=255)
                                                                                                                              cGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCtg                    <  3:167301/149‑1 (MQ=255)
                                                                                                                              cGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCt                     >  3:141190/1‑148 (MQ=255)
                                                                                                                               gggCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCtga                   <  1:192053/149‑1 (MQ=255)
                                                                                                                                ggCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCtgat                  >  7:46462/1‑149 (MQ=255)
                                                                                                                                 gCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCtgat                  >  1:265667/1‑148 (MQ=255)
                                                                                                                                   aaaCGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTgatgat               >  5:256152/1‑149 (MQ=255)
                                                                                                                                    aaCGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATc              <  3:72206/149‑1 (MQ=255)
                                                                                                                                      gctcaCTTCACCAATGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAgg                            >  3:153477/5‑133 (MQ=255)
                                                                                                                                      gctcaCTTCACCAATGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAgg                            <  4:153477/129‑1 (MQ=255)
                                                                                                                                        tcaCTTCACCAATGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGc           >  7:196963/3‑148 (MQ=255)
                                                                                                                                        tcaCTTCACCAATGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGc           >  2:66655/3‑148 (MQ=255)
                                                                                                                                         caCTTCACCAATGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGc           >  2:282618/2‑147 (MQ=255)
                                                                                                                                                ccATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCGt  <  2:127790/149‑1 (MQ=255)
                                                                                                                                                ccATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCGt  <  1:157529/149‑1 (MQ=255)
                                                                                                                                                ccATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCGt  <  8:46462/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CAAGACAGGGGAAGAGATCCCGATCACCGCACGGCGCGTCGTCACCTTTCGTCCAGGGCAGAAGTTGAAGGCCCGGGTTGAGGCCTATGCTGGAACCAAGCCATAACGACGAACTGCCCCCCATACCGGGCAAACGCTACTTCACCATTGGTGAAGTGAGCGAGTTGTGTGCCGTAAAGCCACACGTGCTGCGGTATTGGGAGCAGGAGTTCGACCAGCTGAACCCTGTGAAGCGGCGGGGGAACCGGCGGTATTACCAGCGCCAGGATGTGCTGATGATCCGCCAGATTCGT  >  NC_002947/2816447‑2816739

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