Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,274,244 T→A 5.6% F19I (TTC→ATC) ‡ metZ → O‑succinylhomoserine sulfhydrylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,274,2440TA 5.6% 910.0 / 27.0 265F19I (TTC→ATC) ‡metZO‑succinylhomoserine sulfhydrylase
Reads supporting (aligned to +/- strand):  ref base T (126/124);  new base A (4/11);  total (130/135)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 278 total aligned reads displayed.

CATCGAGCAGGCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAACGTCTACTCG  >  NC_002947/2274095‑2274393
                                                                                                                                                     |                                                                                                                                                     
 aTCGAGCAGGCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGttt                                                                                                                                                      >  1:702970/1‑150 (MQ=255)
  tCGAGCAGGCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTc                                                                                                                                                     <  2:4692294/150‑1 (MQ=255)
  acGAGCAGGCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTc                                                                                                                                                     <  2:3356019/149‑1 (MQ=255)
      gCAGGCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGaca                                                                                                                                                 <  2:403871/150‑1 (MQ=255)
         ggCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGAAACCGACAccc                                                                                                                                              <  1:4379097/150‑1 (MQ=255)
          gCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGAAACCGACACCCt                                                                                                                                             <  1:5449940/150‑1 (MQ=255)
          gCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGAAACCGACACCCt                                                                                                                                             <  1:283579/150‑1 (MQ=255)
            ccGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTgg                                                                                                                                           >  1:5550118/1‑150 (MQ=255)
            ccGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTgg                                                                                                                                           >  2:2934148/1‑150 (MQ=255)
            ccGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTgg                                                                                                                                           >  1:3059976/1‑150 (MQ=255)
             cGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGGCCTCGAGGGTGTCGGTTTCGACACCCTGGc                                                                                                                                          <  2:2205165/150‑1 (MQ=255)
             cGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGAAACCGACACCCtcga                                                                                                                                          <  1:3328515/150‑4 (MQ=255)
                    aCTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTAcgc                                                                                                                                   >  2:1332204/1‑150 (MQ=255)
                     cTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTAcgcg                                                                                                                                  >  1:191617/1‑150 (MQ=255)
                      ttGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACgcgc                                                                                                                                 <  2:2399516/150‑1 (MQ=255)
                      ttGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACgcgc                                                                                                                                 >  2:5316771/1‑150 (MQ=255)
                       tGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGTAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCt                                                                                                                                <  2:5527776/150‑1 (MQ=255)
                        ggCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTg                                                                                                                               <  1:2586002/150‑1 (MQ=255)
                        ggCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTg                                                                                                                               >  1:883460/1‑150 (MQ=255)
                        ggCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTg                                                                                                                               <  2:5135339/150‑1 (MQ=255)
                        ggCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTg                                                                                                                               <  1:1259293/150‑1 (MQ=255)
                           cAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGcc                                                                                                                            >  1:4185486/1‑150 (MQ=255)
                           cAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGcc                                                                                                                            <  2:3558535/150‑1 (MQ=255)
                                tCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAAc                                                                                                                       >  1:3962997/1‑150 (MQ=255)
                                 cAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAAcc                                                                                                                      <  1:1337944/150‑1 (MQ=255)
                                  aaGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCg                                                                                                                     >  2:2597904/1‑150 (MQ=255)
                                  aaGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCg                                                                                                                     >  1:5281755/1‑150 (MQ=255)
                                  aaGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCg                                                                                                                     >  1:4864561/1‑150 (MQ=255)
                                  aaGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGGTTCGACACCCTGGCGGTACGCGCTGGCCACAACCg                                                                                                                     >  2:3507290/1‑150 (MQ=255)
                                    gAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTa                                                                                                                   >  1:1274082/1‑150 (MQ=255)
                                    gAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTa                                                                                                                   <  2:1778930/150‑1 (MQ=255)
                                     aaCCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTAc                                                                                                                  <  2:275803/150‑1 (MQ=255)
                                       ccAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTAccc                                                                                                                >  2:5001573/1‑150 (MQ=255)
                                         aGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCg                                                                                                              >  2:631002/1‑150 (MQ=255)
                                              gTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAgg                                                                                                         >  1:300488/1‑150 (MQ=255)
                                                 aGCGCGATCATCTACCTTTACAACAACTGATTTTGGGAGCGTCTGCATGACGGATCATTGGGTTGCCGGGCGACTGGTCAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCACGGAAGGCGa                                                                                                      <  2:2563144/150‑1 (MQ=255)
                                                 aGCGCGATCATCGACCTTTACACCAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGa                                                                                                      <  1:5019631/150‑1 (MQ=255)
                                                 aGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGa                                                                                                      <  2:2753871/150‑1 (MQ=255)
                                                    gcgATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGca                                                                                                   <  1:1208895/150‑1 (MQ=255)
                                                      gATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGcaca                                                                                                 >  1:4217831/1‑150 (MQ=255)
                                                      gATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGcaca                                                                                                 >  2:3100605/1‑150 (MQ=255)
                                                            cGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAg                                                                                           >  2:36144/1‑150 (MQ=255)
                                                                cTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCt                                                                                       >  2:4733506/1‑150 (MQ=255)
                                                                 tttACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTg                                                                                      <  1:171742/150‑1 (MQ=255)
                                                                  ttACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGt                                                                                     <  2:2064623/150‑1 (MQ=255)
                                                                    acaacaACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTc                                                                                   >  2:4623552/1‑150 (MQ=255)
                                                                      aacaacTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCt                                                                                 >  2:1667241/1‑150 (MQ=255)
                                                                      aacaacTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCt                                                                                 >  2:4145707/1‑150 (MQ=255)
                                                                       acaacTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTg                                                                                <  1:2714986/150‑1 (MQ=255)
                                                                         aacTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAGGGCGAGCACAGCGAAGCGCTGTTCCTGAc                                                                              >  1:1606671/1‑150 (MQ=255)
                                                                         aacTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGAc                                                                              >  1:4384611/1‑150 (MQ=255)
                                                                         aacTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGAc                                                                              >  2:1113736/1‑150 (MQ=255)
                                                                         aacTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGAc                                                                              >  2:383197/1‑150 (MQ=255)
                                                                            tGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGt                                                                                     >  1:1073484/1‑140 (MQ=255)
                                                                            tGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTc                                                                           >  2:3595464/1‑150 (MQ=255)
                                                                            tGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTc                                                                           >  1:2610779/1‑150 (MQ=255)
                                                                            tGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTc                                                                           >  2:4272871/1‑150 (MQ=255)
                                                                             gATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTcc                                                                          <  2:4247462/150‑1 (MQ=255)
                                                                             gATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTcc                                                                          <  2:4246924/150‑1 (MQ=255)
                                                                              tttttGTGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGTCTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCa                                                                         <  1:3807855/149‑1 (MQ=255)
                                                                                  tGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTa                                                                     >  2:4343229/1‑150 (MQ=255)
                                                                                  tGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTa                                                                     <  2:278597/150‑1 (MQ=255)
                                                                                   gggAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTAt                                                                    <  2:1108813/150‑1 (MQ=255)
                                                                                     gAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTAtgt                                                                  >  2:3664947/1‑150 (MQ=255)
                                                                                        cGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGtt                                                               >  2:391378/1‑150 (MQ=255)
                                                                                         gACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTc                                                              <  1:4182985/150‑1 (MQ=255)
                                                                                          aCGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTcc                                                             >  1:3202640/1‑150 (MQ=255)
                                                                                              cATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCAc                                                         <  2:3856330/150‑1 (MQ=255)
                                                                                                tGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACgg                                                       >  1:3579467/1‑150 (MQ=255)
                                                                                                tGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACgg                                                       <  1:1971189/150‑1 (MQ=255)
                                                                                                 gACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGc                                                      >  2:1053799/1‑150 (MQ=255)
                                                                                                    ggATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGc                                                   >  1:3766665/1‑150 (MQ=255)
                                                                                                     gATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGcc                                                  <  1:3421233/150‑1 (MQ=255)
                                                                                                        cAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGAt                                               >  1:5118690/1‑150 (MQ=255)
                                                                                                          aTGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGCTCCGCACGGCAGCCGAtgc                                             <  1:4618118/150‑1 (MQ=255)
                                                                                                                tGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGAtgctgctgc                                       <  1:5601433/150‑1 (MQ=255)
                                                                                                                    gggCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTgcgcgc                                   >  1:4966344/1‑150 (MQ=255)
                                                                                                                           tGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCg                            >  2:1244019/1‑150 (MQ=255)
                                                                                                                             gACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGc                          >  1:139747/1‑150 (MQ=255)
                                                                                                                               cAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGa                        >  1:3670814/1‑150 (MQ=255)
                                                                                                                               cAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGa                        >  1:840957/1‑150 (MQ=255)
                                                                                                                               cAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGa                        <  2:1607195/150‑1 (MQ=255)
                                                                                                                                     ccTCGAGGGGGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAACACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGcc                  <  1:5312176/150‑1 (MQ=255)
                                                                                                                                      cTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCg                 <  2:489958/150‑1 (MQ=255)
                                                                                                                                      cTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCg                 >  1:3498878/1‑150 (MQ=255)
                                                                                                                                       tCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCgg                <  2:5434215/150‑1 (MQ=255)
                                                                                                                                       tCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCgg                <  1:1113736/150‑1 (MQ=255)
                                                                                                                                       tCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCgg                <  2:1274082/150‑1 (MQ=255)
                                                                                                                                       tCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCgg                >  1:3840791/1‑150 (MQ=255)
                                                                                                                                       gccaGGGTGTCGAAACCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCgg                >  2:3763324/4‑150 (MQ=255)
                                                                                                                                        cGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCggg               <  2:903318/150‑1 (MQ=255)
                                                                                                                                         naGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGc              >  1:2389497/2‑150 (MQ=255)
                                                                                                                                         gAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGc              >  1:1756347/1‑150 (MQ=255)
                                                                                                                                         gAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCGGACCTCCAGCTATGTGTGCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGc              >  2:3282018/1‑150 (MQ=255)
                                                                                                                                           gggTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGTTTCGCCggcaaa            <  1:2889600/150‑5 (MQ=255)
                                                                                                                                              tgtCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAACGt         <  2:5550118/150‑1 (MQ=255)
                                                                                                                                                tCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAACGTCt       >  1:722894/1‑150 (MQ=255)
                                                                                                                                                    tttCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGc                              <  2:847341/123‑1 (MQ=255)
                                                                                                                                                    tttCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAACGTCTACTc   >  2:522292/1‑150 (MQ=255)
                                                                                                                                                     ttCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCACCGTCTACTcg  >  2:5158798/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATCGAGCAGGCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAACGTCTACTCG  >  NC_002947/2274095‑2274393

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