Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,708,280 T→C 6.5% D310G (GAT→GGT)  cysI ← sulfite reductase hemoprotein subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,708,2800TC 6.5% 179.2 / 2.6 61D310G (GAT→GGT) cysIsulfite reductase hemoprotein subunit beta
Reads supporting (aligned to +/- strand):  ref base T (29/28);  new base C (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

GTCGGCCACGGCGTCCAGCTGCTTGTCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTGCCGCCACGCAGGTGGACCATTTCGGCCTCGACCTTTTCGGCGAACACTTC  >  NC_002947/2708132‑2708428
                                                                                                                                                    |                                                                                                                                                    
gTCGGCCACGGCGTCCAGCTGCTTGTCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCAt                                                                                                                                                      >  6:109606/1‑149 (MQ=255)
 tCGGCCACGGCGTCCAGCTGCTTGTCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATc                                                                                                                                                     <  6:368804/149‑1 (MQ=255)
  cGGCCACGGCGTCCAGCTGCTTGTCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATcc                                                                                                                                                    <  5:31474/149‑1 (MQ=255)
     ccACGGCGTCCAGCTGCTTGTCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCTTCCAgg                                                                                                                                                 >  3:261707/1‑149 (MQ=255)
      cACGGCGTCCAGCTGCTTGTCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAggg                                                                                                                                                <  5:109606/149‑1 (MQ=255)
        cGGCGTCCAGCTGCTTGTCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCg                                                                                                                                              <  3:38069/149‑1 (MQ=255)
         ggCGTCCAGCTGCTTGTCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAGCCCTGGATGCCGAGCACCCAgggttc                                                                                                                                             <  6:363765/149‑4 (MQ=255)
                aGCTGCTTGTCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTc                                                                                                                                      <  5:100174/149‑1 (MQ=255)
                aGCTGCTTGTCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTc                                                                                                                                      <  6:312793/149‑1 (MQ=255)
                  ctgctTGTCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGTGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGa                                                                                                                                    <  1:302408/149‑1 (MQ=255)
                     ctTGTCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCg                                                                                                                                              <  1:82647/136‑1 (MQ=255)
                     ctTGTCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCg                                                                                                                                              >  2:82647/1‑136 (MQ=255)
                     ctTGTCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGAc                                                                                                                                   <  6:8693/147‑1 (MQ=255)
                      ttGTCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGtt                                                                                                                                >  6:282833/1‑149 (MQ=255)
                        gTCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGgtg                                                                                                                                           <  2:20094/136‑1 (MQ=255)
                        gTCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGgtg                                                                                                                                           >  1:20094/1‑136 (MQ=255)
                         tCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCAGAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCTCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTc                                                                                                                             <  4:218151/149‑1 (MQ=255)
                         tCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTc                                                                                                                             <  4:223279/149‑1 (MQ=255)
                             tCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGt                                                                                                                         <  5:282833/149‑1 (MQ=255)
                                      ccGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTAc                                                                                                                >  6:205149/1‑149 (MQ=255)
                                                  cGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGAc                                                                                                                                   <  8:141982/118‑1 (MQ=255)
                                                  cGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGAc                                                                                                                                   >  7:141982/1‑118 (MQ=255)
                                                          ttGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTgg                                                                                              <  1:310348/147‑1 (MQ=255)
                                                          ttGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTgg                                                                                              >  2:310348/1‑147 (MQ=255)
                                                          ttGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGgcg                                                                                            >  6:195968/1‑149 (MQ=255)
                                                             aGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGa                                                                                         >  1:133605/1‑149 (MQ=255)
                                                                     ggTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAg                                                                                                                          <  5:245861/108‑1 (MQ=255)
                                                                     ggTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAg                                                                                                                          >  6:245861/1‑108 (MQ=255)
                                                                     ggTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCt                                                                                 >  2:217050/1‑149 (MQ=255)
                                                                       tCAAAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTg                                                                                >  5:356210/1‑148 (MQ=255)
                                                                           aGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAg                                                                                                                    >  7:3599/1‑108 (MQ=255)
                                                                           aGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAg                                                                                                                    <  8:3599/108‑1 (MQ=255)
                                                                           aGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCActt                                                                           >  5:374846/1‑149 (MQ=255)
                                                                              cACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCActtcgg                                                                        >  7:141231/1‑149 (MQ=255)
                                                                                 gTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCt                                                                                 <  1:153592/137‑1 (MQ=255)
                                                                                 gTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCt                                                                                 >  2:153592/1‑137 (MQ=255)
                                                                                     cccGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGt                                                                 >  5:63458/1‑149 (MQ=255)
                                                                                          gcgcTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAggg                                                            <  1:344147/149‑1 (MQ=255)
                                                                                                   ggCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTgc                                                   <  6:374846/149‑1 (MQ=255)
                                                                                                      gcgGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCTTCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTgccgc                                                <  4:261707/149‑1 (MQ=255)
                                                                                                      gcgGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTgccgc                                                >  1:323967/1‑149 (MQ=255)
                                                                                                               gcgcGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGggtggtggt                                                     >  2:302093/1‑135 (MQ=255)
                                                                                                               gcgcGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGggtggtggt                                                     <  1:302093/135‑1 (MQ=255)
                                                                                                                cgcgAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTGCCGCCACGCAGGTg                                      >  7:41803/1‑149 (MQ=255)
                                                                                                                  cgAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGGAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGggtggt                                                        <  5:282125/129‑1 (MQ=255)
                                                                                                                  cgAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGGAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGggtggt                                                        >  6:282125/1‑129 (MQ=255)
                                                                                                                      cgcCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTGCCGCCACGCAGGTGGACCAt                                >  8:4549/1‑149 (MQ=255)
                                                                                                                      cgcCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTGCCGCCACGCAGGTGGACCAt                                >  4:226550/1‑149 (MQ=255)
                                                                                                                          aGCGGGCGAACCCGGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTGCCGCCACGCAGGTGGACCATTTCg                            <  5:205149/149‑1 (MQ=255)
                                                                                                                            cGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTGCCGCCACGCAGGTGGACCATTTCGGc                          <  1:217050/149‑1 (MQ=255)
                                                                                                                                  cgcccTGGATGCCGAGCACCCAGGGTGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTgc                                                   >  6:119697/3‑118 (MQ=255)
                                                                                                                                  cgcccTGGATGCCGAGCACCCAGGGTGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTgc                                                   <  5:119697/116‑1 (MQ=255)
                                                                                                                                  cgcccTGGATGCCGAGCACCCAGGGCGGTGTAGTCGACGTTGTCCTGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTGCCGCCACGCAGGTGGACCATTTCGGCCTCGAc                    >  1:268204/3‑149 (MQ=255)
                                                                                                                                  aaCCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTGCCGCCACGCAGGTGGACCATTTCg                            <  4:65047/141‑1 (MQ=255)
                                                                                                                                  aaCCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTGCCGCCACGCAGGTGGACCATTTCg                            >  3:65047/1‑141 (MQ=255)
                                                                                                                                  aaCCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTGCCGCCACGCAGGTGGACCATTTCGGCCTCGAc                    <  7:4549/149‑1 (MQ=255)
                                                                                                                                       tGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTAAGGGTGGTGGTGCCGCCACGCAGGTGGACCATTTCGGCCTCGACCtttt               >  5:212470/1‑149 (MQ=255)
                                                                                                                                             cTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTGCCGCCACGCAGGTGGACCATTTCGGCCTCGACCTTTTCGGCg          >  2:72368/1‑148 (MQ=255)
                                                                                                                                              tCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTgcc                                                  >  5:100623/1‑107 (MQ=255)
                                                                                                                                              tCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTgcc                                                  <  6:100623/107‑1 (MQ=255)
                                                                                                                                              tCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGAAGGAGTGGCCCGATACGCGCAGCACTTCGGCTTCGGTCAGGGTGGTGGTTCCGCCACGCAGGTGGACCATTTCGTCATCTACCTTTACGGCGaa        >  8:243116/1‑149 (MQ=255)
                                                                                                                                                ggCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTGCCGCCACGCAGGTGGACCATTTCGGCCTCGACCTTTTCGGCgtgtt      >  2:152941/1‑145 (MQ=255)
                                                                                                                                                    tCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTTCACTTCGGCTTCGGTCAGGGTGGTGGTGCCGCCACGCAGGTGGACCATTTCGGCCTCGACCTTTTCGGCGAACACTTc  <  6:356210/149‑1 (MQ=255)
                                                                                                                                                    tCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTGCCGCCACGCAGGTGGAACATTTCGGCCTCGACCTTTTCGGCGAACACTTc  >  6:269606/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GTCGGCCACGGCGTCCAGCTGCTTGTCGGTCACGTCGCCCGGGGCAACACCGGTAGGCTTGAGTGACAGGGTCACAGCCACGTAGCCCGGGCGCTTGTGGGCGCGGGTGTTGCGCGAGCGCCAGCGGGCGAACCCTGGGTGCTCGGCATCCAGGGCGGTGTAGTCGACGTTGTCGAGTGCGAGGTACTCGGGGTCGACGAAGTGGCGCGATACGCGCTGCACTTCGGCTTCGGTCAGGGTGGTGGTGCCGCCACGCAGGTGGACCATTTCGGCCTCGACCTTTTCGGCGAACACTTC  >  NC_002947/2708132‑2708428

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑