Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,293,725 T→C 6.6% D221G (GAT→GGT)  mmsA‑II ← methylmalonate‑semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,293,7250TC 6.6% 220.5 / 4.7 75D221G (GAT→GGT) mmsA‑IImethylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (37/33);  new base C (3/2);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01

AGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGCCATGG  >  NC_002947/5293580‑5293872
                                                                                                                                                 |                                                                                                                                                   
aGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATcc                                                                                                                                                   <  5:88856/148‑1 (MQ=255)
aGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGCGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATTGCATcc                                                                                                                                                   <  2:239570/148‑1 (MQ=255)
 ggCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATccac                                                                                                                                                 <  6:268530/149‑1 (MQ=255)
    agcaCCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATAcaccac                                                                                                                                              >  7:277153/1‑149 (MQ=255)
       accaccGCGTGGTTCTTGGCCCCCATCATTGACTGAACACTCTTGACATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGATGCCATTTGCACccacca                                                                                                                                               <  5:17579/145‑1 (MQ=255)
         caccGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCt                                                                                                                                         >  5:44847/1‑149 (MQ=255)
           ccGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTgc                                                                                                                                      <  4:286974/149‑1 (MQ=255)
             gcgTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTgcc                                                                                                                                     >  2:288080/1‑149 (MQ=255)
                    tCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTg                                                                                                                               >  3:206969/1‑148 (MQ=255)
                    tCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTg                                                                                                                               >  3:115953/1‑148 (MQ=255)
                      ttGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGcac                                                                                                                            >  4:127088/1‑149 (MQ=255)
                         gCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGcaccac                                                                                                                         >  1:123034/1‑149 (MQ=255)
                          cccccATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACg                                                                                                                        <  8:117746/149‑1 (MQ=255)
                            cccATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGtt                                                                                                                      >  3:96944/1‑149 (MQ=255)
                             ccATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTc                                                                                                                     <  1:288080/149‑1 (MQ=255)
                              catcatTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGtt                                                                                                                      >  4:97216/1‑147 (MQ=255)
                                   ttGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGc                                                                                                                                          >  3:69594/1‑122 (MQ=255)
                                   ttGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGc                                                                                                                                          <  4:69594/122‑1 (MQ=255)
                                            cacGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGccggccgg                                                                                                       >  4:68213/1‑148 (MQ=255)
                                            cacGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGt                                                                                                      <  8:136728/149‑1 (MQ=255)
                                               gCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGcacc                                                                                                                           <  7:208257/125‑1 (MQ=255)
                                               gCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGcacc                                                                                                                           >  8:208257/1‑125 (MQ=255)
                                               gCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGcacc                                                                                                                           >  8:204221/1‑125 (MQ=255)
                                               gCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGcacc                                                                                                                           <  7:204221/125‑1 (MQ=255)
                                               gCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTAc                                                                                                    >  2:206219/1‑148 (MQ=255)
                                               gCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTAc                                                                                                    <  1:206219/148‑1 (MQ=255)
                                               gCTTGCCATGCTGGCTGCCGAGGTTGTAAAAGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAAACGCCTTCATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACg                                                                                                   >  4:189074/1‑149 (MQ=255)
                                                    ccATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCgg                                                                                               <  4:115953/148‑1 (MQ=255)
                                                       tGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTc                                                                                           <  6:44847/149‑1 (MQ=255)
                                                          tggctgCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAg                                                                                        >  4:66379/1‑149 (MQ=255)
                                                           ggctgCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAgc                                                                                       <  4:47152/149‑1 (MQ=255)
                                                            gctgCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAgc                                                                                       <  4:242589/148‑1 (MQ=255)
                                                              tgCCGAGGTTGTAAACGTGCGTACCTACTTCAGTAGAGCCCACGAAGGACATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAgcgc                                                                                     >  5:169328/1‑148 (MQ=255)
                                                                        gTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGAcc                                                                           <  4:297481/148‑1 (MQ=255)
                                                                            aCGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGcaccac                                                                                                                         >  5:261219/1‑98 (MQ=255)
                                                                            aCGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGcaccac                                                                                                                         <  6:261219/98‑1 (MQ=255)
                                                                             cGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAg                                                                      >  1:265685/1‑148 (MQ=255)
                                                                                 cGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATg                                                                  <  3:66379/148‑1 (MQ=255)
                                                                                     ccGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGCAGCATGGTCGa                                                             >  6:254462/1‑149 (MQ=255)
                                                                                       gACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTGCAGCGCCAGCTCGACCAGAAGCATGGTCGa                                                             >  7:54388/1‑147 (MQ=255)
                                                                                        aCTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAg                                                          >  8:90699/1‑149 (MQ=255)
                                                                                              gTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATc                                                    >  6:161247/1‑149 (MQ=255)
                                                                                                    cccACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTc                                              >  1:278801/1‑149 (MQ=255)
                                                                                                        cGACGGAAATCGCCTTGATGTCCTGGGGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGAt                                          <  4:191819/149‑1 (MQ=255)
                                                                                                        cGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGAt                                          >  5:22718/1‑149 (MQ=255)
                                                                                                           aGGAAATCGCCTTGCTGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGCATGGt                                       <  7:279875/149‑1 (MQ=255)
                                                                                                                     cTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACg                              >  1:3030/1‑148 (MQ=255)
                                                                                                                        gATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAgg                          >  7:110671/1‑149 (MQ=255)
                                                                                                                         aTGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGgt                         >  1:60600/1‑149 (MQ=255)
                                                                                                                            tCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTGCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGtt                      <  8:54388/149‑1 (MQ=255)
                                                                                                                            tCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGtt                      <  4:96944/149‑1 (MQ=255)
                                                                                                                             ggtggtggATGCCATTTGCACCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTg                     >  1:288149/3‑149 (MQ=255)
                                                                                                                              gtggtggATGCCATTTGCACCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGc                    >  8:200036/2‑149 (MQ=255)
                                                                                                                              cTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCCTGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAg                           <  5:87136/142‑1 (MQ=255)
                                                                                                                              cTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAg                           >  6:87136/1‑142 (MQ=255)
                                                                                                                              cTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGc                    <  3:68213/149‑1 (MQ=255)
                                                                                                                                ggtggatgCCATTTGCACCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCg                  >  4:183559/1‑149 (MQ=39)
                                                                                                                                 gtgggtgCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGcc                   <  2:278801/147‑1 (MQ=255)
                                                                                                                                   gggtgCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCg                  >  6:76716/1‑146 (MQ=255)
                                                                                                                                      tgCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGa                                   >  6:121646/1‑126 (MQ=255)
                                                                                                                                      tgCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGa                                   <  5:121646/126‑1 (MQ=255)
                                                                                                                                       gCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGa           >  1:301629/1‑149 (MQ=255)
                                                                                                                                        cAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGAt          <  2:60600/149‑1 (MQ=255)
                                                                                                                                         aaaTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAg                           >  8:66614/1‑131 (MQ=255)
                                                                                                                                         aaaTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAg                           <  7:66614/131‑1 (MQ=255)
                                                                                                                                         aaaTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATg         >  3:179733/1‑149 (MQ=255)
                                                                                                                                          aaTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGc                                 <  1:154572/124‑1 (MQ=255)
                                                                                                                                          aaTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGc                                 >  2:154572/1‑124 (MQ=255)
                                                                                                                                           aTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGc       >  5:203258/1‑149 (MQ=255)
                                                                                                                                             tgCACCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGcc      <  2:288149/147‑1 (MQ=255)
                                                                                                                                             ggCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGCCa     >  8:39365/1‑149 (MQ=255)
                                                                                                                                              gCATCCACCACCTGCTTGGCGCCGTGCACGACGTTCAGCACGCCTGCCTGTACGCCGTCTTCCAGCGCCAGCTCGACCAGAATCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACTAAGGTGTTGCCTCAGACGATGGCCAt    <  4:179733/149‑1 (MQ=255)
                                                                                                                                               cATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGCCATg   <  8:110671/149‑1 (MQ=255)
                                                                                                                                                aTCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGCCATgg  <  3:97216/149‑1 (MQ=255)
                                                                                                                                                aTCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGCCATgg  <  6:169328/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGCCATGG  >  NC_002947/5293580‑5293872

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: