Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,504,921 G→C 6.2% Q9H (CAG→CAC)  tusE → sulfur transfer protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,504,9210GC 6.2% 203.9 / 3.6 65Q9H (CAG→CAC) tusEsulfur transfer protein
Reads supporting (aligned to +/- strand):  ref base G (28/33);  new base C (2/2);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

GCCTGTTCGCCCTGGACGAGGATGTACAGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAAT  >  NC_002947/4504774‑4505069
                                                                                                                                                   |                                                                                                                                                    
gcCTGTTCGCCCTGGACGAGGATGTACAGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAgg                                                                                                                                                     <  6:13127/149‑1 (MQ=255)
        gCCCTGGACGAGGATGTACAGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGcc                                                                                                                                             >  8:248847/1‑149 (MQ=255)
         cccTGGACGAGGATGTACAGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGcca                                                                                                                                            >  3:67824/1‑148 (MQ=255)
          ccTGGACGAGGATGTACAGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCTTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCt                                                                                                                                           <  3:111490/149‑1 (MQ=255)
          ccTGGACGAGGATGTACAGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCt                                                                                                                                           <  2:371384/149‑1 (MQ=255)
              gACGAGGATGTACAGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGa                                                                                                                                        <  6:301985/148‑1 (MQ=255)
              gACGAGGATGTACAGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGAc                                                                                                                                       <  7:248847/149‑1 (MQ=255)
              gACGAGGATGTACAGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCACCTGATCGccaacg                                                                                                                                         <  6:368359/147‑5 (MQ=255)
                 gAGGATGTACAGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCACCTGATCGccaa                                                                                                                                           >  8:78652/1‑140 (MQ=255)
                 gAGGATGTACAGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCACCTGATCGccaa                                                                                                                                           <  7:78652/142‑3 (MQ=255)
                  aGGATGTACAGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAgg                                                                                                                                   <  4:280719/149‑1 (MQ=255)
                       gTACAGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGc                                                                                                                              >  7:190620/1‑149 (MQ=255)
                          cAGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTc                                                                                                                           <  6:32650/149‑1 (MQ=255)
                          cAGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTc                                                                                                                           <  6:38069/149‑1 (MQ=255)
                                    cgGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATc                                                                                                                 <  5:125268/149‑1 (MQ=255)
                                     ggTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCt                                                                                                                                           >  7:103838/1‑122 (MQ=255)
                                     ggTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCt                                                                                                                                           <  8:103838/122‑1 (MQ=255)
                                     ggTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCt                                                                                                                >  3:43078/1‑149 (MQ=255)
                                     ggTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCt                                                                                                                >  5:196688/1‑149 (MQ=255)
                                       tCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGGCGCTGCAATACGACAAGGGCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGc                                                                                                              <  5:247164/149‑1 (MQ=255)
                                                aCCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGt                                                                                                     <  8:190620/149‑1 (MQ=255)
                                                     gCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCac                                                                                                >  3:172487/1‑149 (MQ=255)
                                                     gCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCac                                                                                                <  2:55002/149‑1 (MQ=255)
                                                     gCAAAGGCCGTGGACTACCCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCac                                                                                                <  3:99797/149‑1 (MQ=255)
                                                       aaaGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCt                                                                                                                                           >  3:188112/1‑104 (MQ=255)
                                                       aaaGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCt                                                                                                                                           <  4:188112/104‑1 (MQ=255)
                                                           gCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGAcgcc                                                                                          >  4:117106/1‑149 (MQ=255)
                                                                ggACTACGCCGGGTTCGTCTAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGa                                                                                     <  2:189319/149‑1 (MQ=255)
                                                                 gACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAg                                                                                    >  3:218621/1‑149 (MQ=255)
                                                                      cGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCttt                                                                               <  5:94763/149‑1 (MQ=255)
                                                                         cGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTgg                                                                             >  1:264123/1‑148 (MQ=255)
                                                                         cGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGc                                                                            <  7:349000/149‑1 (MQ=255)
                                                                             ttCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAg                                                                        >  4:294551/1‑149 (MQ=255)
                                                                                gtcgAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGGAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAgcgc                                                                     >  5:319414/1‑149 (MQ=255)
                                                                                gtcgAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAgcgc                                                                     <  6:64824/149‑1 (MQ=255)
                                                                                  cgAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGAt                                                                                                                  <  8:192847/102‑1 (MQ=255)
                                                                                  cgAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGAt                                                                                                                  >  7:192847/1‑102 (MQ=255)
                                                                                   gAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGa                                                                   <  6:319414/148‑1 (MQ=255)
                                                                                    aaCTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAgg                                                                 <  2:182179/149‑1 (MQ=255)
                                                                                    aaCTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAgg                                                                 >  8:129038/1‑149 (MQ=255)
                                                                                    aaCTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAgg                                                                 <  4:172487/149‑1 (MQ=255)
                                                                                         tCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTa                                                              <  6:256977/147‑1 (MQ=255)
                                                                                           gCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATccc                                                          >  7:242457/1‑149 (MQ=255)
                                                                                           gCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATcc                                                           >  8:97917/1‑148 (MQ=255)
                                                                                           gCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATcc                                                           <  7:97917/148‑1 (MQ=255)
                                                                                              gCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGtt                                                       >  7:181655/1‑149 (MQ=255)
                                                                                                 aTACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGAc                                                    <  3:294551/149‑1 (MQ=255)
                                                                                                     gACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAg                                                                  <  1:216644/131‑1 (MQ=255)
                                                                                                     gACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAg                                                                  >  2:216644/1‑131 (MQ=255)
                                                                                                           gTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTGTGGCCACGCGCGAGGGTATCCCGTTGACGGCGGACCAc                                          >  2:75864/1‑149 (MQ=255)
                                                                                                              aaCAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGAAAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTTACGGCGGACCACTgg                                       <  2:264123/149‑1 (MQ=255)
                                                                                                                      gCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGTACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCt                               <  5:42790/149‑1 (MQ=255)
                                                                                                                        tATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTc                              >  6:4957/1‑148 (MQ=255)
                                                                                                                        tATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTc                              >  4:275486/1‑148 (MQ=255)
                                                                                                                        tATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCACCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTc                              >  6:36417/1‑148 (MQ=255)
                                                                                                                             aTACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTg                        >  2:314101/1‑149 (MQ=255)
                                                                                                                             aTACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTg                        >  5:64535/1‑149 (MQ=255)
                                                                                                                              tACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGc                       >  3:85901/1‑149 (MQ=255)
                                                                                                                                   tCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATc                                 <  3:251207/134‑1 (MQ=255)
                                                                                                                                   tCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATc                                 >  4:251207/1‑134 (MQ=255)
                                                                                                                                   tCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGt                  >  2:375609/1‑149 (MQ=255)
                                                                                                                                     aCCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTc                 <  4:218621/148‑1 (MQ=255)
                                                                                                                                     aCCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCa                <  7:129038/149‑1 (MQ=255)
                                                                                                                                       ccagggCGATCACCTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGt              >  2:229827/5‑149 (MQ=255)
                                                                                                                                       cGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGt              <  8:181655/149‑1 (MQ=255)
                                                                                                                                            gCGATCTGGTGAGCGCCCTGGACTAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTAc         <  4:67824/149‑1 (MQ=255)
                                                                                                                                                   ggTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAAt  <  5:4957/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCCTGTTCGCCCTGGACGAGGATGTACAGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAAT  >  NC_002947/4504774‑4505069

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