Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,002,533 T→C 8.9% Y225H (TAC→CAC)  allE → S‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,002,5330TC 8.9% 282.3 / 13.4 103Y225H (TAC→CAC) allES‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base T (53/40);  new base C (4/5);  total (58/45)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

Only 100 of 106 total aligned reads displayed.

GAAAACCACCCGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCT  >  NC_002947/4002388‑4002681
                                                                                                                                                 |                                                                                                                                                    
gAAAACCACCCGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGTGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCg                                                                                                                                                  >  1:1373715/1‑150 (MQ=255)
 aaaaCCACCCGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGc                                                                                                                                                 <  2:1645424/150‑1 (MQ=255)
    accaccCGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTACACCGCCtt                                                                                                                                              <  2:1738005/150‑2 (MQ=255)
     ccaccCGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGgtgt                                                                                                                                                      >  1:386818/1‑141 (MQ=255)
     ccaccCGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGgtgt                                                                                                                                                      <  2:386818/141‑1 (MQ=255)
     ccaccCGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTACACCGCCttg                                                                                                                                             <  2:1350682/150‑3 (MQ=255)
       accCGTTTTGTCGACATGGCCGACCTGCGCCACGACAGGCACGTGAACATCGTCACCTTCCAGCCGGGCGGCGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGGACGTGCGGGAAGGCAAGGCGGTGTACCGCCTGAAc                                                                                                                                           >  1:1507073/1‑150 (MQ=255)
        cccGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTACGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAAcc                                                                                                                                          >  1:1217566/1‑150 (MQ=255)
          cGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACTTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAg                                                                                                                                        >  2:1160855/1‑150 (MQ=255)
                gTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTg                                                                                                                                   <  1:647207/149‑1 (MQ=255)
                gTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTg                                                                                                                                   >  2:647207/1‑149 (MQ=255)
                  cGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGt                                                                                                                                <  2:1217566/150‑1 (MQ=255)
                   gACATGGCCGACATTCGCCACGACATGCACGTGAACATCGGCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGAAGTACCGCCTGAACCAGGACTGCGTc                                                                                                                               >  1:2300291/1‑150 (MQ=255)
                   gACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTAGGTGCTGGAAGGCAAGGCGGTGAACCGCCTGAACCAGGACTGGGTc                                                                                                                               >  1:775746/1‑150 (MQ=255)
                    aCATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCg                                                                                                                              >  2:382389/1‑150 (MQ=255)
                                    ccACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGaa                                                                                                                      >  1:943576/1‑142 (MQ=255)
                                    ccACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGa                                                                                                              >  2:1904312/1‑150 (MQ=255)
                                            tGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCAtgt                                                                                                      >  1:1045237/1‑150 (MQ=255)
                                                 gTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTg                                                                                                 >  1:1694782/1‑150 (MQ=255)
                                                 gTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTg                                                                                                 >  2:2227372/1‑150 (MQ=255)
                                                 gTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTg                                                                                                 >  2:1550715/1‑150 (MQ=255)
                                                      cATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTACACCGCCtt                                                                                                                                              >  2:1770930/1‑99 (MQ=255)
                                                      cATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTACACCGCCtt                                                                                                                                              <  1:1770930/100‑2 (MQ=255)
                                                       aTCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGcc                                                                                           >  2:2368457/1‑150 (MQ=255)
                                                           tCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTgg                                                                                                                                  >  2:1298458/1‑107 (MQ=255)
                                                           tCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTgg                                                                                                                                  <  1:1298458/107‑1 (MQ=255)
                                                           tCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCt                                                                                       >  1:1981701/1‑150 (MQ=255)
                                                           tCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCAGGGTACCCCCTCAACCAGTACTCGGTCGAAGTGGAAGCCGGCCACTTCCTGCGGCTGCCCGCAttgc                                                                                       >  2:747344/1‑148 (MQ=255)
                                                           tCACATTCCAGCCGGGCGGTGTGATTCCGTTCGTAAAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACGAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCt                                                                                       <  1:707678/150‑1 (MQ=255)
                                                               cTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGccc                                                                                   >  1:186121/1‑150 (MQ=255)
                                                                ttCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGaa                                                                                                                      >  2:698139/1‑114 (MQ=255)
                                                                ttCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGaa                                                                                                                      <  1:698139/114‑1 (MQ=255)
                                                                ttCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCg                                                                                  >  2:1088208/1‑150 (MQ=255)
                                                                     gCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGccc                                                                                   <  2:1612671/144‑1 (MQ=255)
                                                                     gCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGccc                                                                                   >  1:1612671/1‑144 (MQ=255)
                                                                     gCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCG‑‑GTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCt                                                                           >  1:1626609/1‑150 (MQ=255)
                                                                              tgtgATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGc                                                                    >  2:2122605/1‑150 (MQ=255)
                                                                               gtgATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTAGCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACggc                                                                   <  2:1686810/150‑3 (MQ=255)
                                                                                  aTTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCcggc                                                                >  2:2003872/1‑150 (MQ=255)
                                                                                   ttCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGc                                                                             <  2:1081788/136‑1 (MQ=255)
                                                                                   ttCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGc                                                                             >  1:1081788/1‑136 (MQ=255)
                                                                                   ttCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCcggcg                                                               <  2:1373715/150‑1 (MQ=255)
                                                                                    tCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGGCTAGGTCGAAGTGGAAGCCCCCGACGCTATCTGTCTTGGCTCCTTCTGCCCGCAAGCCTGCTACGCcggcgg                                                              <  1:1023474/150‑1 (MQ=255)
                                                                                       gTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGccc                                                           >  1:1656536/1‑150 (MQ=255)
                                                                                         tCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAg                                                         <  1:382389/150‑1 (MQ=255)
                                                                                           gCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCt                                                                                          <  2:715060/115‑1 (MQ=255)
                                                                                           gCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCt                                                                                          >  1:715060/1‑115 (MQ=255)
                                                                                             aGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCaa                                                     >  2:631383/1‑150 (MQ=255)
                                                                                             aGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCaa                                                     >  1:1338693/1‑150 (MQ=255)
                                                                                              gAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTa                                                                       >  1:2227998/1‑131 (MQ=255)
                                                                                              gAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTa                                                                       <  2:2227998/131‑1 (MQ=255)
                                                                                                 aCCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACtt                                                                                                           <  2:2052894/92‑1 (MQ=255)
                                                                                                 aCCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACtt                                                                                                           >  1:2052894/1‑92 (MQ=255)
                                                                                                 aCCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACttt                                                 <  1:1531163/150‑1 (MQ=255)
                                                                                                  cccACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTa                                                <  1:2227372/150‑1 (MQ=255)
                                                                                                   ccACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAg                                               >  1:1121699/1‑150 (MQ=255)
                                                                                                   ccACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAg                                               >  2:766527/1‑150 (MQ=255)
                                                                                                     aCGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCt                                                                                                                                               >  2:1406857/1‑52 (MQ=255)
                                                                                                     aCGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCt                                                                                                                                               <  1:1406857/52‑1 (MQ=255)
                                                                                                     aCGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGcc                                                                            >  2:665099/1‑119 (MQ=255)
                                                                                                     aCGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGcc                                                                            <  1:665099/119‑1 (MQ=255)
                                                                                                     aCGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGAGGCCCAGGCAACTTTAGCt                                             >  2:180315/1‑150 (MQ=255)
                                                                                                       gTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTAc                                           <  1:608865/150‑1 (MQ=255)
                                                                                                       gTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGAGGCCCAGGCAACTTTAGCTAc                                           <  1:180315/150‑1 (MQ=255)
                                                                                                            ggAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGcc                                                                                    <  2:513751/104‑1 (MQ=255)
                                                                                                            ggAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGcc                                                                                    >  1:513751/1‑104 (MQ=255)
                                                                                                            ggAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCACGCAAACCTGCTAAGCCGGCGTCTCAGGGAACTTTAGCTACgtgct                                      >  2:1345340/1‑150 (MQ=255)
                                                                                                              aaCACGGCCTGTACGTGCTGGAAGGCACTGCGTTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGt                                    >  1:1642866/1‑150 (MQ=255)
                                                                                                              aaCACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGt                                    >  1:179800/1‑150 (MQ=255)
                                                                                                               acacGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTa                                   >  1:135437/1‑150 (MQ=255)
                                                                                                               acacGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTa                                   <  1:766527/150‑1 (MQ=255)
                                                                                                                  cGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACaa                                <  2:1694782/150‑1 (MQ=255)
                                                                                                                  cGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACaa                                <  1:2122605/150‑1 (MQ=255)
                                                                                                                      cagTACGTGCTGGACCGCAACGCGGTTTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGAc                            <  2:1121699/148‑1 (MQ=255)
                                                                                                                      cTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCt                                                                                          >  2:2395708/1‑88 (MQ=255)
                                                                                                                      cTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCt                                                                                          <  1:2395708/88‑1 (MQ=255)
                                                                                                                        gTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGAGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGt                          <  2:1045237/150‑1 (MQ=255)
                                                                                                                         tACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTg                         >  2:681242/1‑150 (MQ=255)
                                                                                                                         tACGTGCTGGAAGGCAAGGCG‑‑GTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGaa                       <  2:1626609/150‑1 (MQ=255)
                                                                                                                           cGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGAc                            >  2:865439/1‑145 (MQ=255)
                                                                                                                           cGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGAc                            <  1:865439/145‑1 (MQ=255)
                                                                                                                           cGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGaa                       <  1:1160855/150‑1 (MQ=255)
                                                                                                                             tGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCcgg                                                                 >  2:1603102/1‑106 (MQ=255)
                                                                                                                             tGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCcgg                                                                 <  1:1603102/106‑1 (MQ=255)
                                                                                                                             tGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAAcc                     >  2:627585/1‑150 (MQ=255)
                                                                                                                               cTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCACACCAGTCACGTAAGCGTGCTACGCCGGCGGCGCAGGCACCTTTCGCTACCTGCTGTACAAGGACGTGAACCGt                   >  1:1251057/1‑150 (MQ=255)
                                                                                                                                tGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTc                  >  2:1159905/1‑150 (MQ=255)
                                                                                                                                 ggAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGa                             >  2:182790/1‑138 (MQ=255)
                                                                                                                                 ggAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGa                             <  1:182790/138‑1 (MQ=255)
                                                                                                                                 ggAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGAc                            <  2:903580/139‑1 (MQ=255)
                                                                                                                                 ggAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGAc                            >  1:903580/1‑139 (MQ=255)
                                                                                                                                 ggAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCa                 >  2:1967174/1‑150 (MQ=255)
                                                                                                                                    ggttcaggcGGTACACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGt              >  1:1822275/6‑150 (MQ=255)
                                                                                                                                    aggcaaggcGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGt              >  1:1170320/1‑150 (MQ=255)
                                                                                                                                     gttcaggcGGTACACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGAc                            <  1:4027/131‑1 (MQ=255)
                                                                                                                                     gttcaggcGGTACACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGAc                            >  2:4027/5‑135 (MQ=255)
                                                                                                                                      ttcaggcGGTACACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTac            >  1:1635558/4‑150 (MQ=255)
                                                                                                                                       caaggcGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTaca           <  2:186121/150‑1 (MQ=255)
                                                                                                                                            cGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGa      <  2:1507073/150‑1 (MQ=255)
                                                                                                                                                gtACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCt  <  2:919919/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GAAAACCACCCGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCT  >  NC_002947/4002388‑4002681

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