Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,829,573 T→G 5.3% Q20P (CAG→CCG)  PP_3382 ← gluconate 2‑dehydrogenase cytochrome c subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,829,5730TG 5.3% 286.4 / 6.6 95Q20P (CAG→CCG) PP_3382gluconate 2‑dehydrogenase cytochrome c subunit
Reads supporting (aligned to +/- strand):  ref base T (44/46);  new base G (2/3);  total (46/49)
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.59e-01

GGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCAGCCAC  >  NC_002947/3829425‑3829718
                                                                                                                                                    |                                                                                                                                                 
ggtggAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCtg                                                                                                                                                  <  2:1657556/150‑1 (MQ=255)
 gtggAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGAAGGCAAGACGTCATCCGCGATTCGCACTTGTGCCGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCtgt                                                                                                                                                 <  1:1198231/150‑1 (MQ=255)
    gAGTACACCGTGCCAATCGGTGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCa                                                                                                                                              <  2:1639309/150‑1 (MQ=255)
    gAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCa                                                                                                                                              >  2:1702226/1‑150 (MQ=255)
       tACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGcc                                                                                                                                           >  1:1047175/1‑150 (MQ=255)
       tACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGcc                                                                                                                                           <  1:1256624/150‑1 (MQ=255)
         cacCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGc                                                                                                                                         >  1:1255057/1‑150 (MQ=255)
          acCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCa                                                                                                                                        >  2:1311383/1‑150 (MQ=255)
          acCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCGGTGcggctgca                                                                                                                                        >  1:271586/1‑143 (MQ=255)
          acCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCGGTGCGGCTGCa                                                                                                                                        >  1:272508/1‑150 (MQ=255)
              tGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCGGTGCGGCTGCACAgg                                                                                                                                    <  2:272508/150‑1 (MQ=255)
              tGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCGGTGCGGCTGCACAgg                                                                                                                                    <  2:271586/150‑1 (MQ=255)
               gCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCGGTGCGGCTGCACAGgc                                                                                                                                   <  2:584943/150‑1 (MQ=255)
                  aaTCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGcc                                                                                                                                >  1:850550/1‑150 (MQ=255)
                   aTCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGccc                                                                                                                               <  2:355862/150‑1 (MQ=255)
                      ggCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGc                                                                                                                            <  1:1588806/150‑1 (MQ=255)
                       gCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCa                                                                                                                           <  1:253960/150‑1 (MQ=255)
                       gCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCa                                                                                                                           <  2:92771/150‑1 (MQ=255)
                         gTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCAcc                                                                                                                         >  1:1530958/1‑150 (MQ=255)
                           ttCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAg                                                                                                                       <  2:588055/150‑1 (MQ=255)
                               aTCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTg                                                                                                                   <  2:1918768/150‑1 (MQ=255)
                                tCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTgg                                                                                                                  <  1:1830408/150‑1 (MQ=255)
                                   gCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGa                                                                                                               >  1:1345653/1‑150 (MQ=255)
                                     aaCCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCt                                                                                                                                                   >  2:1590323/1‑112 (MQ=255)
                                     aaCCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCt                                                                                                                                                   <  1:1590323/112‑1 (MQ=255)
                                        ccgccgGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAAc                                                                                                          >  1:1299033/1‑150 (MQ=255)
                                        ccgccgGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAAc                                                                                                          <  1:1297756/150‑1 (MQ=255)
                                                aaaGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCtt                                                                                                  >  1:1063237/1‑150 (MQ=255)
                                                 aaGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTc                                                                                                 <  1:1339980/150‑1 (MQ=255)
                                                   gggTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCAt                                                                                               >  2:752441/1‑150 (MQ=255)
                                                     gTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGGCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGc                                                                                             <  2:1255057/150‑1 (MQ=255)
                                                     gTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGc                                                                                             <  2:575778/150‑1 (MQ=255)
                                                     gTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGc                                                                                             <  2:1530958/150‑1 (MQ=255)
                                                     gTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGc                                                                                             <  2:1063237/150‑1 (MQ=255)
                                                          ccgccTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCa                                                                                                                                              >  2:1724186/1‑96 (MQ=255)
                                                          ccgccTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCa                                                                                                                                              <  1:1724186/96‑1 (MQ=255)
                                                          ccgccTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGCTGGCGATCAACAGTGTCTTCAtccccctc                                                                                        >  2:705466/1‑143 (MQ=255)
                                                           cgccTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCg                                                                                       >  2:286089/1‑150 (MQ=255)
                                                            gccTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGc                                                                                      >  2:1090608/1‑150 (MQ=255)
                                                              cTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCCGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAg                                                                                                                       <  1:627364/115‑1 (MQ=255)
                                                              cTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAg                                                                                                                       >  2:627364/1‑115 (MQ=255)
                                                                     ggTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGc                                                                             <  1:705466/150‑1 (MQ=255)
                                                                         tGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCATGCCTgg                                                                         <  1:752441/150‑1 (MQ=255)
                                                                          ggCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGa                                                                        <  1:286089/150‑1 (MQ=255)
                                                                           gCAGGCCACGCAGTTACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCt                                                                                            <  1:1824586/129‑1 (MQ=255)
                                                                           gCAGGCCACGCAGTTACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCt                                                                                            >  2:1824586/1‑129 (MQ=255)
                                                                           gCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGAt                                                                                                              <  1:1196869/111‑1 (MQ=255)
                                                                           gCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGCTCTGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGAc                                                                       >  1:621604/1‑150 (MQ=255)
                                                                           gCAGGCCACGCAGTCACCGGCGCGGGCCAGGCACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGAt                                                                                                              >  2:1196869/1‑111 (MQ=255)
                                                                            cAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGAcc                                                                      >  1:872750/1‑150 (MQ=255)
                                                                             aGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTc                                                                                                    <  1:729216/119‑1 (MQ=255)
                                                                             aGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTc                                                                                                    >  2:729216/1‑119 (MQ=255)
                                                                             aGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCa                                                                     >  1:491709/1‑150 (MQ=255)
                                                                             aGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACATGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCa                                                                     >  1:1320036/1‑150 (MQ=255)
                                                                                ccACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAgtg                                                                                                      <  2:1450740/114‑1 (MQ=255)
                                                                                ccACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAgtg                                                                                                      >  1:1450740/1‑114 (MQ=255)
                                                                                  aCGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGAt                                                                                                              <  1:6272/104‑1 (MQ=255)
                                                                                  aCGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGAt                                                                                                              >  2:6272/1‑104 (MQ=255)
                                                                                   cGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGAc                                                               <  2:1320036/150‑1 (MQ=255)
                                                                                          aCCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCCGGGAGGCGAGCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGtt                                                        >  2:751061/1‑150 (MQ=255)
                                                                                          aCCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCAt                                                                                         <  1:850684/117‑1 (MQ=255)
                                                                                          aCCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCAt                                                                                         >  2:850684/1‑117 (MQ=255)
                                                                                                      cAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTg                                            >  1:1361254/1‑150 (MQ=255)
                                                                                                       aGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCTCCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGTTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGc                                           <  1:782227/150‑1 (MQ=255)
                                                                                                         gTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGc                                                                                             >  2:3439/1‑98 (MQ=255)
                                                                                                         gTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGc                                                                                             <  1:3439/98‑1 (MQ=255)
                                                                                                           aCTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGcc                                                                                                                                <  2:1500700/61‑1 (MQ=255)
                                                                                                           aCTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGcc                                                                                                                                >  1:1500700/1‑61 (MQ=255)
                                                                                                           aCTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCAc                                       <  2:1047175/150‑1 (MQ=255)
                                                                                                            cTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGAc                                                                       >  1:1320479/1‑117 (MQ=255)
                                                                                                            cTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGAc                                                                       <  2:1320479/117‑1 (MQ=255)
                                                                                                            cTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACg                                      >  2:503309/1‑150 (MQ=255)
                                                                                                            cTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCCAACAGCACCAGGACATCAATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACg                                      >  1:1691691/1‑150 (MQ=255)
                                                                                                             tCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACgg                                     >  2:904904/1‑150 (MQ=255)
                                                                                                              cGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGa                                    <  1:1702226/150‑1 (MQ=255)
                                                                                                                cccAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACgg                                     >  2:1596724/1‑147 (MQ=255)
                                                                                                                cccAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACgg                                     <  1:1596724/147‑1 (MQ=255)
                                                                                                                 ccAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCATCAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATgg                                 >  1:1143788/1‑150 (MQ=255)
                                                                                                                 ccAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATgg                                 <  1:503309/150‑1 (MQ=255)
                                                                                                                     gCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGctc                                                                                    >  2:911689/1‑95 (MQ=255)
                                                                                                                     gCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGctc                                                                                    <  1:911689/95‑1 (MQ=255)
                                                                                                                           cTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTTCTCAAGc                                                                             <  2:522937/96‑1 (MQ=255)
                                                                                                                           cTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTTCTCAAGc                                                                             >  1:522937/1‑96 (MQ=255)
                                                                                                                            tGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGAc                      <  1:950289/150‑1 (MQ=255)
                                                                                                                            tGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGAc                      <  1:1311383/150‑1 (MQ=255)
                                                                                                                            tGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGAc                      >  1:733582/1‑150 (MQ=255)
                                                                                                                             gcgcATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGttttt                                                     >  1:1502496/1‑118 (MQ=255)
                                                                                                                             gcgcATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGttttt                                                     <  2:1502496/118‑1 (MQ=255)
                                                                                                                               gcATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCgggg                                                          >  1:1245169/1‑111 (MQ=255)
                                                                                                                               gcATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCgggg                                                          <  2:1245169/111‑1 (MQ=255)
                                                                                                                               gcATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGCCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGAcc                     <  1:1535892/148‑1 (MQ=255)
                                                                                                                               gcATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGCCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGAcc                     >  2:1535892/1‑148 (MQ=255)
                                                                                                                                 aTCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTa                 >  2:81568/1‑150 (MQ=255)
                                                                                                                                    gTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGt              <  2:1186857/150‑1 (MQ=255)
                                                                                                                                           gCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCa       <  1:904904/150‑1 (MQ=255)
                                                                                                                                               cGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGTGCGGGCCGAc                      <  1:1579538/131‑1 (MQ=255)
                                                                                                                                               cGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGTGCGGGCCGAc                      >  2:1579538/1‑131 (MQ=255)
                                                                                                                                                ggCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCAGccac  <  1:538848/150‑1 (MQ=255)
                                                                                                                                                ggCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCAGccac  >  1:160764/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGTGGAGTACACCGTGCCAATCGGCGTTTCCATCGGCAACCCGCCGGCAAAGGGTTTGCCGCCTTTGGCGGTATGGCAGGCCACGCAGTCACCGGCGCGGGCCAGGTACTCGCCCAGGCGTACCTGCGCATCGTCATTCGCTACGGCCTGTGCAGCCGCACCGGCGCCCAGCACCAGGGTGGCGATCAACAGTGTCTTCATGCTCATCGCTCCTCAAGCCTGGACCAACGGACCGGGGTTTTTCAGGTACTGCTCACGGATGGCGCGGGCCGACCAGTACGTGAGCGCAGCCAC  >  NC_002947/3829425‑3829718

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