Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,269,250 G→A 5.2% G18D (GGT→GAT)  metZ → O‑succinylhomoserine sulfhydrylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,269,2500GA 5.2% 358.1 / 5.3 116G18D (GGT→GAT) metZO‑succinylhomoserine sulfhydrylase
Reads supporting (aligned to +/- strand):  ref base G (56/53);  new base A (2/4);  total (58/58)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

Only 100 of 122 total aligned reads displayed.

ATCGAGCAGGCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAACGTCTAC  >  NC_002947/2269104‑2269398
                                                                                                                                                   |                                                                                                                                                    
  cGAGCAGGCCCGTAACGACTTGGCCAAGGTCAGAGAACCAGCCCGTCAGCGCGATCATCGACCTTTACAACAACTCATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGCCCACTGGACAGTGAACACGAGGGTGTCGtttt                                                                                                                                                   >  1:206828/1‑149 (MQ=255)
   gAGCAGGCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTc                                                                                                                                                  <  2:106248/148‑1 (MQ=255)
         gCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGAAACCGACAccc                                                                                                                                           <  4:214537/149‑1 (MQ=255)
           ccGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGCCCCCCTg                                                                                                                                         >  3:286560/1‑149 (MQ=255)
            cGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTgg                                                                                                                                        <  1:75406/149‑1 (MQ=255)
              tAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGAAACCGACACCCtcgag                                                                                                                                      <  4:289967/149‑5 (MQ=255)
              tAACGACCTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCg                                                                                                                                      <  8:72448/149‑1 (MQ=255)
                   aCTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTAcg                                                                                                                                 >  2:217298/1‑149 (MQ=255)
                   aCTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGAGGCCGGGCGACTGGACAGTGACCTCGAGGGTGGCGGTTTCGACACCCTGGCGGTAcg                                                                                                                                 >  8:213080/1‑149 (MQ=255)
                    cTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTAcgc                                                                                                                                >  4:75153/1‑149 (MQ=255)
                       ggCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCt                                                                                                                             >  4:49579/1‑149 (MQ=255)
                        gCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTg                                                                                                                            <  7:128087/149‑1 (MQ=255)
                        gCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTg                                                                                                                            <  5:59636/149‑1 (MQ=255)
                        gCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTg                                                                                                                            >  2:50330/1‑149 (MQ=255)
                         ccAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGAGTTTGGGAGCGACGGCAGGACGGAGCACTGGGATGCCGTGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGCCACCCGGGCGGTACGAGCTgg                                                                                                                           >  5:168282/1‑149 (MQ=255)
                              gTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTc                                                                                                                                                  <  1:147777/121‑1 (MQ=255)
                              gTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTc                                                                                                                                                  >  2:147777/1‑121 (MQ=255)
                              gTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCaa                                                                                                                       >  1:170997/1‑148 (MQ=255)
                                 aaGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGCCGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACCCCCTGGCGGTACGCGCTGGCCAAAAcc                                                                                                                   >  3:231743/1‑149 (MQ=255)
                                 aaGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAAcc                                                                                                                   >  3:108650/1‑149 (MQ=255)
                                 aaGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAAcc                                                                                                                   <  3:147824/149‑1 (MQ=255)
                                   aGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCg                                                                                                                  <  6:217624/149‑1 (MQ=255)
                                    gAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCg                                                                                                                  >  4:291759/1‑148 (MQ=255)
                                      aCCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTAc                                                                                                               <  5:85735/149‑1 (MQ=255)
                                       ccAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTAcc                                                                                                              <  4:211042/149‑1 (MQ=255)
                                         aGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTAcc                                                                                                              >  7:58166/1‑147 (MQ=255)
                                          ggCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCg                                                                                                           >  2:209657/1‑149 (MQ=255)
                                          ggCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGCCGGTACGCGCTGGCCAAAACCGTACCCCg                                                                                                           >  4:168328/1‑149 (MQ=255)
                                                 cgcggGATCAGCGCACTTTCCCACACCTGCTTTTGGGCGCGACGGCCTTACGGCTCTCTTGGATGCCGGGCTACTGGACATTGACCTCGAGGGTGGCGTTTTCGACACACTGGCGGTACGCGCTGGCCAACACCGTCCCCCGGACGGc                                                                                                     <  6:328500/147‑1 (MQ=255)
                                                  gcgcGATCATCGACCTTTAGAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGa                                                                                                   <  8:122053/149‑1 (MQ=255)
                                                   cgcgATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAg                                                                                                  <  1:217298/149‑1 (MQ=255)
                                                      gATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGcac                                                                                               <  7:302237/149‑1 (MQ=255)
                                                          atcGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGc                                                                                            >  5:316816/1‑148 (MQ=255)
                                                           tcGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCg                                                                                           >  3:11153/1‑148 (MQ=255)
                                                           tcGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCg                                                                                           >  4:212202/1‑148 (MQ=255)
                                                               ccTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTgg                                                                                                                           >  5:172471/1‑112 (MQ=255)
                                                               ccTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTgg                                                                                                                           <  6:172471/112‑1 (MQ=255)
                                                               ccTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAgcg                                                                                      <  2:245350/149‑1 (MQ=255)
                                                               ccTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAgcg                                                                                      <  8:79784/149‑1 (MQ=255)
                                                                 tttACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCt                                                                                    >  5:83704/1‑149 (MQ=255)
                                                                 tttACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCt                                                                                    >  8:156611/1‑149 (MQ=255)
                                                                  ttACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACAGTACCCCCGAAGGCGAGCACAGCGAAGCGCTg                                                                                   >  6:232984/1‑149 (MQ=255)
                                                                    tcatgatcTGATTTTGGTTGCGACGGCATGAAGGATCAATGGGATGCCGGGCGACTGGACAGTGAGCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGt                                                                                  <  8:72764/141‑1 (MQ=255)
                                                                    acaacaACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGt                                                                                  <  4:300935/148‑1 (MQ=255)
                                                                     caacaaCTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTc                                                                                <  6:83704/149‑1 (MQ=255)
                                                                      aacaacTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTTGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTcc                                                                               >  5:247848/1‑149 (MQ=255)
                                                                      aacaacTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGACGGCGAGCACAGCGAAGCGCTGTTcc                                                                               >  5:23200/1‑149 (MQ=255)
                                                                      aacaacTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTcc                                                                               >  5:263533/1‑149 (MQ=255)
                                                                      aacaacTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTcc                                                                               >  3:316791/1‑149 (MQ=255)
                                                                      aacaacTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTcc                                                                               >  8:288915/1‑149 (MQ=255)
                                                                      aacaacTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGAGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACAGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTcc                                                                               >  4:18405/1‑149 (MQ=255)
                                                                         aacTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGa                                                                            >  3:28310/1‑149 (MQ=255)
                                                                              aTTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTcc                                                                       <  1:50330/149‑1 (MQ=255)
                                                                              aTTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTc                                                                        <  4:28310/148‑1 (MQ=255)
                                                                              aTTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTc                                                                        >  5:291297/1‑148 (MQ=255)
                                                                               ttttGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTcc                                                                       <  3:291759/148‑1 (MQ=255)
                                                                                   gggAGCGAGGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCt                                                                   <  3:244058/148‑1 (MQ=255)
                                                                                    ggAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTTGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTAt                                                                 <  6:247848/149‑1 (MQ=255)
                                                                                    ggAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTAt                                                                 <  3:168328/149‑1 (MQ=255)
                                                                                    ggAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTAt                                                                 <  4:108650/149‑1 (MQ=255)
                                                                                    ggAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTAt                                                                 <  3:212202/149‑1 (MQ=255)
                                                                                      aGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGACAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTAtgt                                                               <  5:232984/149‑1 (MQ=255)
                                                                                          aCGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTc                                                           <  7:288915/149‑1 (MQ=255)
                                                                                             gCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCg                                                         >  7:88389/1‑148 (MQ=255)
                                                                                                tGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACg                                                     <  3:75153/149‑1 (MQ=255)
                                                                                                     gATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGc                                                <  8:88389/149‑1 (MQ=255)
                                                                                                       tCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGcc                                               >  2:54443/1‑148 (MQ=255)
                                                                                                          aTGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGAtg                                           <  7:93345/149‑1 (MQ=255)
                                                                                                            gggATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTg                                                                             >  3:234448/1‑113 (MQ=255)
                                                                                                              gATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACgg                                                    <  5:271126/136‑1 (MQ=255)
                                                                                                              gATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACgg                                                    >  6:271126/1‑136 (MQ=255)
                                                                                                              gATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGAtgctgc                                       >  2:246354/1‑149 (MQ=255)
                                                                                                                   cGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGc                                                <  4:97470/135‑1 (MQ=255)
                                                                                                                   cGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGc                                                >  3:97470/1‑135 (MQ=255)
                                                                                                                      gCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCt                               >  6:54549/1‑149 (MQ=255)
                                                                                                                      gCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCGGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCt                               <  5:252763/149‑1 (MQ=255)
                                                                                                                          cTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTg                                                                             <  1:327734/99‑1 (MQ=255)
                                                                                                                          cTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTg                                                                             >  2:327734/1‑99 (MQ=255)
                                                                                                                           tGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGcc                          >  8:273192/1‑149 (MQ=255)
                                                                                                                             gACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCgg                        >  5:215967/1‑149 (MQ=255)
                                                                                                                             gACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCg                         <  3:18405/148‑1 (MQ=255)
                                                                                                                               cAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCg                      <  4:111077/149‑1 (MQ=255)
                                                                                                                                      cTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGcc               >  1:133218/1‑149 (MQ=255)
                                                                                                                                        ccaGGGTGTCGAAACCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTcc                                                                               <  6:78721/81‑1 (MQ=38)
                                                                                                                                        ccaGGGTGTCGAAACCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTcc                                                                               >  5:78721/3‑83 (MQ=38)
                                                                                                                                         gAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCggg            >  3:186852/1‑149 (MQ=255)
                                                                                                                                         gAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCgg             <  6:263533/148‑1 (MQ=255)
                                                                                                                                         caGGGTGTCGAAACCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCggg            >  8:141901/2‑149 (MQ=255)
                                                                                                                                           gggTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCgg             >  2:195840/1‑146 (MQ=255)
                                                                                                                                           gggTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCgg             <  1:195840/146‑1 (MQ=255)
                                                                                                                                           gggTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGc           >  3:34715/1‑148 (MQ=255)
                                                                                                                                           gggTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGc           >  1:268055/1‑148 (MQ=255)
                                                                                                                                              tgtCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAACg       <  2:170997/149‑1 (MQ=255)
                                                                                                                                                tCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAACGTc     <  6:291297/149‑1 (MQ=255)
                                                                                                                                                 cGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAg                                                                                        <  7:155206/65‑1 (MQ=255)
                                                                                                                                                  ggTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAACGTCTa   >  7:222379/1‑149 (MQ=255)
                                                                                                                                                  ggTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAACGTCTa   <  5:54549/149‑1 (MQ=255)
                                                                                                                                                   gTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAACGTCTAc  >  8:198998/1‑149 (MQ=255)
                                                                                                                                                   gTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAACGTCTAc  <  2:268055/149‑1 (MQ=255)
                                                                                                                                                   gTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGGTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAACGTCTAc  <  7:213080/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATCGAGCAGGCCCGTAACGACTTGGCCAAGGTCAAGAACCAGGCCGTCAGCGCGATCATCGACCTTTACAACAACTGATTTTGGGAGCGACGGCATGACGGATCAATGGGATGCCGGGCGACTGGACAGTGACCTCGAGGGTGTCGGTTTCGACACCCTGGCGGTACGCGCTGGCCAAAACCGTACCCCGGAAGGCGAGCACAGCGAAGCGCTGTTCCTGACCTCCAGCTATGTGTTCCGCACGGCAGCCGATGCTGCTGCGCGCTTTGCCGGCGAAACGCCGGGCAACGTCTAC  >  NC_002947/2269104‑2269398

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