Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,212,955 T→G 6.2% K253Q (AAG→CAG)  PP_4594 ← cystathionine gamma‑synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,212,9550TG 6.2% 206.9 / 4.7 65K253Q (AAG→CAG) PP_4594cystathionine gamma‑synthase
Reads supporting (aligned to +/- strand):  ref base T (30/31);  new base G (2/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

CTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGCCACCACCACACCGCCGGTGAGGTCGGAGTGGCCACCGATGTACTTGCT  >  NC_002947/5212807‑5213102
                                                                                                                                                    |                                                                                                                                                   
cTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCtt                                                                                                                                                     <  6:187628/149‑1 (MQ=255)
   cgcATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCtt                                                                                                                                                     >  3:23228/1‑146 (MQ=255)
              tttGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCTCACATCCAGGGTCTGAAGAccctggat                                                                                                                                         <  1:17463/147‑7 (MQ=255)
              tttGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTGAAGAccctggat                                                                                                                                         <  1:297535/147‑7 (MQ=255)
              tttGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTGAAGAccctggat                                                                                                                                         >  2:17463/1‑141 (MQ=255)
              tttGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTGAAGAccctggat                                                                                                                                         >  2:297535/1‑141 (MQ=255)
                   acaGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAgg                                                                                                                                  >  1:288586/1‑149 (MQ=255)
                   acaGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAgg                                                                                                                                  <  1:100161/149‑1 (MQ=255)
                     aGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAgg                                                                                                                                  >  5:219206/1‑147 (MQ=255)
                      gTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAg                                                                                                                               <  6:146376/149‑1 (MQ=255)
                        tCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAg                                                                                                                               >  5:77859/1‑147 (MQ=255)
                            gctgcGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTc                                                                                                                          >  6:134344/1‑148 (MQ=255)
                            gctgcGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTc                                                                                                                          >  5:6282/1‑148 (MQ=255)
                                cGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGaa                                                                                                                       <  5:204073/147‑1 (MQ=255)
                                    tGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGccc                                                                                                                 <  4:23228/149‑1 (MQ=255)
                                      gCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAgg                                                                                                                                        >  1:265798/1‑124 (MQ=255)
                                      gCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAgg                                                                                                                                        <  2:265798/124‑1 (MQ=255)
                                      gCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCt                                                                                                               >  5:62218/1‑149 (MQ=255)
                                       cTTTCCAGACCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTg                                                                                                              >  2:233719/1‑149 (MQ=255)
                                        tttCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTg                                                                                                              >  7:38367/1‑148 (MQ=255)
                                           ccAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTgcac                                                                                                           <  6:62218/148‑1 (MQ=255)
                                            cAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTgcaccg                                                                                                         >  6:261482/1‑149 (MQ=255)
                                             aGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTgcaccgc                                                                                                        <  2:288586/149‑1 (MQ=255)
                                               ccccGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTgcaccgcac                                                                                                      >  7:117627/1‑149 (MQ=255)
                                                cccGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTgcaccgcacc                                                                                                     <  6:77859/149‑1 (MQ=255)
                                                  cGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGCGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCa                                                                                                    >  4:258700/1‑148 (MQ=255)
                                                    gctAATACACCTGCTCGACCTGCGCGAGACCTTCCAGGAAGCGCGCCACCTGCAGCGCATTGGCGCACTGGCGCTCCATGCGCCCATCCACGGTCATCAGACCGCGCACGGCCAGGGAGCAGTCGAACGGCCCCTCCAAAGCACCAATc                                                                                                 >  8:158907/3‑149 (MQ=255)
                                                         tACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCAt                                                                                            >  4:206280/1‑149 (MQ=255)
                                                         tACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCAt                                                                                            <  5:261482/149‑1 (MQ=255)
                                                           cacCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGc                                                                                          <  3:206280/149‑1 (MQ=255)
                                                           cacCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGc                                                                                          <  6:219206/149‑1 (MQ=255)
                                                            acCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCt                                                                                         >  3:77760/1‑149 (MQ=255)
                                                              cTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGa                                                                                       >  4:260398/1‑149 (MQ=255)
                                                                 cTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCTGACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATa                                                                                     <  7:169239/148‑1 (MQ=255)
                                                                      ccTGCGGGTTACCTTCCACGATTCGCTGCACCTTCATGTCATTTGCGCACTGGCGCTCCATTCGCAAATCCAGGGTCTTCATACCGCGCATGGCCATGTATCATTCGAACTGCCCTTGCACCGCACCAATCGCCATGCTGATACGgcgc                                                                               <  8:239762/149‑1 (MQ=255)
                                                                          cgcgTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAgg                                                                            <  2:364892/148‑1 (MQ=255)
                                                                               gACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCa                                                                                                                             <  5:100198/94‑1 (MQ=255)
                                                                               gACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCa                                                                                                                             >  6:100198/1‑94 (MQ=255)
                                                                                aCCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCa                                                                     <  6:211648/149‑1 (MQ=255)
                                                                                  cTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGt                                                                         <  5:116369/143‑1 (MQ=255)
                                                                                  cTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGt                                                                         >  6:116369/1‑143 (MQ=255)
                                                                                      cAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTgc                                                                                                             <  1:380559/103‑1 (MQ=255)
                                                                                      cAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTgc                                                                                                             >  2:380559/1‑103 (MQ=255)
                                                                                       aGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACg                                                                                   <  7:306073/128‑1 (MQ=255)
                                                                                       aGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACg                                                                                   >  8:306073/1‑128 (MQ=255)
                                                                                           aGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCg                                                          >  5:317889/1‑149 (MQ=255)
                                                                                             cgcgcCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGtt                                                        <  6:317889/149‑1 (MQ=255)
                                                                                               cgcCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTgg                                                      <  4:77760/149‑1 (MQ=255)
                                                                                                gcCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGgc                                                     <  8:117627/149‑1 (MQ=255)
                                                                                                       gCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGccacc                                              >  7:303928/1‑149 (MQ=255)
                                                                                                       gCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGccac                                               >  1:54443/1‑148 (MQ=255)
                                                                                                        cAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGgcc                                                    >  7:281369/1‑142 (MQ=255)
                                                                                                        cAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGgcc                                                    <  8:281369/142‑1 (MQ=255)
                                                                                                         aGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGccaccac                                            <  6:71208/149‑1 (MQ=255)
                                                                                                         aGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGccaccac                                            <  5:134344/149‑1 (MQ=255)
                                                                                                             cATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAAGGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGCCACCACcaca                                        >  5:175747/1‑149 (MQ=255)
                                                                                                                      aCTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGCCACCACCACACCGCCGGTg                               <  8:38367/149‑1 (MQ=255)
                                                                                                                      aCTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGCCACCACCACACCGCCGGTg                               <  2:54443/149‑1 (MQ=255)
                                                                                                                        tGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGCCACCACCACACCGCCGGTg                               <  7:199127/147‑1 (MQ=255)
                                                                                                                        tGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGCCACCACCACACCGCCGGTg                               >  8:199127/1‑147 (MQ=255)
                                                                                                                          gcgcTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGCCACCACCACACCGCCGGTGAGGt                           <  2:122316/149‑1 (MQ=255)
                                                                                                                                cATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGCCACCACCACACCGCCGGTGAGGTCGGa                       <  3:260398/147‑1 (MQ=255)
                                                                                                                                     gcACATCCAGGGTCTTCAGACCGCACAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATAGTTGGCCGCCACCACCACACCGCCGGTGAGGTCGGAGTGGccc                >  7:239371/1‑148 (MQ=255)
                                                                                                                                       acaTCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGCCACCACCACACCGCCGGTGAGGTCGGAGTGGCCAcc              >  6:359266/1‑149 (MQ=255)
                                                                                                                                        caTCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGCCACCACCACACCGCCGGTGAGGTCGGAGTGGCCACCg             <  6:342120/149‑1 (MQ=255)
                                                                                                                                           ccAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGCCACCACCACACCGCCGGTGAGGTCGGAGTGGCCACCGATg          >  3:84986/1‑149 (MQ=255)
                                                                                                                                                  cTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGCCACCACCACACCGCCGGTGAGGTCGGAGTGGCCACCGATGTACTTGc   >  1:320868/1‑149 (MQ=255)
                                                                                                                                                   ttCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATg                                                                                           >  7:124662/1‑60 (MQ=255)
                                                                                                                                                   ttCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATg                                                                                           <  8:124662/60‑1 (MQ=255)
                                                                                                                                                   ttCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGCCACCACCACACCGCCGGTGAGGTCGGAGTGGCCACCGATGTACTTGCt  >  5:382354/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTGCGCATCTGGCGTTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGCCACCACCACACCGCCGGTGAGGTCGGAGTGGCCACCGATGTACTTGCT  >  NC_002947/5212807‑5213102

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: ‑