Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,814,593 C→A 5.3% intergenic (+93/‑20) frmC → / → ispF S‑formylglutathione hydrolase/S‑lactoylglutathione hydrolase/2‑C‑methyl‑D‑erythritol 2,4‑cyclodiphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,814,5930CA 5.3% 245.6 / 2.9 76intergenic (+93/‑20)frmC/ispFS‑formylglutathione hydrolase/S‑lactoylglutathione hydrolase/2‑C‑methyl‑D‑erythritol 2,4‑cyclodiphosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (33/39);  new base A (2/2);  total (35/41)
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.81e-01

CCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGACGGTGACGTGCTGCTGCACGCCC  >  NC_002947/1814445‑1814742
                                                                                                                                                    |                                                                                                                                                           
ccAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGAc                                                                                                                                                             >  7:188276/1‑149 (MQ=255)
    cTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTTCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTc                                                                                                                                                           <  7:238469/147‑1 (MQ=255)
    cTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTc                                                                                                                                                           >  8:238469/1‑147 (MQ=255)
        aTCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCa                                                                                                                                                      <  3:183587/148‑1 (MQ=255)
         tCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAgg                                                                                                                                                    <  7:261671/149‑1 (MQ=255)
         tCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAgg                                                                                                                                                    <  1:30712/149‑1 (MQ=255)
         tCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAgg                                                                                                                                                    <  1:48004/149‑1 (MQ=255)
         tCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGAATGAGTCAgg                                                                                                                                                    >  3:140390/1‑149 (MQ=255)
         tCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGAATGAGTCAgg                                                                                                                                                    >  1:63234/1‑149 (MQ=255)
           gaggagCACCTGCGCCATCGTGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGc                                                                                                                                                  >  5:21143/1‑149 (MQ=255)
           gaggagCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGAATGAGTCAGGGc                                                                                                                                                  <  2:63234/149‑1 (MQ=255)
            aggagCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGAATGAGTCAGGGCg                                                                                                                                                 <  4:140390/149‑1 (MQ=255)
              gagCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGtt                                                                                                                                               >  4:65212/1‑149 (MQ=255)
                 cACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTc                                                                                                                                            >  4:344462/1‑149 (MQ=255)
                 cACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTc                                                                                                                                            <  6:65063/149‑1 (MQ=255)
                 cACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTc                                                                                                                                            <  2:364159/149‑1 (MQ=255)
                  cccTGCGCCATCATGCGGTGGCACGGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCt                                                                                                                                           <  7:291660/148‑1 (MQ=255)
                   ccTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCtt                                                                                                                                          <  1:164164/149‑1 (MQ=255)
                       cgcCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATg                                                                                                                                       <  5:172586/148‑1 (MQ=255)
                          catcatGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTa                                                                                                                                   <  5:83619/149‑1 (MQ=255)
                          catcatGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTa                                                                                                                                   >  3:51153/1‑149 (MQ=255)
                                  ggTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGcc                                                                                                                             >  8:175947/1‑147 (MQ=255)
                                  ggTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCAc                                                                                                                           >  1:342396/1‑149 (MQ=255)
                                    tGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACgg                                                                                                                         <  4:181110/149‑1 (MQ=255)
                                      gCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCt                                                                                                                       <  7:222483/149‑1 (MQ=255)
                                        aCTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAgg                                                                                                                                                    >  7:160763/1‑118 (MQ=255)
                                        aCTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAgg                                                                                                                                                    <  8:160763/118‑1 (MQ=255)
                                        aCTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTAc                                                                                                                     >  5:80529/1‑149 (MQ=255)
                                               cGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGc                                                                                                              <  7:175947/149‑1 (MQ=255)
                                                   tttAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCg                                                                                                          >  5:293274/1‑149 (MQ=255)
                                                      aaCGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCa                                                                                                             <  5:55416/143‑1 (MQ=255)
                                                      aaCGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCa                                                                                                             >  6:55416/1‑143 (MQ=255)
                                                               tgtgGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACgg                                                                                         >  8:345494/1‑149 (MQ=255)
                                                                 tgGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGt                                                                                        <  2:67116/148‑1 (MQ=255)
                                                                 tgGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGt                                                                                        <  3:65212/148‑1 (MQ=255)
                                                                 tgGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTg                                                                                       <  6:21143/149‑1 (MQ=255)
                                                                      gCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTc                                                                                  <  2:342396/149‑1 (MQ=255)
                                                                              gCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCt                                                                          >  2:225930/1‑149 (MQ=255)
                                                                                ccGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTgg                                                                        >  2:18396/1‑149 (MQ=255)
                                                                                ccGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTgg                                                                        >  5:51866/1‑149 (MQ=255)
                                                                                    ttggggCTGCTGTGCAGCCCTAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCgg                                                                   >  3:341527/4‑149 (MQ=255)
                                                                                        ggCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTg                                                               >  1:369053/1‑149 (MQ=255)
                                                                                          cTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCg                                                             <  6:80529/149‑1 (MQ=255)
                                                                                            gCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTa                                                           >  2:102276/1‑149 (MQ=255)
                                                                                             cacaGCAGGCGCAATCGCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCggc                                                                  <  8:88051/141‑2 (MQ=255)
                                                                                             cacaGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTAt                                                          >  5:33758/1‑149 (MQ=255)
                                                                                             cacaGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTCACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCCCGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCggc                                                                  >  7:88051/1‑140 (MQ=255)
                                                                                             cacaGCAGCCCCAATCCCACAAATCCACTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTAt                                                          >  6:191071/1‑149 (MQ=255)
                                                                                             cacaGCAGCCCCAATCCCACAAATCCACTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTAt                                                          >  6:185710/1‑149 (MQ=255)
                                                                                               cagcagCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATcc                                                        >  2:113479/1‑149 (MQ=255)
                                                                                                 gcagcCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATcccc                                                      <  4:51153/149‑1 (MQ=255)
                                                                                                     ccccAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCcaca                                                  <  1:102276/149‑1 (MQ=255)
                                                                                                     ccccAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCcac                                                   >  6:268272/1‑148 (MQ=255)
                                                                                                      cccAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCcaca                                                  >  8:23499/1‑148 (MQ=255)
                                                                                                       ccAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACaaa                                                >  4:284053/1‑149 (MQ=255)
                                                                                                       ccAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACaaa                                                >  8:258430/1‑149 (MQ=255)
                                                                                                             ccACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGCGCGGGGTGCGTATCCCCCACAAATACg                                            <  1:209120/148‑1 (MQ=255)
                                                                                                              cacaAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGcc                                         >  4:365008/1‑149 (MQ=255)
                                                                                                               acaAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCg                                                                    >  1:69165/1‑121 (MQ=255)
                                                                                                               acaAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCg                                                                    <  2:69165/121‑1 (MQ=255)
                                                                                                               acaAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCct                                        >  3:192490/1‑149 (MQ=255)
                                                                                                                  aaTCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCctgct                                     <  3:284053/149‑1 (MQ=255)
                                                                                                                       ccTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGccc                                <  7:345494/149‑1 (MQ=255)
                                                                                                                           cAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCAcc                                                                                                           <  4:170225/76‑1 (MQ=255)
                                                                                                                               gCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCCTGGGGGTTTTCTGAGACGGGTGATTTAATTACCATGGGCGGGGTGCTTATCCCCAACAAATACGCCCTGCTGGCCCACt                             >  1:307984/1‑149 (MQ=255)
                                                                                                                               gCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGAc                        <  7:258430/149‑1 (MQ=255)
                                                                                                                                cAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGACg                       <  5:268272/149‑1 (MQ=255)
                                                                                                                                cAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGACg                       <  1:225930/149‑1 (MQ=255)
                                                                                                                                cAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGACg                       <  5:185710/149‑1 (MQ=255)
                                                                                                                                cAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGACg                       <  5:191071/149‑1 (MQ=255)
                                                                                                                                   gTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGACGGTg                    >  6:45516/1‑149 (MQ=255)
                                                                                                                                      gAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGACGGTGAc                  <  6:293274/148‑1 (MQ=255)
                                                                                                                                         tCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGcc                                 <  7:214586/130‑1 (MQ=255)
                                                                                                                                         tCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGcc                                 >  8:214586/1‑130 (MQ=255)
                                                                                                                                         tCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGACGGTGACGtgc              <  1:18396/149‑1 (MQ=255)
                                                                                                                                           aCGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGACGGTGACGtgctg            <  5:74692/149‑1 (MQ=255)
                                                                                                                                               ccTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGACGGTGACGtgctgctgc        <  3:365008/149‑1 (MQ=255)
                                                                                                                                                    cTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTAACCCCCACAAATACGGCCTGCTGGCCCACTCCGACGGTGAAGTGCTCCAGCACGcc   >  5:225906/1‑149 (MQ=255)
                                                                                                                                                    cTCATTCAGGGCGTTTTCTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGACGGTGACGTGCAGCAGCACgtca  <  8:188276/149‑4 (MQ=255)
                                                                                                                                                    |                                                                                                                                                           
CCAGCTTCATCGAGGAGCACCTGCGCCATCATGCGGTGGCACTGGGGCGGGTTTAACGCCTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCACAGCAGCCCCAATCCCACAAATCCCCTGCAAAGCAGGTAGAATCACGCCCTGACTCATTCAGGGCGTTTTCTTATGCGTATTGGCCACGGCTACGATGTGCACCG‑‑‑‑TTTCTGCGACGGTGATTTCATTACCCTGGGCGGGGTGCGTATCCCCCACAAATACGGCCTGCTGGCCCACTCCGACGGTGACGTGCTGCTGCACGCCC  >  NC_002947/1814445‑1814742

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: