Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 564,453 T→C 6.2% intergenic (+105/‑68) rplQ → / → katA 50S ribosomal protein L17/catalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947564,4530TC 6.2% 189.8 / 2.1 64intergenic (+105/‑68)rplQ/katA50S ribosomal protein L17/catalase
Reads supporting (aligned to +/- strand):  ref base T (35/25);  new base C (2/2);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCACCTGTAGCCGACAACCAGAACTCCCGCTCCGCCGGCCCGCGCGGTCC  >  NC_002947/564307‑564600
                                                                                                                                                    |                                                                                                                                                   
gTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTaa                                                                                                                                                   >  1:263256/1‑149 (MQ=255)
gTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTaa                                                                                                                                                   <  3:131398/149‑1 (MQ=255)
           ggTCGGCGGTGCTGTATTAGCTGGTTAGTAAGAAGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTTTTTATACATCTATTGTTTACTTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGt                                                                                                                                         <  1:291039/149‑1 (MQ=255)
           ggTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTaa                                                                                                                                                   >  3:377242/1‑138 (MQ=255)
           ggTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTaa                                                                                                                                                   <  4:377242/138‑1 (MQ=255)
            gTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCa                                                                                                                                       <  6:75707/149‑1 (MQ=255)
                gcggTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATAc                                                                                                                                   >  8:147918/1‑149 (MQ=255)
                gcggTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATAc                                                                                                                                   >  8:139301/1‑149 (MQ=255)
                        gTAGAAGCTGCTGAGTATGAAGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctctct                                                                                                                            <  7:324917/148‑1 (MQ=255)
                           gAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGCAATCTTGCTCGTCAATACTCTCTCTa                                                                                                                           <  7:341479/146‑1 (MQ=255)
                           gAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGCAATCTTGCTCGTCAATACTCTCTCTa                                                                                                                           >  8:341479/1‑146 (MQ=255)
                            aaGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCAAGTGAACGGGTTTTTGTTTATCCATCTATTGTTGATTTCGAGTGATGTAAGTCATTTAATTAATTTGTCCCTGTAATCTTGCTCGTCACTACTCTCGCTACacc                                                                                                                       >  4:243156/1‑147 (MQ=255)
                                 gctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGCAATCTTg                                                                                                                                              <  8:187376/121‑1 (MQ=255)
                                 gctgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGCAATCTTg                                                                                                                                              >  7:187376/1‑121 (MQ=255)
                                   tgAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAg                                                                                                                >  3:80896/1‑149 (MQ=255)
                                      gTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATATATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGc                                                                                                               <  8:115835/147‑1 (MQ=255)
                                      gTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGc                                                                                                               >  7:115835/1‑147 (MQ=255)
                                       tAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCgg                                                                                                            >  7:142781/1‑149 (MQ=255)
                                       tAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCgg                                                                                                            >  6:40644/1‑149 (MQ=255)
                                         aGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTATCTACATCCGATTAGCCGGcc                                                                                                          >  6:218578/1‑149 (MQ=255)
                                          gACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATAc                                                                                                                                   <  1:102333/123‑1 (MQ=255)
                                          gACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATAc                                                                                                                                   >  2:102333/1‑123 (MQ=255)
                                          gACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGccc                                                                                                         >  8:19378/1‑149 (MQ=255)
                                                       cTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGa                                                                                              <  7:147918/147‑1 (MQ=255)
                                                        aaCAAAAAACCGGGCACTAGGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGAcc                                                                                                >  7:333428/2‑144 (MQ=255)
                                                        aaCAAAAAACCGGGCACTAGGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGAcc                                                                                                <  8:333428/143‑1 (MQ=255)
                                                           aaaaaaCCGGGCACTAGGCCCGGTTTTCTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATg                                                                                          >  5:77427/3‑147 (MQ=255)
                                                           aaaaaaCCGGGCACTAGGCCCGGGTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGATCGGCAATACTATCTCTACAGCAGATTAGCCGGCCCTGGAAGACCGAAATg                                                                                          >  2:93147/3‑147 (MQ=255)
                                                           aGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGga                                                                                        <  7:19378/149‑1 (MQ=255)
                                                           aGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTAAATACTCTCTCTACAGCCGAGTAGCCGGCCGTTGAAGACCGAAATg                                                                                          >  8:97442/1‑147 (MQ=255)
                                                            gAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTc                                                                                                                                           <  7:254854/97‑1 (MQ=255)
                                                            gAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTc                                                                                                                                           >  8:254854/1‑97 (MQ=255)
                                                            gAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCATGTAATATTGATCGTCTATACTGTCTCTACAGCAGATTATAAGGAACTTGTAGTCCTAAATGgag                                                                                       >  5:292104/1‑149 (MQ=255)
                                                             aaaaCCGGGCCTAGTGGCCGGTTTTTTGTTTATAGATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGgaga                                                                                      <  8:84613/149‑1 (MQ=255)
                                                               aaCCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGc                                                                                    >  3:351634/1‑149 (MQ=255)
                                                                 ccGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGccc                                                                                  <  5:40644/149‑1 (MQ=255)
                                                                 ccGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGcc                                                                                   >  5:168206/1‑148 (MQ=255)
                                                                   gggCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCAt                                                                                <  1:154578/149‑1 (MQ=255)
                                                                            tGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGa                                                                        <  5:310575/149‑1 (MQ=255)
                                                                              cccGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGAt                                                                       <  8:142781/148‑1 (MQ=255)
                                                                               ccGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTACTGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGACCAAGAtt                                                                      <  7:97442/148‑1 (MQ=255)
                                                                                    tttttGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGCCCGAAATGGAGCGCCCATGAGCAAGATTCTcacc                                                                >  8:364937/1‑149 (MQ=255)
                                                                                                cATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCgg                                                     <  5:218578/148‑1 (MQ=255)
                                                                                                    tATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCAc                                                <  7:139301/149‑1 (MQ=255)
                                                                                                      ttGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCACCt                                              >  3:111401/1‑149 (MQ=255)
                                                                                                      ttGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCACCt                                              >  8:169379/1‑149 (MQ=255)
                                                                                                                 cTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTggt                                                 >  7:10992/1‑134 (MQ=255)
                                                                                                                 cTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTggt                                                 <  8:10992/136‑3 (MQ=255)
                                                                                                                  tATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCACCTGTAGCCGACaa                                   >  4:264990/1‑149 (MQ=255)
                                                                                                                    ttGATTTAAGTCATGTAATGAATTTGCCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCACCTGTAGCCGACAAc                                  <  6:292104/148‑1 (MQ=255)
                                                                                                                          tAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCACCTGTAGCCGACAACCAGAAc                            <  2:263256/148‑1 (MQ=255)
                                                                                                                             gTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCACCTGTAGCCGACAACCAGAACTccc                        >  2:228722/1‑149 (MQ=255)
                                                                                                                              tCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCACCTGTAGCCGACAACCAGAACTCCCg                       >  6:20944/1‑149 (MQ=255)
                                                                                                                                   tAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCACCTGTAGCCGACAACCAGAACTCCCGCTccg                  >  1:206778/1‑149 (MQ=255)
                                                                                                                                     aTGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCACCTGTAGCCGACAACCAGAACTCCCGCTccgcc                >  7:297548/1‑149 (MQ=255)
                                                                                                                                     aTGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCACCTGTAGCCGACAACCAGAACTCCCGCTacgcc                >  7:292243/1‑149 (MQ=255)
                                                                                                                                       gAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACa                                                                                                                         <  4:287159/42‑1 (MQ=255)
                                                                                                                                       gAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCACCTGTAGCCGACAACCAGAACTCCCGCTCCGCCgg              >  2:383040/1‑149 (MQ=255)
                                                                                                                                             gTCCCTGTAATTTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCACCTGTAGCCGACAACCAGAACTCCCGCTCCGCCGGCCcgcg        >  8:259803/1‑149 (MQ=255)
                                                                                                                                               cccTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCACCTGTAGCCGACAACCAGAACTCCCGCTCCGCCGGCCcgcgc       <  4:80896/148‑1 (MQ=255)
                                                                                                                                                 cTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCACCTGTAGCCGACAACCAGAACTccc                        <  2:383239/129‑1 (MQ=255)
                                                                                                                                                 cTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCACCTGTAGCCGACAACCAGAACTccc                        >  1:383239/1‑129 (MQ=255)
                                                                                                                                                 cTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCACCTGTAGCCGACAACCAGAACTCCCGCTCCGCCGGCCCGCGCGGt    >  5:378010/1‑149 (MQ=255)
                                                                                                                                                   gTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCACCTGTAGCCGACAACCAGAACTCCCGCTCCGCCGGCCCGCGCGGTcc  >  5:260938/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGAGCAAGATTCTCACCACCGCTAGCGGTGCACCTGTAGCCGACAACCAGAACTCCCGCTCCGCCGGCCCGCGCGGTCC  >  NC_002947/564307‑564600

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: