Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,708,586 C→T 6.1% P222P (CCG→CCA 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,5860CT 6.1% 242.3 / 3.4 81P222P (CCG→CCAcysIsulfite reductase hemoprotein subunit beta
Reads supporting (aligned to +/- strand):  ref base C (38/38);  new base T (3/2);  total (41/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CATTTCGGCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTCCTGCGA  >  NC_002947/2708439‑2708735
                                                                                                                                                   |                                                                                                                                                     
cATTTCGGCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTGCTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACgg                                                                                                                                                     <  1:548806/150‑1 (MQ=255)
cATTTCGGCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACgg                                                                                                                                                     >  1:249114/1‑150 (MQ=255)
  tttCGGCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTTTTGGCGTCCGTAACGTTTGTACACGCCCAGGATGTCGTCGAGGTAGCTGATCAGGTCCTGCCGCGTCa                                                                                                                                                   >  2:1521385/1‑150 (MQ=255)
  tttCGGCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCa                                                                                                                                                   <  1:1782831/150‑1 (MQ=255)
    tCGGCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAgg                                                                                                                                                 <  2:1436246/150‑1 (MQ=255)
         cTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTc                                                                                                                                            >  2:1942459/1‑150 (MQ=255)
          tCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGACCAGGTCCTGCCACGGCAGGAACTCg                                                                                                                                           >  2:151754/1‑150 (MQ=255)
           cGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCCAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGt                                                                                                                                          <  1:1571854/150‑1 (MQ=255)
            gACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGtt                                                                                                                                         <  2:1230044/150‑1 (MQ=255)
            gACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGtt                                                                                                                                         <  2:41475/150‑1 (MQ=255)
            gACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGtt                                                                                                                                         >  1:1630196/1‑150 (MQ=255)
                        gAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCGTGGCAGGacctgat                                                                                                                                          >  1:956854/1‑131 (MQ=255)
                        gAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCGTGGCAGGacctgat                                                                                                                                          <  2:956854/137‑7 (MQ=255)
                          acacTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCcca                                                                                                                           <  2:1630196/150‑1 (MQ=255)
                              ttCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGTTCCTGCCGTGGCAGGacc                                                                                                                                              >  2:1357018/1‑125 (MQ=255)
                              ttCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGTTCCTGCCGTGGCAGGacc                                                                                                                                              <  1:1357018/127‑3 (MQ=255)
                              ttCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCcaccac                                                                                                                       >  2:1566621/1‑150 (MQ=255)
                              ttCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCcaccac                                                                                                                       >  1:1044921/1‑150 (MQ=255)
                                ccGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCg                                                                                                                     <  2:249114/150‑1 (MQ=255)
                                            ccTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGTGTCCACCACCGGGTTACGGCGCa                                                                                                         <  2:1044921/150‑1 (MQ=255)
                                            ccTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGACCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCa                                                                                                         <  1:151754/150‑1 (MQ=255)
                                             cTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGGAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAg                                                                                                        >  2:1397356/1‑150 (MQ=255)
                                             cTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAg                                                                                                        >  1:1397445/1‑150 (MQ=255)
                                                  ccAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAAATCGTTCATGAACTAACCCACCGCCGGGGTACGGCCCAGGccac                                                                                                   >  1:1580441/1‑150 (MQ=255)
                                                   cAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACg                                                                                                                                                      >  2:42672/1‑98 (MQ=255)
                                                   cAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACg                                                                                                                                                      <  1:42672/98‑1 (MQ=255)
                                                   cAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGccacc                                                                                                  >  1:1911999/1‑150 (MQ=255)
                                                   cAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGccacc                                                                                                  <  2:1397445/150‑1 (MQ=255)
                                                   cAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGccacc                                                                                                  >  1:1508737/1‑150 (MQ=255)
                                                    aGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAg                                                                                                                               >  1:99643/1‑120 (MQ=255)
                                                    aGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAg                                                                                                                               <  2:99643/120‑1 (MQ=255)
                                                    aGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGccacca                                                                                                 <  2:883348/150‑1 (MQ=255)
                                                    aGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCATACCATCAGGGGCGCTGTTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGccacca                                                                                                 >  2:1020922/1‑150 (MQ=255)
                                                     ggATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGccacca                                                                                                 >  2:160646/1‑149 (MQ=255)
                                                     ggATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGccacca                                                                                                 <  1:160646/149‑1 (MQ=255)
                                                          ttGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTAcc                                                                                           <  2:1656912/150‑1 (MQ=255)
                                                           tGATCCGCGCCTTGTACTTTTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCa                                                                                          <  1:595710/150‑1 (MQ=255)
                                                            gATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCa                                                                                                                                                   <  2:1395307/92‑1 (MQ=255)
                                                            gATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCa                                                                                                                                                   >  1:1395307/1‑92 (MQ=255)
                                                                cgcgCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGcaca                                                                                     <  2:1593939/150‑1 (MQ=255)
                                                                            ttgttgTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCCGGTCGCGACCTACCAGCACACGGAAGCCCAGt                                                                         <  1:1020922/150‑1 (MQ=255)
                                                                                  tCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCg                                                                   >  2:1010656/1‑150 (MQ=255)
                                                                                          aCCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGtt                                                           <  2:1911999/150‑1 (MQ=255)
                                                                                            cGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGcc                                                                             <  2:385667/130‑1 (MQ=255)
                                                                                            cGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGcc                                                                             >  1:385667/1‑130 (MQ=255)
                                                                                                aCGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCCGCGTTGCGCAc                                                     >  2:816885/1‑150 (MQ=255)
                                                                                                   gTTGTACACGCGCAGGATGGGGCCTGGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCgg                                                  <  2:1415146/150‑1 (MQ=255)
                                                                                                   gTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCgg                                                  >  1:1725798/1‑150 (MQ=255)
                                                                                                         cacGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAg                                            >  2:235564/1‑150 (MQ=255)
                                                                                                            gcgcAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGcc                                         >  2:1369955/1‑150 (MQ=255)
                                                                                                            gcgcAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCAACTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGcc                                         >  2:337667/1‑150 (MQ=255)
                                                                                                                aGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTc                                                                       >  2:849489/1‑116 (MQ=255)
                                                                                                                aGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTc                                                                       <  1:849489/116‑1 (MQ=255)
                                                                                                                  gATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTc                                   >  2:926226/1‑150 (MQ=255)
                                                                                                                  gATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCACCGCGGCAATCTGCATGGGCGAACAATCCCGGGTCGCCGGCAGCGTTGCGCACCGGCTCAAGGGCAAtttt                                   >  1:716907/1‑147 (MQ=255)
                                                                                                                                 gCTGTTCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATggc                    <  1:1942459/150‑1 (MQ=255)
                                                                                                                                   tGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTgc                                                         <  2:620395/111‑1 (MQ=255)
                                                                                                                                   tGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTgc                                                         >  1:620395/1‑111 (MQ=255)
                                                                                                                                     aTCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCcaccac                                                                                                <  2:1550914/70‑1 (MQ=255)
                                                                                                                                     aTCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCcaccac                                                                                                >  1:1550914/1‑70 (MQ=255)
                                                                                                                                     aTCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGgcggcg                >  1:207074/1‑150 (MQ=255)
                                                                                                                                      tCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCg                                                   >  1:695038/1‑114 (MQ=255)
                                                                                                                                      tCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCg                                                   <  2:695038/114‑1 (MQ=255)
                                                                                                                                      tCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAg                                            <  1:569506/121‑1 (MQ=255)
                                                                                                                                      tCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAg                                            >  2:569506/1‑121 (MQ=255)
                                                                                                                                          gTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTgcg                                                        <  2:1186273/105‑1 (MQ=255)
                                                                                                                                          gTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTgcg                                                        >  1:1186273/1‑105 (MQ=255)
                                                                                                                                          gTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCaa                                             <  2:1274273/116‑1 (MQ=255)
                                                                                                                                          gTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCaa                                             >  1:1274273/1‑116 (MQ=255)
                                                                                                                                             cTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAg                                                                                                        <  2:1130006/54‑1 (MQ=255)
                                                                                                                                             cTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAg                                                                                                        >  1:1130006/1‑54 (MQ=255)
                                                                                                                                             cTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCgg                                                                  <  2:756300/92‑1 (MQ=255)
                                                                                                                                             cTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCgg                                                                  >  1:756300/1‑92 (MQ=255)
                                                                                                                                             cTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTc        <  1:1566621/150‑1 (MQ=255)
                                                                                                                                             atGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTc        >  2:1326448/2‑150 (MQ=255)
                                                                                                                                                ccACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTc        <  1:1413789/147‑1 (MQ=255)
                                                                                                                                                ccACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTc        >  2:1413789/1‑147 (MQ=255)
                                                                                                                                                ccACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTCCTg     <  1:1010656/150‑1 (MQ=255)
                                                                                                                                                 cACGGCAGGAACTCGTTGATGAACGAGCCCAACCCCGCGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTCCTGc    <  1:1397356/150‑1 (MQ=255)
                                                                                                                                                  aCGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGAt                       <  1:148788/130‑1 (MQ=255)
                                                                                                                                                  aCGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGAt                       >  2:148788/1‑130 (MQ=255)
                                                                                                                                                   cGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTCCTGCGa  <  2:1725798/150‑1 (MQ=255)
                                                                                                                                                   cGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTCCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTCCGATGGCGGCGCGGTCTTCCTGCGa  <  1:1521385/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CATTTCGGCCTCGACCTTTTCGGCGAACACTTCCGGGGTAAGGGCCTTGACCAGGATCTTGATCCGCGCCTTGTACTTGTTGTCACGGCGACCGTAACGGTTGTACACGCGCAGGATGGCGTCGAGGTAGCTGATCAGGTCCTGCCACGGCAGGAACTCGTTGATGAACGAGCCCACCACCGGGGTACGGCCCAGGCCACCACCTACCAGCACACGGAAGCCCAGTTCGCCGGCAGCGTTGCGCACCGGCTCAAGGCCAATGTCGTGCACTTCGATGGCGGCGCGGTCTTCCTGCGA  >  NC_002947/2708439‑2708735

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