Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,293,717 A→C 8.1% C224G (TGC→GGC)  mmsA‑II ← methylmalonate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,293,7170AC 8.1% 224.0 / 9.9 74C224G (TGC→GGC) mmsA‑IImethylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (34/34);  new base C (3/3);  total (37/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.42e-01

GTTGGCATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATG  >  NC_002947/5293571‑5293865
                                                                                                                                                  |                                                                                                                                                    
gTTGGCATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCaaa                                                                                                                                                    <  5:197837/149‑1 (MQ=255)
 ttGGCATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAAt                                                                                                                                                   >  2:337498/1‑149 (MQ=255)
   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                                                                                                                                 <  2:139016/149‑1 (MQ=255)
                   accGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCtt                                                                                                                                 <  7:128791/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                                                 >  6:137325/1‑149 (MQ=255)
                                                                                                   ttCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCg                                                 >  5:324415/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                <  2:298888/145‑1 (MQ=255)
                                                                                                                                      ggtggtggATGCCATTTGCACCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGgtgt                >  1:298888/3‑147 (MQ=255)
                                                                                                                                      ggtggtggATGCCATTTGCACCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTg              >  5:329432/3‑149 (MQ=255)
                                                                                                                                       cTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGa                            >  6:302292/1‑134 (MQ=255)
                                                                                                                                       cTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGa                            <  5:302292/134‑1 (MQ=255)
                                                                                                                                       cTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGc             >  4:98416/1‑149 (MQ=255)
                                                                                                                                       cTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGc             <  6:79419/149‑1 (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)
                                                                                                                                                  |                                                                                                                                                    
GTTGGCATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATG  >  NC_002947/5293571‑5293865

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: