Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,137,332 C→A 7.2% L90L (CTG→CTT)  PP_4522 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,137,3320CA 7.2% 216.1 / 6.5 69L90L (CTG→CTT) PP_4522LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (34/30);  new base A (2/3);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-01

GGGCTGGACACCGATAGATGCAGGTGCACCGCCGGGTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAG‑‑CCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGG  >  NC_002947/5137192‑5137476
                                                                                                                                            |                                                                                                                                                   
gggCTGGACACCGATAGATGCAGGTGCACCGCCGGGTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTGAAGggcgaa                                                                                                                                             <  7:27586/149‑3 (MQ=255)
 ggCTGGACACCGATAGATGCAGGTGCACCGCCGGGTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAgg                                                                                                                                                  >  4:223875/1‑143 (MQ=255)
 ggCTGGACACCGATAGATGCAGGTGCACCGCCGGGTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAgg                                                                                                                                                  <  3:223875/143‑1 (MQ=255)
 ggCTGGACACCGATAGATGCAGGTGCACCGCCGGGTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTGAAGggcgaac                                                                                                                                            <  4:295240/149‑4 (MQ=255)
  gCTGGACACCGATAGATGCAGGTGCACCGCCGGGTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTGGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTGAAGggcgaact                                                                                                                                           <  5:266261/149‑5 (MQ=255)
  gCTGGACACCGATAGATGCAGGTGCACCGCCGGGTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGggcggcgg                                                                                                                                           <  8:306355/149‑1 (MQ=255)
      gACACCGATAGATGCAGGTGCACCGCCGGGTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGCCGGCGGCATa                                                                                                                                       >  2:203786/1‑149 (MQ=255)
      gACACCGATAGATGCAGGTGCACCGCCGGGTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGGTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATa                                                                                                                                       >  3:48713/1‑149 (MQ=255)
       acacCGATAGATGCAGGAGCACCGCCGGGTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATAc                                                                                                                                      >  3:184686/1‑149 (MQ=255)
             aTAGATGCAGGTGCACCGCCGGGTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGggcggcg                                                                                                                                            >  1:166007/1‑137 (MQ=255)
             aTAGATGCAGGTGCACCGCCGGGTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGggcggcg                                                                                                                                            <  2:166007/137‑1 (MQ=255)
                gATGCAGGTGCACCGCCGGGTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGaaa                                                                                                                             <  3:202209/149‑1 (MQ=255)
                      ggTGCACCGCCGGGTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTc                                                                                                                        >  3:260901/1‑148 (MQ=255)
                       gTGCACCGCCGGGTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTc                                                                                                                        >  7:48319/1‑147 (MQ=255)
                        tGCACCGCCGGGTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAg                                                                                                                     >  4:269760/1‑149 (MQ=255)
                           aCCGCAGGGTGATCCTGGCTGTAGGCACCGATGGATTCTGCCAACGGCAAGGCGCGTTCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCATGGCGGCGGCATACTGCTCGAAACCGTCCAGCTc                                                                                                                  <  2:193926/149‑1 (MQ=255)
                              gccgGGTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGc                                                                                                               <  4:260901/149‑1 (MQ=255)
                                   gTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCg                                                                                                                 <  2:168886/142‑1 (MQ=255)
                                   gTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCg                                                                                                                 >  1:168886/1‑142 (MQ=255)
                                   gTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCcagca                                                                                                          <  4:111074/149‑1 (MQ=255)
                                    tGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCcagcag                                                                                                         >  6:243227/1‑149 (MQ=255)
                                       tCCTGGCTGTAGGCACCGATGGCTTCGGACAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCCTACTGCTCGAAACCGTCCAGCTCGGCCAGCAGgcg                                                                                                      >  3:50690/1‑149 (MQ=255)
                                                tAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTc                                                                                             >  8:184740/1‑149 (MQ=255)
                                                   gCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCt                                                                                           >  2:14137/1‑148 (MQ=255)
                                                   gCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCt                                                                                           >  3:300998/1‑148 (MQ=255)
                                                   gCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTg                                                                                          <  1:221376/149‑1 (MQ=255)
                                                     aCCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGAt                                                                                        <  2:63853/149‑1 (MQ=255)
                                                     aCCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGAt                                                                                        <  4:184686/149‑1 (MQ=255)
                                                      ccGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATg                                                                                       <  5:243227/149‑1 (MQ=255)
                                                                 ggCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGc                                                                           <  1:203786/149‑1 (MQ=255)
                                                                         gCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTg                                                                   <  4:300998/149‑1 (MQ=255)
                                                                           aaGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTACCTGATGG‑AACAGCTCGCCCTTGCTGGt                                                                 <  1:187722/149‑1 (MQ=255)
                                                                            aGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTc                                                                <  1:14137/149‑1 (MQ=255)
                                                                             ggCGCTATCACCGACGGTGGAGGCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCa                                                               <  7:269617/149‑1 (MQ=255)
                                                                             ggCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGt                                                                 >  7:36701/1‑147 (MQ=255)
                                                                              gcgcGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAg                                                              >  1:293024/1‑149 (MQ=255)
                                                                                gcgATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGc                                                            <  3:269760/149‑1 (MQ=255)
                                                                                 cgATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCt                                                           >  8:18286/1‑149 (MQ=255)
                                                                                    tCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑aa                                                                                   <  8:2159/122‑1 (MQ=255)
                                                                                    tCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑aa                                                                                   >  7:2159/1‑122 (MQ=255)
                                                                                    tCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGaa                                                        <  3:216368/149‑1 (MQ=255)
                                                                                       ccGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAG‑‑cc                                                   >  2:90993/1‑149 (MQ=255)
                                                                                        cGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAG‑‑ccc                                                  >  3:19796/1‑149 (MQ=255)
                                                                                            ggTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAG‑‑ccccccc                                              <  2:293024/149‑1 (MQ=255)
                                                                                                 aGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑GACAGCTCGCCCTTGCTGGTCAGGCTGAAG‑‑ccccccccgacc                                         >  1:92684/1‑144 (MQ=255)
                                                                                                     gATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGAt                                                                                        <  4:266935/101‑1 (MQ=255)
                                                                                                     gATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGAt                                                                                        >  3:266935/1‑101 (MQ=255)
                                                                                                       tCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAGC‑CCCCCCCGACCCCGGTg                                    <  7:184740/149‑1 (MQ=255)
                                                                                                         aCCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAGCCCCCCCCCGACCCCGGTgg                                   <  4:48713/149‑1 (MQ=255)
                                                                                                           cccAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAGC‑CCCCCCCGACCCCGGTGGca                                 >  6:120323/1‑148 (MQ=255)
                                                                                                           cccAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAG‑‑CCCCCCCGACCCCGGTGGcac                                >  3:54258/1‑148 (MQ=255)
                                                                                                             cAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAGC‑CCCCCCCGACCCCGGTGGCacaa                              >  8:72149/1‑149 (MQ=255)
                                                                                                               aGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAGC‑CCCCCCCGACCCCGGTGGCacaaca                            <  7:72149/149‑1 (MQ=255)
                                                                                                                ggTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAGC‑CCCCCCCGACCCCGGTGGCACAacac                           >  3:188281/1‑149 (MQ=255)
                                                                                                                 gTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGGCAACAGCTCGCCCTTGCTGGTCAGGCTGAAGC‑CCCCCCCGACCCCGGTGGCACAacac                           <  1:90993/149‑1 (MQ=255)
                                                                                                                 gTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAGC‑CCCCCCCGCCCCCGGTGGCACACCACg                          >  4:306428/1‑149 (MQ=255)
                                                                                                                  ttGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAGC‑CCCCCCCGACCCCGGTGGCACAACACGa                         >  4:235325/1‑149 (MQ=255)
                                                                                                                      gggTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAGC‑CCCCCCCGACCCCGGTGGCACAACACGATAcc                     >  8:301885/1‑149 (MQ=255)
                                                                                                                      gggTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAGC‑CCCCCCCGACCCCGGTGGCACAACACGATAcc                     <  8:48319/149‑1 (MQ=255)
                                                                                                                           ccACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAGC‑CCCCCCCGCCCCCGGTGGCACAACCCGATACCCAAgg                >  5:12747/1‑149 (MQ=255)
                                                                                                                            cACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAGC‑CCCCCCCGACCCCGGTGGCACAACACGATACCCAAGgc               <  3:235325/149‑1 (MQ=255)
                                                                                                                             aCGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAGC‑CCCCCCCGACCCCGGTGGCACAACACGATACCCAAGgcg              >  5:157995/1‑149 (MQ=255)
                                                                                                                               gCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAGC‑CCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGcc            >  8:173657/1‑149 (MQ=255)
                                                                                                                                 ccgccGCCCTGAAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAGC‑CCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCt          >  5:129635/5‑149 (MQ=255)
                                                                                                                                  cgccGCCCTGAAGGGCGGCGGCAAACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAGC‑CCCCCCCGCCCCCGGTGGCACAACCCGAAACCCAAGGCGCCCTc         >  8:186093/4‑149 (MQ=255)
                                                                                                                                     cGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAGC‑CCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAg      >  5:219112/1‑149 (MQ=255)
                                                                                                                                      gCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAGC‑CCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGc     <  7:301885/149‑1 (MQ=255)
                                                                                                                                       cccTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAG‑‑CCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTg   >  6:156706/1‑149 (MQ=255)
                                                                                                                                         cTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑GACAGCTCGCCCTTGCTGGTCAGGCTGAAGC‑CCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTgg  <  2:92684/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
GGGCTGGACACCGATAGATGCAGGTGCACCGCCGGGTGCTCCTGGCTGTAGGCACCGATGGCTTCGGCCAACGGCAAGGCGCGATCACCGACGGTGGAGTCGATCACCCCAAGGTTGAGGGTGCCACGCAGTTCGCCCTTCAGGGCGGCGGCATACTGCTCGAAACCGTCCAGCTCGGCCAGCAGGCGCAGGGTTTCCTGATGG‑AACAGCTCGCCCTTGCTGGTCAGGCTGAAG‑‑CCCCCCCGACCCCGGTGGCACAACACGATACCCAAGGCGCCCTCCAGCTGG  >  NC_002947/5137192‑5137476

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