Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,293,721 G→T 7.5% A222A (GCC→GCA mmsA‑II ← methylmalonate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,293,7210GT 7.5% 229.4 / 6.1 76A222A (GCC→GCAmmsA‑IImethylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (36/34);  new base T (3/3);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

GGCATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGC  >  NC_002947/5293574‑5293867
                                                                                                                                                   |                                                                                                                                                  
ggCATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATg                                                                                                                                                    <  5:147736/148‑1 (MQ=255)
  cATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCa                                                                                                                                                 >  6:160571/1‑149 (MQ=255)
        gcagcaCCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATccac                                                                                                                                            >  1:196084/1‑148 (MQ=255)
        gcagcaCCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCAACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATccac                                                                                                                                            >  5:79419/1‑148 (MQ=255)
          agcaCCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCcaccac                                                                                                                                         >  3:346502/1‑149 (MQ=255)
          agcaCCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGATGCCATTTGCACccaccag                                                                                                                                         <  5:343883/149‑2 (MQ=255)
          agcaCCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGATGCCATTTGCACccac                                                                                                                                            <  3:172106/146‑1 (MQ=255)
          agcaCCACCGCGTGGTTCTTGGCCCCCATAATTGACTGCACACGCTTGCCATGCTGGCTGCGGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGATGCCATTTGCACccac                                                                                                                                            >  4:172106/1‑146 (MQ=255)
                accGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCtt                                                                                                                                   <  7:128791/149‑1 (MQ=255)
                accGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCtt                                                                                                                                   <  2:139016/149‑1 (MQ=255)
                accGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCtt                                                                                                                                   <  3:267913/149‑1 (MQ=255)
                accGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCtt                                                                                                                                   <  3:1828/149‑1 (MQ=255)
                  cgcgTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTgc                                                                                                                                 >  4:195218/1‑149 (MQ=255)
                   gcgTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTgcc                                                                                                                                >  5:306704/1‑149 (MQ=255)
                       ggTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGccgccg                                                                                                                            >  4:83570/1‑149 (MQ=255)
                        gTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGt                                                                                                                           <  1:119566/149‑1 (MQ=255)
                              ggCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGcacca                                                                                                                     >  2:241042/1‑149 (MQ=255)
                               gCCCCCCTCCTTCACCGCCCCCGCTGGCCATGCTGGCTCCCGAGGTTGTAAACGTGCGTACCGACTTCACTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGGGCAAATGGCATCCACCACCTGCTTGCCGCCGTGcaccac                                                                                                                    <  4:230419/149‑1 (MQ=255)
                                cccccATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACg                                                                                                                   >  5:96197/1‑149 (MQ=255)
                                 ccccATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGt                                                                                                                  <  5:286071/149‑1 (MQ=255)
                                 ccccATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGt                                                                                                                  >  4:336238/1‑149 (MQ=255)
                                 ccccATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGt                                                                                                                  <  4:25545/149‑1 (MQ=255)
                                     atcatTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAg                                                                                                              >  6:272138/1‑149 (MQ=255)
                                     atcatTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAg                                                                                                              <  4:346502/149‑1 (MQ=255)
                                      tcatTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGc                                                                                                             <  8:268649/149‑1 (MQ=255)
                                         ttGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACg                                                                                                          >  4:210011/1‑149 (MQ=255)
                                            aCTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACgccg                                                                                                       <  6:156521/149‑1 (MQ=255)
                                                cacaCGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAATTCGCCTTGATGTCCTGGTGGGTGCAAATGGCATGGACCACCTGCTTGCCGCCTTGCACCACGTTCAGCACgccggc                                                                                                     <  4:153954/147‑1 (MQ=255)
                                                cacaCGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTCCAGCACgccggccg                                                                                                   >  4:123831/1‑149 (MQ=255)
                                                   acGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTa                                                                                                >  3:259126/1‑149 (MQ=255)
                                                     gCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACg                                                                                              >  2:319764/1‑149 (MQ=255)
                                                       ttGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGcc                                                                                            >  2:84303/1‑149 (MQ=255)
                                                        tGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCg                                                                                           >  7:342957/1‑149 (MQ=255)
                                                            aTGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCtt                                                                                       <  8:342957/149‑1 (MQ=255)
                                                              gctggctgCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTACACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTc                                                                                      >  7:214679/1‑148 (MQ=255)
                                                                      ccGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCa                                                                              <  3:195218/148‑1 (MQ=255)
                                                                        gAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCt                                                                           <  7:55176/149‑1 (MQ=255)
                                                                            ttGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGAc                                                                       <  1:241042/149‑1 (MQ=255)
                                                                                 aaCGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGa                                                                   <  2:207910/148‑1 (MQ=255)
                                                                                 aaCGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTCCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGa                                                                   <  3:123831/148‑1 (MQ=255)
                                                                                     tGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCAt                                                              <  1:337498/149‑1 (MQ=255)
                                                                                         tACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTc                                                          <  2:196084/149‑1 (MQ=255)
                                                                                         tACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTc                                                          <  1:33199/149‑1 (MQ=255)
                                                                                         tACCGACTTCAGTAGAGCCAACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTc                                                          >  8:65720/1‑149 (MQ=255)
                                                                                           ccGACTTCACTAGAGCCAACGAAGCAAATCGCCTTGATGTCCTGGTGGGTGCAAAAGGCATCCACCACATGCTTGAAGCCCAGCACAAAGTTCACAACGCCGGCCGGTACGACGGATTCAAGAGCCAGCTCGACCAGAACCATGTTCg                                                         >  4:362792/1‑148 (MQ=255)
                                                                                                ttCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGaa                                                                  >  6:150375/1‑134 (MQ=255)
                                                                                                ttCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGaa                                                                  <  5:150375/134‑1 (MQ=255)
                                                                                                ttCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCg                                                   >  5:324415/1‑149 (MQ=255)
                                                                                                ttCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCg                                                   >  6:137325/1‑149 (MQ=255)
                                                                                                   aGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGa                                                 <  5:160571/148‑1 (MQ=255)
                                                                                                            cACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCg                                        >  1:89059/1‑148 (MQ=255)
                                                                                                                   gAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACa                                                      <  8:164588/127‑1 (MQ=255)
                                                                                                                   gAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACa                                                      >  7:164588/1‑127 (MQ=255)
                                                                                                                     aaTCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCCCGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGa                              <  7:177757/149‑1 (MQ=255)
                                                                                                                     aaTCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGttt                                >  1:109868/1‑147 (MQ=255)
                                                                                                                       tCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGc                            <  3:336238/149‑1 (MQ=255)
                                                                                                                           cTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGa                        <  3:83570/149‑1 (MQ=255)
                                                                                                                               aTGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAgg                     <  6:306704/148‑1 (MQ=255)
                                                                                                                               aTGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGgt                    <  6:317656/149‑1 (MQ=255)
                                                                                                                                 gTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGa                                      >  7:190717/1‑129 (MQ=255)
                                                                                                                                 gTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGa                                      <  8:190717/129‑1 (MQ=255)
                                                                                                                                  tCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGtt                 >  5:64533/1‑149 (MQ=255)
                                                                                                                                   ggtggtggATGCCATTTGCACCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGgtgt                  >  1:298888/3‑147 (MQ=255)
                                                                                                                                   ggtggtggATGCCATTTGCACCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGgtgt                  <  2:298888/145‑1 (MQ=255)
                                                                                                                                   ggtggtggATGCCATTTGCACCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTg                >  5:329432/3‑149 (MQ=255)
                                                                                                                                    cTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGa                              <  5:302292/134‑1 (MQ=255)
                                                                                                                                    cTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGa                              >  6:302292/1‑134 (MQ=255)
                                                                                                                                    cTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGc               <  6:79419/149‑1 (MQ=255)
                                                                                                                                    cTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGc               >  4:98416/1‑149 (MQ=255)
                                                                                                                                    cTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTACACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGc               <  8:214679/149‑1 (MQ=255)
                                                                                                                                          ggtgCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGa         >  7:342403/1‑149 (MQ=255)
                                                                                                                                            tgCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCTGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGAc        >  6:61839/1‑148 (MQ=255)
                                                                                                                                            tgCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGAc        >  5:374351/1‑148 (MQ=255)
                                                                                                                                               aaaTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATg    <  5:272138/149‑1 (MQ=255)
                                                                                                                                               aaaTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATg    >  3:154060/1‑149 (MQ=255)
                                                                                                                                                aaTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATgg   >  6:128207/1‑149 (MQ=255)
                                                                                                                                                 aTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGc  >  3:325865/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGCATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGC  >  NC_002947/5293574‑5293867

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