Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 697,214 T→C 12.2% intergenic (+141/‑610) mmsA‑I → / → PP_16SD methylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947697,2140TC12.2% 225.6 / 20.3 90intergenic (+141/‑610)mmsA‑I/PP_16SDmethylmalonate‑semialdehyde dehydrogenase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (46/33);  new base C (4/7);  total (50/40)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.17e-01

CTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGC  >  NC_002947/697068‑697360
                                                                                                                                                  |                                                                                                                                                   
cTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCg                                                                                                                                                   >  3:107536/1‑149 (MQ=255)
 tGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATc                                                                                                                                                    >  3:144122/1‑147 (MQ=255)
 tGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATc                                                                                                                                                    >  2:262173/1‑147 (MQ=255)
 tGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTGAGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATc                                                                                                                                                    >  1:288299/1‑147 (MQ=255)
   tAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAgc                                                                                                                                                <  5:370566/149‑1 (MQ=255)
         aaGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCgccgcc                                                                                                                                          <  1:151746/149‑1 (MQ=255)
              gATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCcgcgc                                                                                                                                     >  6:376486/1‑149 (MQ=255)
              gATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCcgcgc                                                                                                                                     >  6:265187/1‑149 (MQ=255)
              gATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACGGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCcgcgc                                                                                                                                     >  8:70863/1‑149 (MQ=255)
                    gCGGGGAGGCCTTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCt                                                                                                                               <  3:349596/149‑1 (MQ=255)
                     cGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                              <  5:265187/149‑1 (MQ=255)
                       gggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCga                                                                                                                            <  4:107536/149‑2 (MQ=255)
                        ggAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGTGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATTGAGCGCCGCGCGGGCGGCGCTCgat                                                                                                                           <  4:259577/149‑3 (MQ=255)
                         gAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGtc                                                                                                                          <  5:124402/149‑1 (MQ=255)
                          aGGCCGTCAGGCCCCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCGAtct                                                                                                                         <  3:213441/149‑1 (MQ=21)
                                   ggCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCAGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGgcgc                                                                                                                >  5:379354/1‑149 (MQ=255)
                                    gCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGTAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAgg                                                                                                                   <  2:64156/145‑1 (MQ=255)
                                    gCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGAGCCGCCCGCGCGGCGCTCGGTCTCACAgg                                                                                                                   >  1:64156/1‑145 (MQ=255)
                                       tCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGa                                                                                                             >  3:157558/1‑148 (MQ=255)
                                       tCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGa                                                                                                             >  3:258054/1‑148 (MQ=255)
                                        ccccGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGaaa                                                                                                           <  2:288299/149‑1 (MQ=255)
                                              tttGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCg                                                                                                     >  5:273457/1‑149 (MQ=255)
                                               ttGTTTTTGGGGGTTGGTTGTTGGTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGGCACGACCATTGCAGTGCCTGTGTGATGGAGGGCGGGCCGCGCGGCGCTCGGTGTCACATGCGCTGAATCTGTCGt                                                                                                    >  6:10701/1‑149 (MQ=255)
                                                 gTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTgg                                                                                                  >  4:269106/1‑149 (MQ=255)
                                                 gTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTg                                                                                                   <  5:169293/148‑1 (MQ=255)
                                                 gTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTg                                                                                                   <  4:144122/148‑1 (MQ=255)
                                                  tttttGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCgcgg                                                                                                                                   >  2:373698/1‑115 (MQ=255)
                                                    tttGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCg                                                                                                >  6:216905/1‑148 (MQ=255)
                                                      tGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACGACAGTTTCAGCGCCTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCACTGCAAGGGTTGTGGCgga                                                                                              <  5:3527/148‑3 (MQ=14)
                                                       gggggTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGACTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAac                                                                                             >  8:188526/1‑148 (MQ=255)
                                                          ggttggttGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACtt                                                                                        <  7:188526/149‑1 (MQ=255)
                                                            ttggttgTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTgg                                                                                      <  5:376486/149‑1 (MQ=255)
                                                              ggttgTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCa                                                                                    >  8:198161/1‑149 (MQ=255)
                                                                  gttgCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGccc                                                                                >  6:225881/1‑149 (MQ=255)
                                                                   ttgCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGcccc                                                                               <  6:379354/149‑1 (MQ=255)
                                                                   ttgCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGccc                                                                                <  5:216905/148‑1 (MQ=255)
                                                                     gCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGCCGAACACTTGGCAGCCCCGa                                                                             >  8:109820/1‑149 (MQ=255)
                                                                            gggAGGAAGTGTTTTGAGAGGGTGTTGGCTGGGAGATGTGTCCGCCACAACCCTTTCAGTGCCTGTGAGATCGAGCGCCGCCGGCGCTGCGCTCGGTCTCACAGGCGCTGAATCTGTGGTGTCGAACACTTGGCAGCCCCGATACGGtg                                                                      <  5:66056/149‑2 (MQ=25)
                                                                            gggAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGt                                                                       <  5:225881/148‑1 (MQ=255)
                                                                            gggAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTc                                                                      >  2:23462/1‑149 (MQ=255)
                                                                             ggaggaACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTgg                                                                                      >  4:189171/1‑132 (MQ=255)
                                                                              gaggaACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCtt                                                                    >  3:253885/1‑149 (MQ=255)
                                                                                     tGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCt                                                             >  8:93736/1‑149 (MQ=255)
                                                                                       ttttGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCtt                                                            <  6:273457/148‑1 (MQ=255)
                                                                                       ttttGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTta                                                           >  3:302355/1‑149 (MQ=255)
                                                                                       ttttGAGAGGGGGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACATTGGCAGCCCCGATACGGTCTTTTcccc                                                              <  1:23462/147‑1 (MQ=255)
                                                                                         ttGACAGGCTATTGGCTCTGAGATGTGTCCGCCACAACCCTTCCAGTGCCTGTGAGATCGAGCGACGCCCCCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTtata                                                         >  4:162436/1‑149 (MQ=21)
                                                                                            agagGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATAtt                                                       >  5:97109/1‑148 (MQ=255)
                                                                                               gggTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCAc                                                   >  1:24440/1‑149 (MQ=255)
                                                                                                ggTGTTGGCTGTGAGATGTGCCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACt                                                  <  4:267764/149‑1 (MQ=255)
                                                                                                    ttGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGACCGCCGCCCGCGCGGCGCTCGGTCTCACATGCGCTGAAACTTTCGTGGCGAACACTTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTc                                               <  1:262173/149‑1 (MQ=25)
                                                                                                    tgggCTGTGAGATGGGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTc                                               <  3:200674/146‑1 (MQ=255)
                                                                                                      ggCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCGATTCCCCTTATATTTCACTTTCtt                                             >  3:42168/1‑148 (MQ=255)
                                                                                                       gCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTg                                            >  8:193022/1‑148 (MQ=255)
                                                                                                       gCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTg                                            >  6:33040/1‑148 (MQ=255)
                                                                                                         tgtgAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGGTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGGCTTTTCCCCTTATATTTCACTTTCTTGaaa                                         >  2:244811/1‑149 (MQ=255)
                                                                                                           tcacaTGTGACCGCCACAACCCCTGCACACCCAGACAGATCCAGCGCCCCCCGCTCGGCGCTCGGTCTCACAGGCGCAGAAACTCTCCTGGCGAACACTTGGCAGCCCCGATACGGTCTATTCCCCTTATATATCACATTCATCaaa                                         <  5:10701/143‑1 (MQ=12)
                                                                                                              gaTGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTaaa                                    <  7:93736/149‑1 (MQ=255)
                                                                                                              gaTGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTaaa                                    <  5:33040/149‑1 (MQ=255)
                                                                                                              gaTGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGAGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTaaa                                    <  7:70863/149‑1 (MQ=255)
                                                                                                               aTGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAg                                   <  4:302355/149‑1 (MQ=255)
                                                                                                                 gtgtCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGt                                 >  2:214534/1‑149 (MQ=255)
                                                                                                                 gtgtCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGt                                 >  8:133786/1‑149 (MQ=255)
                                                                                                                  tgtCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGtt                                >  6:136224/1‑149 (MQ=255)
                                                                                                                      cGCCACAACCCTTGCAGTGCCTGTGTGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGaca                            <  7:27418/149‑1 (MQ=255)
                                                                                                                       gCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCtt                                             <  4:353220/131‑1 (MQ=255)
                                                                                                                       gCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCtt                                             >  3:353220/1‑131 (MQ=255)
                                                                                                                        ccACAACCCTTGCAGTGCCTGGGAGATGGTGCGCCGCCCGCGCGGCGGTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATTCGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACAcc                          <  4:42168/149‑1 (MQ=18)
                                                                                                                          acaACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACAccc                         >  7:23468/1‑148 (MQ=255)
                                                                                                                          acaACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCct                        >  1:154384/1‑149 (MQ=255)
                                                                                                                            aaCCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGCACACATGGCAGCCCCGATACGGTCATTTCCCCTTATATTTGAATTTCTTGaaa                                         >  6:50854/1‑130 (MQ=18)
                                                                                                                            aaCCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGaaa                                         >  4:152918/1‑130 (MQ=255)
                                                                                                                                cTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTAGGCAGCCCCGCTACGGTCTTTTCCCCTTTTATTTCACTTTCTTGAAATTAAAGGATGACACCGTCTGGaa                  >  3:237020/1‑149 (MQ=17)
                                                                                                                                 ttGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTa                                      <  7:281092/128‑1 (MQ=255)
                                                                                                                                 ttGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTa                                      >  8:281092/1‑128 (MQ=255)
                                                                                                                                   gCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATc                <  7:193022/148‑1 (MQ=255)
                                                                                                                                    cAGCGCCTATGAGACCGAGCGCCGCCCGCGCGGCGCTCGGTTTCATAGGCGGTGAAtgtctcaagccggacgcgggcggcgctcgatctcacaggcgctgcaagggtgttggcggtgagatgtgtccgccacaacccttgc                       <  4:68113‑M2/141‑86 (MQ=255)
                                                                                                                                    cAGCGCCTATGAGACCGAGCGCCGCCCGCGCGGCGCTCGGTTTCATAGGCGGTGAAtgtctcaagccggacgcgggcggcgctcgatctcacaggcgctgcaagggtgttggcggtgagatgtgtccgccacaacccttgc                       >  3:68113‑M2/1‑56 (MQ=255)
                                                                                                                                         ccTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTta          <  8:23468/148‑1 (MQ=255)
                                                                                                                                          cTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCACTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATAttt                                                      >  7:47812/1‑103 (MQ=11)
                                                                                                                                          cTGTGAGACCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCACTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATAttt                                                      <  8:47812/103‑1 (MQ=11)
                                                                                                                                            atgaGACCGAGCGCCGCCCGCGCGGCGCTCGGTTTCATAGGCGGTGAAtgtctcaagccggacgcgggcggcgctcgatctcacaggcgctgcaagggtgttggcggtgagatgtgtccgccacaacccttg                        >  1:142827‑M2/2‑48 (MQ=255)
                                                                                                                                            atgaGACCGAGCGCCGCCCGCGCGGCGCTCGGTTTCATAGGCGGTGAAtgtctcaagccggacgcgggcggcgctcgatctcacaggcgctgcaagggtgttggcggtgagatgtgtccgccacaacccttg                        <  2:142827‑M2/131‑85 (MQ=255)
                                                                                                                                             tgAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCACCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATg     >  2:206194/1‑149 (MQ=255)
                                                                                                                                               agaTCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACtt                                                                                        <  1:380306/64‑1 (MQ=255)
                                                                                                                                               agaTCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACtt                                                                                        >  2:380306/1‑64 (MQ=255)
                                                                                                                                               agaTCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCGCCTTATAATgcg   >  7:316355/1‑149 (MQ=255)
                                                                                                                                               agaTCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATaa       >  6:338064/1‑145 (MQ=255)
                                                                                                                                               agaTCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATaa       <  5:338064/145‑1 (MQ=255)
                                                                                                                                                gaTCGAGCGCCGCCCGCGCGGCGCGCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTGGAATCCCCTTAGACTgcgc  >  7:63322/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTGTAAGGGAAGTAGATGAAGCGGGGAGGCCGTCAGGCCTCCCCGCTTTGTTTTTGGGGGTTGGTTGTTGCTGTTGGGGAGGAACTGTTTTGAGAGGGTGTTGGCTGTGAGATGTGTCCGCCACAACCCTTGCAGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGGTCTCACAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCCCCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACACCCTCTTGAATCCCCTTATAATGCGC  >  NC_002947/697068‑697360

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