Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,351,519 C→T T48T (ACG→ACA)  dcuS ← sensor histidine kinase DcuS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,351,5190CT97.7% 139.2 / NA 43T48T (ACG→ACA) dcuSsensor histidine kinase DcuS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (22/20);  minor base G (0/1);  total (22/21)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGA  >  NC_000913/4351275‑4351751
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                         
tCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCat                                                                                                                                                                                                                                      <  1:399000/250‑1 (MQ=255)
tCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCat                                                                                                                                                                                                                                      <  1:219784/250‑1 (MQ=255)
              aTGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCg                                                                                                                                                                                                                        >  1:411944/1‑250 (MQ=255)
                             aaTGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGATAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATg                                                                                                                                                                                                        >  1:377762/1‑251 (MQ=255)
                               tGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGc                                                                                                                                                                                                       >  1:236166/1‑250 (MQ=255)
                                            tgCTGGACTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGACGTAAATCAGATGCACCAC‑CAATAGCAc                                                                                                                                                                                        >  2:104170/2‑251 (MQ=255)
                                              cTGGGCTTCAGGATCCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAAGCAGATGCACCAC‑CAATAGCACCg                                                                                                                                                                                      >  1:325892/1‑251 (MQ=255)
                                                 ggtttCAGGATTCGAGTAGCGAAGTCTTTGTATATCGGTAACTGCAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGGGAGGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAATCGCGGAGGGATTTCCGGCGAGTCGGCAAGGGTACGCGCCACTGACGGTGCCTTGTTGGCTAGCCCAGGTCGTGTCATATCAGTGATTTGCGAGAGGTAAATCAGAGGCACTCC‑CAATAGAACCgaga                                                                                                                                                                                   <  2:460963/248‑1 (MQ=255)
                                                       aGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGGGAAGTAAATCAGATGAACCAC‑CAATAGCACCGAGAACAGt                                                                                                                                                                              >  2:466371/1‑250 (MQ=255)
                                                        ggATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATA‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGCTTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTAc                                                                                                                                                                            >  1:237194/1‑251 (MQ=255)
                                                          aTGCGAGTAGCGAAGACATTGCATATCGGTAACGACAATAAACAGCAGATC‑GTTTCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGAAGGATTTCCGGCGAGTCGGAGAGGGTACGCGACACTGCCAGTGCATTGTTGGCTAGCCCAACTCGTGTAATATCACTGATTTGAGAGAAGTAAATCAGATACAC‑ACACAATAGCACACAGAACAggaacg                                                                                                                                                                          >  1:12478/1‑246 (MQ=255)
                                                                     gAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCAtt                                                                                                                                                               <  2:421760/251‑1 (MQ=255)
                                                                     gAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCAtt                                                                                                                                                               <  1:107510/251‑1 (MQ=255)
                                                                           tttGCATATCGGTAACGACAATAAACAGCAGATA‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATAAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGAt                                                                                                                                                          >  1:151670/1‑250 (MQ=255)
                                                                              gcatgacGGTAACGACAATAAACAGCAG‑TCTGTTGCGTTTGCGTCCGGTGTCCGGGTGGGCCTGGATGCTCCTCTCCTGCGGTTTTTGCGGCAAGCCCGGCCGGATGTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTACCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGGACCCC‑CAATAGCACCGACAACAGTACCGCCCTGATCATTAATATCAc                                                                                                                                                       <  2:3964/244‑1 (MQ=255)
                                                                              gCATATCGGTAACGACAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCTCTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCAc                                                                                                                                                       <  1:60700/250‑1 (MQ=255)
                                                                                 aatCGGTAACGACAATAAACAGCAGTTC‑GTTGCGTTGGCTTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGTACCGAGAACAGTACCGCACTGATCAGTAAGATCACtgtg                                                                                                                                                   <  2:381563/250‑1 (MQ=255)
                                                                                            aCAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAAttt                                                                                                                                        <  1:466371/251‑1 (MQ=255)
                                                                                               aTAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCAt                                                                                                                                     >  1:213750/1‑251 (MQ=255)
                                                                                                  aaCAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATGTCATCgg                                                                                                                                  <  1:332705/251‑1 (MQ=255)
                                                                                                            c‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGt                                                                                                                        >  1:109957/1‑251 (MQ=255)
                                                                                                                   gtttgcgtACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATg                                                                                                                  >  1:334517/1‑251 (MQ=255)
                                                                                                                                     cgATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGt                                                                                                >  2:94987/1‑251 (MQ=255)
                                                                                                                                      gATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGtc                                                                                               >  1:70117/1‑251 (MQ=255)
                                                                                                                                       aTGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGCGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGttt                                                                                              >  2:119843/1‑249 (MQ=255)
                                                                                                                                                   ccACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATc                                                                                                                                    >  1:221283/1‑201 (MQ=255)
                                                                                                                                                   ccACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATc                                                                                                                                    <  2:221283/201‑1 (MQ=255)
                                                                                                                                                    cACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGCGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGCTCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCtt                                                                                 >  2:179058/1‑251 (MQ=255)
                                                                                                                                                    cACTCTCCTGCGGTTTTTTATGCAAGCCGTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCCGTGCCTTGTTGGCTAGCACATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CACTAGCACCGAGAACAGTACCGCACTGCTCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGCATGTCTCATCAGCTTCCtt                                                                                 >  2:120164/1‑251 (MQ=255)
                                                                                                                                                            tGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGAc                                                                          >  2:112129/1‑250 (MQ=255)
                                                                                                                                                                      tCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATTCGGTAGTGCAATTAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTa                                                                <  2:70117/250‑1 (MQ=255)
                                                                                                                                                                             gCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTAt                                                        <  1:87614/251‑1 (MQ=255)
                                                                                                                                                                              cccTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACACATTTCTTAAGCATTAtc                                                       >  2:118325/1‑251 (MQ=255)
                                                                                                                                                                                  gACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTAtctc                                                     <  2:1039/249‑1 (MQ=255)
                                                                                                                                                                                    cGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTGGTTGGCTATCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAAACAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCtgat                                                 <  1:255892/251‑1 (MQ=255)
                                                                                                                                                                                            cGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCggg                                         <  2:213750/251‑1 (MQ=255)
                                                                                                                                                                                                gAGTCGGCGAGGGAACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACCGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGCACTCAATTTCATCGGACGTTTGCGGAACATGCGGTAGGGCAGTGAATGTCTCATCAGCTTCCTTGTGTGCCAAATTTCTTAATCATTATCTCTGATGAGGCggtgatt                                     >  2:255220/1‑246 (MQ=255)
                                                                                                                                                                                                          gggTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGGAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAg                           <  1:291258/251‑1 (MQ=255)
                                                                                                                                                                                                              aCGCGCCGCTGCCAGTGCCTTGTGGGCTAGCCCATCTCGGGTCATATCACTGATTTGCTAGAAGTAAAGCAGATGCACCAC‑CAATGGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACTGTTGCGTAACATGCTGTAGGGCAATGAATGTCTCATCAGCTTCATTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTaa                        <  1:35847/250‑1 (MQ=255)
                                                                                                                                                                                                                 cgcCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAAt                    <  1:162220/251‑1 (MQ=255)
                                                                                                                                                                                                                        gCCAGTGCCTTGTTGGCTAGCCCATCTCGTGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTTAAAGGGAGTAAGAATGATTggg             >  1:345581/1‑250 (MQ=255)
                                                                                                                                                                                                                           aGTGCCTTGTTGGCTAGCCCATCTCGTTTCATATCTCTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCTCTGATCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCTGTAGGTCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCtat          <  2:357608/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gttgGATAGACCATCTCGTGTCATATCACGGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGATCATTAAGATCACTGTGGTACGCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGa  <  1:347136/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                         
TCAGCGCTTTAAGGATGTCATCACCTTTAAATGGCTGACCAATACGCTGGGCTTCAGGATGCGAGTAGCGAAGACTTTGCATATCGGTAACGACAATAAACAGCAGATC‑GTTGCGTTTGCGTACGGCTTCCGCGATGGCCTGGATGCCACTCTCCTGCGGTTTTTTCTGCAAGCCCTGACGGATTTCCGGCGAGTCGGCGAGGGTACGCGCCACTGCCAGTGCCTTGTTGGCTAGCCCATCTCGCGTCATATCACTGATTTGCGAGAAGTAAATCAGATGCACCAC‑CAATAGCACCGAGAACAGTACCGCACTGACCATTAAGATCACTGTGGTACTCAATTTCATCGGACGTTTGCGTAACATGCGGTAGGGCAATGAATGTCTCATCAGCTTCCTTGTGTGACAAATTTCTTAAGCATTATCTCTGATGAGGCGGGTAATTCAAAGGGAGTAAGAATGATTGGCTATATAGGGGA  >  NC_000913/4351275‑4351751

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑