Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,293,721 G→T 6.6% 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 6.6% 222.4 / 3.9 70A222A (GCC→GCAmmsA‑IImethylmalonate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (37/28);  new base T (3/2);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

CATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGCCA  >  NC_002947/5293576‑5293869
                                                                                                                                                 |                                                                                                                                                    
cATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCa                                                                                                                                                   >  7:136728/1‑149 (MQ=255)
    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:115953/1‑148 (MQ=255)
                        tCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTg                                                                                                                            >  3:206969/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                                                                                                                                       <  4:69594/122‑1 (MQ=255)
                                       ttGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGc                                                                                                                                       >  3:69594/1‑122 (MQ=255)
                                                cacGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGccggccgg                                                                                                    >  4:68213/1‑148 (MQ=255)
                                                cacGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGt                                                                                                   <  8:136728/149‑1 (MQ=255)
                                                   gCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGcacc                                                                                                                        >  8:208257/1‑125 (MQ=255)
                                                   gCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGcacc                                                                                                                        <  7:208257/125‑1 (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                                                                                                                      <  6:261219/98‑1 (MQ=255)
                                                                                aCGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGcaccac                                                                                                                      >  5:261219/1‑98 (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                        <  7:66614/131‑1 (MQ=255)
                                                                                                                                             aaaTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAg                        >  8:66614/1‑131 (MQ=255)
                                                                                                                                             aaaTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATg      >  3:179733/1‑149 (MQ=255)
                                                                                                                                              aaTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGc                              >  2:154572/1‑124 (MQ=255)
                                                                                                                                              aaTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGc                              <  1:154572/124‑1 (MQ=255)
                                                                                                                                               aTGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGc    >  5:203258/1‑149 (MQ=255)
                                                                                                                                                 ggCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGCCa  >  8:39365/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CATCAGGCAGCACCACCGCGTGGTTCTTGGCCCCCATCATTGACTGCACACGCTTGCCATGCTGGCTGCCGAGGTTGTAAACGTGCGTACCGACTTCAGTAGAGCCCACGAAGGAAATCGCCTTGATGTCCTGGTGGGTGCAAATGGCATCCACCACCTGCTTGCCGCCGTGCACCACGTTCAGCACGCCGGCCGGTACGCCGGCTTCCAGCGCCAGCTCGACCAGAAGCATGGTCGACAGCGGATCCTGCTCGGATGGTTTGAGCACGAAGGTGTTGCCGCAGACGATGGCCA  >  NC_002947/5293576‑5293869

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: