Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,251,589:1 +G 7.7% coding (253/603 nt) PP_1094 ← electron transport complex protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,251,5891.G 7.7% 251.9 / 16.0 65coding (253/603 nt)PP_1094electron transport complex protein
Reads supporting (aligned to +/- strand):  ref base . (24/36);  new base G (2/3);  total (26/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

CCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAGCG  >  NC_002947/1251442‑1251737
                                                                                                                                                     |                                                                                                                                                    
ccAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCt                                                                                                                                                     >  8:28297/1‑149 (MQ=255)
  aGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGc                                                                                                                                                   <  2:268490/149‑1 (MQ=255)
      gATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATa                                                                                                                                               >  2:296033/1‑149 (MQ=255)
        tCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGt                                                                                                                                             <  2:628/149‑1 (MQ=255)
          gCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTgg                                                                                                                                           <  7:173451/149‑1 (MQ=255)
           cccAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCGTGCCactatgg                                                                                                                                           >  6:79705/1‑143 (MQ=255)
            ccAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCa                                                                                                                                         >  2:137527/1‑149 (MQ=255)
             cAGCCCGGGCCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCAc                                                                                                                                        >  4:214155/1‑149 (MQ=255)
              aGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATa                                                                                                                                               <  5:206799/141‑1 (MQ=255)
              aGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATa                                                                                                                                               >  6:206799/1‑141 (MQ=255)
                 ccGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGTGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTa                                                                                                                                    <  4:215091/149‑1 (MQ=255)
                  tggggCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATa                                                                                                                                               <  2:217754/136‑1 (MQ=255)
                  tggggCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATa                                                                                                                                               >  1:217754/2‑137 (MQ=255)
                   ggggCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCGTGCCACTATGGCAGGTAAg                                                                                                                                  <  5:79705/149‑1 (MQ=255)
                     ggCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCGTGCCACTATGGCAGGTAAgcc                                                                                                                                <  5:290138/149‑1 (MQ=255)
                                  ttCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGc                                                                                                                   <  5:250643/149‑1 (MQ=255)
                                  ttCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGc                                                                                                                   >  6:250643/1‑149 (MQ=255)
                                     gACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTAt                                                                                                                >  1:1942/1‑149 (MQ=255)
                                        aaCCTGGCTGCCTATCAGCCGGCCGTCCGTGGTGATGGCGACGGTCAGCACGATCGTCCCGGCATACCCGTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTgg                                                                                                             <  5:218486/149‑1 (MQ=255)
                                        aaCCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTgg                                                                                                             >  3:300211/1‑149 (MQ=255)
                                         aCCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTggg                                                                                                            <  3:214155/149‑1 (MQ=255)
                                          ccTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCTACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTggg                                                                                                            >  1:70373/1‑149 (MQ=255)
                                                tGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGc                                                                                                           <  2:128070/143‑1 (MQ=255)
                                                tGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGc                                                                                                           >  1:128070/1‑143 (MQ=255)
                                                  ccGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTgc                                                                                                   >  8:91613/1‑149 (MQ=255)
                                                      tCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTgcggcg                                                                                               >  6:245912/1‑149 (MQ=255)
                                                        aGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGcggcgg                                                                                              <  4:130322/148‑1 (MQ=255)
                                                        aGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGcggcgg                                                                                              <  2:70373/148‑1 (MQ=255)
                                                          ccggccgTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGc                                                                                            >  8:113233/1‑148 (MQ=255)
                                                            ggccgTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCtt                                                                                         <  3:199156/149‑1 (MQ=255)
                                                                 tccgtGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGa                                                                                                                         >  1:25031/1‑112 (MQ=255)
                                                                 tccgtGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGa                                                                                                                         <  2:25031/112‑1 (MQ=255)
                                                                     tggtgATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCggg                                                                                <  7:91613/149‑1 (MQ=255)
                                                                     tggtgATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCggg                                                                                <  4:300211/149‑1 (MQ=255)
                                                                          aTGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTc                                                                           <  1:296033/149‑1 (MQ=255)
                                                                              cGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAgg                                                                       <  3:147964/149‑1 (MQ=255)
                                                                               gACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAggg                                                                      >  5:168377/1‑149 (MQ=255)
                                                                                 cGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGc                                                                    <  7:28297/149‑1 (MQ=255)
                                                                                   cTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTg                                                                  <  1:121649/149‑1 (MQ=255)
                                                                                     cAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTc                                                                           >  5:13597/1‑138 (MQ=255)
                                                                                     cAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTc                                                                           <  6:13597/138‑1 (MQ=255)
                                                                                      aGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTgg                                                                 <  6:167496/147‑1 (MQ=255)
                                                                                      aGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGtt                                                               <  2:1942/149‑1 (MQ=255)
                                                                                       gCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTg                                                              <  7:113233/149‑1 (MQ=255)
                                                                                        cACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCCGGATAATGGCGCTTGGCTTACCTGCCATAGGGGCACGGTAAGCCGCCAGGATACGCCTATTGGGCCTCTGTGCGGCGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTc                                                            <  4:93907/149‑1 (MQ=255)
                                                                                             aTCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATa                                                         <  5:245912/149‑1 (MQ=255)
                                                                                             aTCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATa                                                         <  5:127155/149‑1 (MQ=255)
                                                                                                   ccGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTAcc                                                   >  2:204689/1‑149 (MQ=255)
                                                                                                            cccTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCGATAGTGGCACGGTAAGCCGCCAGGATACGGGTATTGGGCCTCTGTGCGGCGGGCCTTGGCCGCGGGCTTTCCAGGGGCTGGTTGTCATGGCTACCGTCGGgcag                                          >  5:150874/1‑149 (MQ=255)
                                                                                                             ccTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTAcc                                                   <  8:128435/139‑1 (MQ=255)
                                                                                                             ccTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTAcc                                                   >  7:128435/1‑139 (MQ=255)
                                                                                                                    cTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCa                                  >  1:99688/1‑149 (MQ=255)
                                                                                                                       tgtgATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGCCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGc                                <  8:23609/148‑1 (MQ=255)
                                                                                                                           aTCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCt                            <  1:204689/148‑1 (MQ=255)
                                                                                                                            tCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTgg                          >  8:80106/1‑149 (MQ=255)
                                                                                                                                    aTGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCCAGCGGCtgg                   <  1:137527/148‑3 (MQ=255)
                                                                                                                                      ggCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTg                  >  4:162343/1‑147 (MQ=255)
                                                                                                                                       gcgcTTGGATTACCTGACATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCa                <  6:150874/148‑1 (MQ=255)
                                                                                                                                              gCTTACCGTGCCACTATGGCAGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCAc                                       >  7:34868/1‑118 (MQ=38)
                                                                                                                                              gCTTACCGTGCCACTATGGCAGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCAc                                       <  8:34868/118‑1 (MQ=38)
                                                                                                                                              gCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTgg                                                                 >  6:133700/1‑92 (MQ=255)
                                                                                                                                              gCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTgg                                                                 <  5:133700/92‑1 (MQ=255)
                                                                                                                                               cTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCtttt       <  3:162343/149‑1 (MQ=255)
                                                                                                                                                ttACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCtttt       >  7:46374/1‑148 (MQ=255)
                                                                                                                                                  aCCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCa     <  7:80106/148‑1 (MQ=255)
                                                                                                                                                  aCCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCa     <  6:168377/148‑1 (MQ=255)
                                                                                                                                                    cTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAgcg  >  7:181302/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCAGGCGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAGCG  >  NC_002947/1251442‑1251737

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: