Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,002,532 G→A 8.8% V224V (GTG→GTA allE → S‑ureidoglycine aminohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,002,5320GA 8.8% 278.5 / 13.4 103V224V (GTG→GTAallES‑ureidoglycine aminohydrolase
Reads supporting (aligned to +/- strand):  ref base G (54/39);  new base A (4/5);  total (58/45)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

Only 100 of 107 total aligned reads displayed.

CGCTGGAAAACCACCCGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCT  >  NC_002947/4002383‑4002681
                                                                                                                                                     |                                                                                                                                                     
cGCTGGAAAACCACCCGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGgtg                                                                                                                                                       <  1:875310/150‑1 (MQ=255)
     gAAAACCACCCGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGTGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCg                                                                                                                                                  >  1:1373715/1‑150 (MQ=255)
     gAAAACCACCCGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCg                                                                                                                                                  >  2:1023474/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)
          ccaccCGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTACACCGCCttg                                                                                                                                             <  1:2356879/150‑3 (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)
                    tGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATGCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTg                                                                                                                                   <  1:1620878/150‑1 (MQ=255)
                     gTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTg                                                                                                                                   >  2:647207/1‑149 (MQ=255)
                     gTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTg                                                                                                                                   <  1:647207/149‑1 (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)
                           aTGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGaa                                                                                                                            <  2:31971/150‑1 (MQ=255)
                                         ccACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGaa                                                                                                                      >  1:943576/1‑142 (MQ=255)
                                         ccACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGa                                                                                                              >  2:1904312/1‑150 (MQ=255)
                                                 tGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCAtgt                                                                                                      >  1:1045237/1‑150 (MQ=255)
                                                      gTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTg                                                                                                 >  2:1550715/1‑150 (MQ=255)
                                                      gTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTg                                                                                                 >  2:2227372/1‑150 (MQ=255)
                                                      gTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTg                                                                                                 >  1:1694782/1‑150 (MQ=255)
                                                           cATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTACACCGCCtt                                                                                                                                              <  1:1770930/100‑2 (MQ=255)
                                                            aTCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGcc                                                                                           >  2:2368457/1‑150 (MQ=255)
                                                                tCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTgg                                                                                                                                  <  1:1298458/107‑1 (MQ=255)
                                                                tCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCt                                                                                       >  1:1981701/1‑150 (MQ=255)
                                                                tCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCAGGGTACCCCCTCAACCAGTACTCGGTCGAAGTGGAAGCCGGCCACTTCCTGCGGCTGCCCGCAttgc                                                                                       >  2:747344/1‑148 (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                                                                       <  2:2227998/131‑1 (MQ=255)
                                                                                                   gAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTa                                                                       >  1:2227998/1‑131 (MQ=255)
                                                                                                      aCCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACtt                                                                                                           <  2:2052894/92‑1 (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                                                                                                                                               <  1:1406857/52‑1 (MQ=255)
                                                                                                          aCGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCt                                                                                                                                               >  2:1406857/1‑52 (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                                                                                    >  1:513751/1‑104 (MQ=255)
                                                                                                                 ggAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGcc                                                                                    <  2:513751/104‑1 (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:766527/150‑1 (MQ=255)
                                                                                                                    acacGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTa                                   >  1:135437/1‑150 (MQ=255)
                                                                                                                       cGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACaa                                <  2:1694782/150‑1 (MQ=255)
                                                                                                                       cGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACaa                                <  1:2122605/150‑1 (MQ=255)
                                                                                                                           cagTACGTGCTGGACCGCAACGCGGTTTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGAc                            <  2:1121699/148‑1 (MQ=255)
                                                                                                                           cTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCt                                                                                          <  1:2395708/88‑1 (MQ=255)
                                                                                                                           cTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCt                                                                                          >  2:2395708/1‑88 (MQ=255)
                                                                                                                             gTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGt                          <  2:1642866/150‑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)
                                                                                                                                  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                            >  1:903580/1‑139 (MQ=255)
                                                                                                                                      ggAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGAc                            <  2:903580/139‑1 (MQ=255)
                                                                                                                                      ggAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCa                 >  2:1967174/1‑150 (MQ=255)
                                                                                                                                         ggttcaggcGGTACACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGt              >  1:1822275/6‑150 (MQ=255)
                                                                                                                                         aggcaaggcGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGt              >  1:1170320/1‑150 (MQ=255)
                                                                                                                                          gttcaggcGGTACACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGAc                            >  2:4027/5‑135 (MQ=255)
                                                                                                                                          gttcaggcGGTACACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGAc                            <  1:4027/131‑1 (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)
                                                                                                                                                     |                                                                                                                                                     
CGCTGGAAAACCACCCGTTTTGTCGACATGGCCGACATGCGCCACGACATGCACGTGAACATCGTCACCTTCCAGCCGGGCGGTGTGATTCCGTTCGCAGAAACCCACGTGATGGAACACGGCCTGTACGTGCTGGAAGGCAAGGCGGTGTACCGCCTGAACCAGGACTGGGTCGAAGTGGAAGCCGGCGACTTCATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCT  >  NC_002947/4002383‑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: