Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 353,067 A→G intergenic (+34/+84) hisF → / ← cbcV imidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947353,0670AG98.4% 219.6 / NA 63intergenic (+34/+84)hisF/cbcVimidazole glycerol phosphate synthase subunit HisF/choline/betaine/carnitine ABC transporter ATP binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (28/34);  minor base C (1/0);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAG  >  NC_002947/352923‑353201
                                                                                                                                                |                                                                                                                                       
cTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGCTCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCggg                                                                                                                                       >  6:46230/1‑147 (MQ=255)
cTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCggg                                                                                                                                       <  8:49649/147‑1 (MQ=255)
 tGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCggg                                                                                                                                       <  6:302161/146‑1 (MQ=255)
  ggAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCAc                                                                                                                                    >  7:298832/1‑148 (MQ=255)
  ggAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCAc                                                                                                                                    >  7:347932/1‑148 (MQ=255)
   gAAGGTCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGATTTGCAGCCCTTCGCGGGCa                                                                                                                                     <  8:97624/146‑1 (MQ=255)
   gAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACg                                                                                                                                   <  1:23012/148‑1 (MQ=255)
          aTGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCt                                                                                                                             <  5:398833/147‑1 (MQ=255)
             ccAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCca                                                                                                                         <  1:286339/148‑1 (MQ=255)
             ccAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCca                                                                                                                         >  7:394268/1‑148 (MQ=255)
              cAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCca                                                                                                                         <  6:290216/147‑1 (MQ=255)
                    ggTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCCGCCCTTCGCGGGCACGCCCGCTCCCACAGGtt                                                                                                                   >  6:402947/1‑147 (MQ=255)
                     gTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGttt                                                                                                                  <  4:190201/147‑1 (MQ=255)
                      tGCTGGCGGCCAGCATCTTCCACTTTGGCTAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTAGTTCGCTTATCCATTGACCGGGTCTGCTTTGCAGCCATTCGCGGGCACGCCCGCTCCCACAGGTTTGc                                                                                                                <  7:350448/148‑1 (MQ=39)
                        cTGGCGGCCAGCATCTTCCACTTTGGCGATTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                               <  4:347737/147‑1 (MQ=255)
                         tGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGc                                                                                                                <  8:248097/145‑1 (MQ=255)
                         tGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGc                                                                                                                >  7:248097/1‑145 (MQ=255)
                           gcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCt                                                                                                                             <  7:354551/130‑1 (MQ=255)
                           gcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCt                                                                                                                             >  8:354551/1‑130 (MQ=255)
                           gcggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCGCTCCCACAGGGTTGGTggga                                                                                                          >  2:145981/1‑148 (MQ=39)
                            cggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCt                                                                                                               >  7:205184/1‑143 (MQ=255)
                            cggcCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTgag                                                                                                           >  3:291710/1‑146 (MQ=255)
                             ggcCAGCATCTACCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCCGGGCTGCTTTGCAGCCCTTAGCCGGCACGCCCGCTACCACAGGGTAGCTgagaa                                                                                                         >  1:324161/1‑146 (MQ=255)
                               ccAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGAtctc                                                                                                      <  1:215675/148‑1 (MQ=255)
                                cAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATctct                                                                                                     >  2:294454/1‑148 (MQ=255)
                                 aGCATCTTCCACTTTGGCGAGTACACTGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCTCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtg                                                                                                    <  1:301703/148‑1 (MQ=255)
                                 aGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtg                                                                                                    <  1:378083/148‑1 (MQ=255)
                                  gCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCtg                                                                                                    <  4:291710/147‑1 (MQ=255)
                                        tCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTGCGGAGATCTCTGTGGGAGc                                                                                             >  7:373267/1‑148 (MQ=255)
                                           aCTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCgg                                                                                           <  4:208723/147‑1 (MQ=255)
                                            cTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGGCAGCCCGCACCCACAGGGTTGCTGCGATCTCTGTGGGAGCggg                                                                                          >  6:252296/1‑147 (MQ=255)
                                             tttGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCg                                                                                        >  8:72221/1‑148 (MQ=255)
                                                ggCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                            >  1:196089/1‑125 (MQ=255)
                                                ggCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTga                                                                                                            <  2:196089/125‑1 (MQ=255)
                                                  cGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGcc                                                                                    >  3:305954/1‑147 (MQ=255)
                                                         acaGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCGCTGTGGGAGCGGGCGTGCCCGCGAAg                                                                             >  3:355834/1‑147 (MQ=255)
                                                                  gAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGGCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCCGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTCGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGaa                                                                    <  1:294454/148‑1 (MQ=38)
                                                                   aGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGCAg                                                                   >  3:381619/1‑147 (MQ=255)
                                                                   aGGCAAAGGGCTACATGGCTTGGCGCGGGATCGTCGTTCGCTGATCCATTTACGGGGGCTGTTTTGCAGCCCTTGGCGGGCACGCCGGCTCCCACAGTTTTGCTTAGATCTCTGTGGGAGGGGGCGTGCCCGCGAAGGGCCGCGAAg                                                                   <  4:355834/147‑1 (MQ=255)
                                                                      ccAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCgg                                                                <  4:305954/147‑1 (MQ=255)
                                                                        aaGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGAtc                                                                                                        >  1:11145/1‑105 (MQ=255)
                                                                        aaGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGAtc                                                                                                        <  2:11145/105‑1 (MQ=255)
                                                                        aaGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGa                                                                     <  3:135555/140‑1 (MQ=255)
                                                                        aaGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGGTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGa                                                                     >  4:135555/1‑140 (MQ=255)
                                                                        aaGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTACCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCgggc                                                              >  4:308167/1‑145 (MQ=255)
                                                                           gCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACGGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCa                                                           <  4:114413/147‑1 (MQ=255)
                                                                           gCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCCCTCCCACACGTTTGCGGAGCTCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGCAGCGGCCCCa                                                           >  4:287463/1‑147 (MQ=37)
                                                                                 aTGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGc                                                      <  6:405625/146‑1 (MQ=255)
                                                                                     cTTCGCGCGGGATCGGCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCa                                                  <  3:139359/146‑1 (MQ=255)
                                                                                             gggATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCGGAGATCTCTGTGGGAGCGGGCGGGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTGAACCGTGGTTc                                        >  7:174677/1‑148 (MQ=255)
                                                                                               gATGGTCGTTCGCGTATCCATTGACGGGGGCTGCTTTGCAGCCCTGCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGGGGGGGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCtt                                      <  8:174677/148‑1 (MQ=38)
                                                                                                   gtcgtTCGCTGATCAATTGACGGGGGCTGCTTTGCAGCCCTTCGGGGGCACGCCCGCTCCCACAGGTTTGCTGATATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGc                                    <  3:229581/146‑1 (MQ=38)
                                                                                                        tCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCt                                       >  3:41975/1‑138 (MQ=255)
                                                                                                        tCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCt                                       <  4:41975/138‑1 (MQ=255)
                                                                                                             gATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTgag                                                                                                           <  3:53713/65‑1 (MQ=255)
                                                                                                             gATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTgag                                                                                                           >  4:53713/1‑65 (MQ=255)
                                                                                                               tCCATTGACCGGGGCGGCTTTCCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGAGCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGt                        <  5:252296/146‑1 (MQ=255)
                                                                                                                 cATTGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTa                    <  8:394268/148‑1 (MQ=255)
                                                                                                                   ttGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCg                                                                                            >  3:66306/1‑74 (MQ=255)
                                                                                                                   ttGACCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCg                                                                                            <  4:66306/74‑1 (MQ=255)
                                                                                                                       ccGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGc                >  7:91535/1‑146 (MQ=255)
                                                                                                                           ggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCt               <  6:409517/143‑1 (MQ=255)
                                                                                                                           ggCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCt               >  5:409517/1‑143 (MQ=255)
                                                                                                                                   tGCCGCCCGGCGCGGGCACCCCCGCCCCCCCCGGTTTGCTGACATCGCGGGGGGCGCGGGGGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCg    <  8:347932/146‑1 (MQ=37)
                                                                                                                                    gCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAg  <  7:72221/147‑1 (MQ=255)
                                                                                                                                            tCGCGGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGTCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGagc                               >  3:37053/1‑110 (MQ=255)
                                                                                                                                                |                                                                                                                                       
CTGGAAGGGCATGCCAGCGCGGTGCTGGCGGCCAGCATCTTCCACTTTGGCGAGTACACAGTGCCGGAGGCCAAGGCCTACATGGCTTCGCGCGGGATCGTCGTTCGCTGATCCATTGACCGGGGCTGCTTTGCAGCCCTTCGCAGGCACGCCCGCTCCCACAGGTTTGCTGAGATCTCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCGAAGCGGCCCCAATTGCTTCAACCGTGGTTCTTGCCGAGCAGCGCGTGGTAAAGCTCGGTATCCCCGAG  >  NC_002947/352923‑353201

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: