Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 660,909 T→C 14.0% intergenic (+41/‑126) speA → / → PP_0568 biosynthetic arginine decarboxylase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947660,9090TC14.0% 137.2 / 13.6 57intergenic (+41/‑126)speA/PP_0568biosynthetic arginine decarboxylase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/30);  new base C (4/4);  total (23/34)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

CATTACCGCGACAAGGTTGCCAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTC  >  NC_002947/660761‑661051
                                                                                                                                                    |                                                                                                                                              
cATTACCGCGACAAGGTTGCCAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTtgt                                                                                                                                                <  6:236327/149‑1 (MQ=255)
 aTTACCGCGACAAGGTTGCCAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTtgtg                                                                                                                                               >  8:64434/1‑149 (MQ=255)
  ttACCGCGACAAGGTTGCCAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTgg                                                                                                                                              <  4:228317/149‑1 (MQ=255)
     ccGCGACAAGGTTGCCAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTTGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAg                                                                                                                                           <  1:155397/149‑1 (MQ=255)
     ccGCGACAAGGTTGCCAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGa                                                                                                                                            <  1:199230/148‑1 (MQ=255)
     ccGCGACAAGGTTGCCAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGa                                                                                                                                            <  1:233037/148‑1 (MQ=255)
      cgcgACAAGGTTGCCAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGc                                                                                                                                          >  1:99570/1‑149 (MQ=255)
        cgACAAGGTTGACAGCGCCAAGAGAACTGCGCGTCAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGTTCTTCGTACCTGTAGTCGTAATTGACGGGGACGCTCCCACAAAAGGCAGCACAGGCTTTTGCAGGAGCgg                                                                                                                                        <  7:128844/149‑1 (MQ=25)
            aaGGTTGCCAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGGCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCtt                                                                                                                                    <  4:323190/149‑1 (MQ=255)
              ggTTGCCAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCAGTGCTGCCTTTTGTGGGAGCGGccaagc                                                                                                                                  >  2:308936/1‑145 (MQ=255)
               gTTGCCAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCg                                                                                                                                 >  6:160728/1‑149 (MQ=255)
                   ccAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTcgc                                                                                                                             <  1:46602/149‑1 (MQ=255)
                     aGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGa                                                                                                                           <  6:24518/149‑1 (MQ=255)
                     aGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGa                                                                                                                           <  6:29794/149‑1 (MQ=255)
                           aaGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAg                                                                                                                     <  8:100681/149‑1 (MQ=255)
                             gATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAggg                                                                                                                   >  5:98381/1‑149 (MQ=255)
                                  cTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAg                                                                                                              <  6:294674/149‑1 (MQ=255)
                                      gcgTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCt                                                                                                          <  7:64434/149‑1 (MQ=255)
                                       cgTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTGGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCt                                                                                                          >  5:311329/1‑148 (MQ=255)
                                         tGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCAc                                                                                                       >  4:28493/1‑149 (MQ=255)
                                                        tCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCg                                                                                         <  8:296263/148‑1 (MQ=255)
                                                            gATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTaaa                                                                                    >  7:78476/1‑149 (MQ=255)
                                                            gATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTaa                                                                                     <  1:1452/148‑1 (MQ=255)
                                                                ccTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAatat                                                                                 <  6:311329/148‑1 (MQ=255)
                                                                  ttGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCg                                                                              >  3:344143/1‑149 (MQ=255)
                                                                    gcgcCTGGGCTTGACCCGGTCTTCGTACGTGTGGTCGTTATTGTCGGGGACGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCgg                                                                             <  1:308936/148‑1 (MQ=255)
                                                                    gcgcCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCgg                                                                             <  3:28493/148‑1 (MQ=255)
                                                                      gcCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGccc                                                                          >  3:267547/1‑149 (MQ=255)
                                                                       ccTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGt                                                                                                 >  1:33048/1‑125 (MQ=255)
                                                                       ccTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGt                                                                                                 <  2:33048/125‑1 (MQ=255)
                                                                        cTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTg                                                                        >  5:132812/1‑149 (MQ=255)
                                                                           ggCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTgg                                                                       <  8:78476/147‑1 (MQ=255)
                                                                                 aCCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCg                                                                <  6:132812/148‑1 (MQ=255)
                                                                                  cccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCa                                                              >  3:270325/1‑149 (MQ=255)
                                                                                      gTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCaa                                                          <  2:99570/149‑1 (MQ=255)
                                                                                            gTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATg                                                        >  7:115320/1‑145 (MQ=255)
                                                                                            gTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATg                                                        <  8:115320/145‑1 (MQ=255)
                                                                                                cTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTAcc                                                >  1:238991/1‑149 (MQ=255)
                                                                                                     gtcgtAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGACGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCCACGCACCAATTCCGCTTccgggg                                            >  7:234412/1‑144 (MQ=255)
                                                                                                        gtAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATg                                        <  2:238991/149‑1 (MQ=255)
                                                                                                            ttGACGGGGCCGCTCCTGCAAAAGCCTGTGCTTGCTTTTGTGGGTGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAt                                    <  6:98381/149‑1 (MQ=255)
                                                                                                             acgcaaggCCGCTCCCACAAAAGGCAGCACAGGCTTTTGCAGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAAt                                    <  7:37405/142‑1 (MQ=25)
                                                                                                                cGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTc                                <  8:167396/149‑1 (MQ=255)
                                                                                                                 aaggCCGCTCCCACAAAAGGCAGCACAGGCTTTTGCAGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCttt                                             >  4:139687/3‑135 (MQ=25)
                                                                                                                 aaggCCGCTCCCACAAAAGGCAGCACAGGCTTTTGCAGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCttt                                             <  3:139687/133‑1 (MQ=25)
                                                                                                                   ggCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCaaa                             <  5:160728/149‑1 (MQ=255)
                                                                                                                    gCCGCTCCCACAAAAGGCAGCACAGGCTTTTGCAGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCCCCAATGCCGCTACCTTTGAATGCAATGTTCAAag                            >  3:305204/1‑149 (MQ=37)
                                                                                                                    gCCGCTCCCACAAAAGGCAGCACAGGCTTTTGCAGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCt                                               <  8:246258/130‑1 (MQ=25)
                                                                                                                    gCCGCTCCCACAAAAGGCAGCACAGGCTTTTGCAGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCt                                               >  7:246258/1‑130 (MQ=25)
                                                                                                                    gCCGCTCCCACAAAAGGCAGCACAGGCTTTTGCAGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAag                            >  5:244220/1‑149 (MQ=37)
                                                                                                                     ccGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAagg                           <  4:267547/149‑1 (MQ=255)
                                                                                                                       gCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAaggag                         >  7:171278/1‑149 (MQ=255)
                                                                                                                          ccaacAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGAc                      <  1:303464/145‑1 (MQ=255)
                                                                                                                               aaaaGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGatat                 <  4:344143/149‑1 (MQ=255)
                                                                                                                                           tgcCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAgg     >  3:153679/1‑149 (MQ=255)
                                                                                                                                            gcCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAAt                                         >  2:292630/1‑112 (MQ=255)
                                                                                                                                            gcCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAAt                                         <  1:292630/112‑1 (MQ=255)
                                                                                                                                             ccTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAgg     <  4:270325/147‑1 (MQ=255)
                                                                                                                                               ttttGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTc  >  3:8958/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CATTACCGCGACAAGGTTGCCAGCGCCAAGATCACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTC  >  NC_002947/660761‑661051

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