Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,844,871:1 +G 32.5% intergenic (+48/‑249) ccoP‑II → / → PP_4259 cbb3‑type cytochrome c oxidase subunit/iron‑sulfur cluster‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,844,8711.G32.5% 274.4 / inf 206intergenic (+48/‑249)ccoP‑II/PP_4259cbb3‑type cytochrome c oxidase subunit/iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (75/64);  new base G (37/30);  total (112/94)
Fisher's exact test for biased strand distribution p-value = 8.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.05e-03

Only 100 of 228 total aligned reads displayed.

GATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACTGCCAT  >  NC_002947/4844721‑4845020
                                                                                                                                                       |                                                                                                                                                     
gATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCt                                                                                                                                                        <  1:5599539/150‑1 (MQ=255)
  tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                                                                                                      <  1:7154051/150‑1 (MQ=255)
  tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                                                                                                      <  2:6811695/150‑1 (MQ=255)
  tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                                                                                                      <  2:7165642/150‑1 (MQ=255)
  tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                                                                                                      <  1:5840147/150‑1 (MQ=255)
  tGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCtgt                                                                                                                                                      <  2:4913918/150‑3 (MQ=255)
    ccGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGt                                                                                                                                                    <  1:319255/150‑1 (MQ=255)
    ccGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGt                                                                                                                                                    <  1:734931/150‑1 (MQ=255)
    ccGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGt                                                                                                                                                    <  2:5724806/150‑1 (MQ=255)
     cGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTa                                                                                                                                                   <  2:3221289/150‑1 (MQ=255)
      ggCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAg                                                                                                                                                  <  2:7009313/150‑1 (MQ=255)
        cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGa                                                                                                                                                <  2:2176076/150‑1 (MQ=255)
        cccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGa                                                                                                                                                <  2:9131885/150‑1 (MQ=255)
          cAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGc                                                                                                                                              >  2:1300789/1‑150 (MQ=255)
          cAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGc                                                                                                                                              >  1:9152530/1‑150 (MQ=255)
          cAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGc                                                                                                                                              >  2:5393288/1‑150 (MQ=255)
          cAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTGGTAGGAGc                                                                                                                                              <  1:3073779/150‑1 (MQ=255)
          cAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGc                                                                                                                                              <  2:265644/150‑1 (MQ=255)
           aGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCg                                                                                                                                             <  1:2958839/150‑1 (MQ=255)
            ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTTCGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                                                                                            <  2:2769115/150‑1 (MQ=255)
            ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                                            <  1:1015791/150‑1 (MQ=255)
            ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                                            <  1:8840727/150‑1 (MQ=255)
            ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCgg                                                                                                                                            <  2:6502556/150‑1 (MQ=255)
            ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                                                                                            <  2:7235474/150‑1 (MQ=255)
            ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                                                                                            <  1:2973141/150‑1 (MQ=255)
            ggAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCgg                                                                                                                                            <  2:2281818/150‑1 (MQ=255)
             gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCggg                                                                                                                                           <  1:1397128/150‑1 (MQ=255)
             gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCggg                                                                                                                                           <  2:4249533/150‑1 (MQ=255)
             gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCggg                                                                                                                                           <  2:62637/150‑1 (MQ=255)
             gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCggg                                                                                                                                           <  2:6148933/150‑1 (MQ=255)
             gAACAACACCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACACCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCAGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCggg                                                                                                                                           <  2:8232981/150‑1 (MQ=255)
               acaacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGtt                                                                                                                                         <  2:414497/150‑1 (MQ=255)
               acaacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGtt                                                                                                                                         <  2:9152530/150‑1 (MQ=255)
                caacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGttt                                                                                                                                        <  1:1177481/150‑1 (MQ=255)
                    tGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTAccc                                                                                                                                    <  2:7783721/150‑1 (MQ=255)
                    tGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTAAACAGGATCCTTGTAGGAGCGGGTTTAccc                                                                                                                                    >  2:1686648/1‑150 (MQ=255)
                        gggCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCGTGTAGGAGCGGGTTTACCCGCGa                                                                                                                                >  2:6946344/1‑150 (MQ=255)
                               gACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCtt                                                                                                                                                      >  1:4247747/1‑121 (MQ=255)
                               gACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGcc                                                                                                                         <  2:5994810/150‑1 (MQ=255)
                                      tGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCCCTCCTACACAGGATCCTTGTAGGAGCGGGGTTACCCGCGAAGAGGCCGGTAcag                                                                                                                  >  1:5988759/1‑148 (MQ=255)
                                           ttGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGcaca                                                                                                             >  1:7596947/1‑150 (MQ=255)
                                             gcgggcggcTTACGTCTACCGCCTGTCGCACCCGGCTGCGCCGGCCACGGCTGAGTACGGCCGGCCGGGCCTCTTCGCGGGTAACCCCGCTCCTCCACAGCAGCCGGGGAGGAGCGGGTTCCCCCGCGAAGAGGCCGGTACTGGCACAGc                                                                                                           <  2:7230056/147‑1 (MQ=255)
                                               tGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGGGAAACCCCGCCAAACCCCGATACGTGTTGGAGCGGGTTTTACCCCGCAGAGGGCGGTACAGGGCCAGCAt                                                                                                         >  1:2643275/1‑150 (MQ=255)
                                                 gcggcTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACCGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCAttt                                                                                                       <  2:7396843/150‑1 (MQ=255)
                                                 gcggcTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAAACCGCTCCTACACAGGATCCGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCAttt                                                                                                       >  2:4984301/1‑150 (MQ=255)
                                                       tACGTCTACAGCCTGTCGCACCAGGCCGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACCCAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCa                                                                                                 <  2:4520690/150‑1 (MQ=255)
                                                            cTACCGCCTGCCCCACCCGGCTGCGCCGGCCACGGCTGAGTACGGCCGGCCTGGCCTCTGCGGGGGTAAACCCGCGCCTACCCCGTTTCCTGGAGGGCGGGGTTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCaa                                                                                            <  2:1162801/150‑1 (MQ=255)
                                                                   cTGTCCCACCCGGCTGCCCCGCCCACCGCTGAGTAAGCCAGGCCTGCCCACTACGCGGTAAACCCCGCCCCTCCACACCATCCTTGTTGGAGGGGGTAAACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCt                                                                                     <  1:746091/150‑1 (MQ=255)
                                                                      tCGCACCAGGCTGAGCCGGCCACGGCTGAGTAAGGCAGGCCGGGCCTCTTCGCGGTTAAACCCGCTCCTACACAGGATCCTTGTAGGTGGGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGa                                                                                  <  2:8488022/150‑1 (MQ=38)
                                                                            cAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGAACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                            <  2:1464333/150‑1 (MQ=255)
                                                                              ggCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTTAACCCCCACCGACACACGAAACCTGGAGGCGCGGGTTTAACCGCCGACAGGCCCGTACAGGCACAGCAATTCTAGCAGGCAACTGGCCAGGTTCAAATGa                                                                          >  1:7756235/1‑150 (MQ=255)
                                                                               gCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAAACCGCTCCTACACAGGATCCGTGTAGGAGCGGGTTTAACCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATAAATTGac                                                                         >  1:6744911/1‑150 (MQ=255)
                                                                                      cGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAAACCGCTCCAACACAGGATACGTTTAGGATCGGGTTTAACCCCGAAGAGGCCGGTAATGGCACAGCATTTATCGCAGGCAAGTTTCCTGGATCAATTGACACCCGcc                                                                  >  2:1793169/1‑150 (MQ=37)
                                                                                                       ggcaggcCTGGCCTCTTCGCGGGTAAACCCGCTCCAACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGGACTGGCACAGCATTTCTCGCAGGCAAGTGACCAGGATCAAATTACACCCGCCATAAAACGTTTCATAcc                                                 >  2:3129832/1‑150 (MQ=37)
                                                                                                        gcaggcCTGGCCTCTTCGCGGGTAAAAACGCTCCGACACAGGATACGTGTTGGAGCGGGTTTTACCGCGCAGAGGCCGGTACTTGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCCCCATAACACGTTTCATACCa                                                >  2:264641/1‑150 (MQ=255)
                                                                                                          aggcCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACg                                              >  2:4341474/1‑150 (MQ=255)
                                                                                                            gcctggcctCTCCGCGGGTAAACCCGCTCCTACACCGGATCCTTGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGaa                                            <  2:7802104/150‑1 (MQ=255)
                                                                                                              ctggcctCTTCGCGGGTAAACCCGCTCCTACACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCAATTCTAGCAGGCAAGTGACCAGGATCAATTGTCACCCGCCATAAAACGTTTCAAACCACGAAcc                                          >  1:6093903/1‑150 (MQ=25)
                                                                                                               tggcctCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAAccc                                         >  2:8288499/1‑150 (MQ=255)
                                                                                                                 gcctCTCCGCGGGAAACCCCGCTCCTACACAGGATCCTTTTTGGAGGGGTTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCaa                                       <  1:8002604/150‑1 (MQ=39)
                                                                                                                    tctTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGCAGAGGCCGGTACTGGCACAGCATTTCTCGCCGGCAAGTGACCGGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAAccc                                    >  2:5840172/1‑148 (MQ=255)
                                                                                                                    tctTCGCGGGTAAACCCCCTCCGACACACGATCCGTGTAGGAGCGGGTTTACCCCCGAAGAGGCCCGTACAGGCACACCAATTCTCGCAGGCCAGTGACCTGGATCAATTGACACCCGCCAAAACACCGTTCATAACACGAAACCAAagc                                    >  2:2483854/1‑150 (MQ=255)
                                                                                                                       tCGCGGGAAAACCCGCTCCTACACAGTATCCTTGTTGGCGCGGTTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGAc                                 <  1:9029188/150‑1 (MQ=39)
                                                                                                                         gcgGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTa                               >  1:3468114/1‑150 (MQ=255)
                                                                                                                           gggTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTaaa                             >  2:4463696/1‑150 (MQ=255)
                                                                                                                            gtaaaaCCCGCTCATAAACTTTATCTTTGTTGGCGGGTTTTAACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAg                            <  1:326519/147‑1 (MQ=255)
                                                                                                                            ggTAAACCCGCTCCTACACAGGATCCTTGGAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGACCCCAACGCGACTAAAg                            >  2:7872733/1‑150 (MQ=255)
                                                                                                                            ggTAAACCCGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATCACACGTTTCATACCACGAACCCAACGCGACTAAAg                            >  2:7684151/1‑150 (MQ=255)
                                                                                                                             gTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAgg                           >  1:656402/1‑150 (MQ=255)
                                                                                                                              tacaCCCGCGCCTACCCCGGTTCGTTGTTGGTGGGGTTTCCCCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGt                          <  1:1132555/147‑1 (MQ=255)
                                                                                                                              tAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGt                          >  1:2676428/1‑150 (MQ=255)
                                                                                                                               aaaCCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTc                         >  1:2428591/1‑150 (MQ=255)
                                                                                                                                 aCCCGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGAt                                                                                 >  2:5437712/1‑92 (MQ=21)
                                                                                                                                  cccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCa                      >  1:5256270/1‑150 (MQ=255)
                                                                                                                                  cccGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCa                      >  2:2232826/1‑150 (MQ=255)
                                                                                                                                   ccGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg                     >  2:2012365/1‑150 (MQ=255)
                                                                                                                                   ccGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg                     >  1:233373/1‑150 (MQ=255)
                                                                                                                                   ccGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg                     >  1:6164492/1‑150 (MQ=255)
                                                                                                                                   ccGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAg                     >  1:1851325/1‑150 (MQ=255)
                                                                                                                                    cGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGt                    <  2:4947086/150‑1 (MQ=255)
                                                                                                                                    cGCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGt                    >  2:4305876/1‑150 (MQ=255)
                                                                                                                                     gCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTg                   >  2:7948487/1‑150 (MQ=255)
                                                                                                                                     gCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTg                   >  1:4645354/1‑150 (MQ=255)
                                                                                                                                     gCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTg                   >  1:4979067/1‑150 (MQ=255)
                                                                                                                                     gCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTg                   >  2:3040202/1‑150 (MQ=255)
                                                                                                                                     gCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTg                   >  2:3793921/1‑150 (MQ=255)
                                                                                                                                     gCTCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTg                   >  2:3786668/1‑150 (MQ=255)
                                                                                                                                     gCTCCTACCAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTg                   >  2:3418547/1‑150 (MQ=255)
                                                                                                                                       tCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCg                 >  1:3349775/1‑150 (MQ=255)
                                                                                                                                       tCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCg                 >  2:1073257/1‑150 (MQ=255)
                                                                                                                                       tCCTACAAGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCg                 >  2:4017761/1‑150 (MQ=255)
                                                                                                                                         cTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGaa               >  2:7548325/1‑150 (MQ=255)
                                                                                                                                           tacaaGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAAc              >  2:3609061/5‑150 (MQ=255)
                                                                                                                                             acaGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCa           >  1:4514763/1‑150 (MQ=255)
                                                                                                                                             acaGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCa           >  1:3913516/1‑150 (MQ=255)
                                                                                                                                              caGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCAt          >  2:8769714/1‑150 (MQ=255)
                                                                                                                                              aaGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCa           >  2:743293/2‑150 (MQ=255)
                                                                                                                                                   tCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACTGc     >  2:5112040/1‑150 (MQ=255)
                                                                                                                                                      ttGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACTGCCAt  >  1:2549667/1‑150 (MQ=255)
                                                                                                                                                      ttGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACTGCCAt  >  1:1471927/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GATGCCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGATCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACGCGACTAAAGGTCGCAGTGCGAACCCATACTGCCAT  >  NC_002947/4844721‑4845020

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: