Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,212,956 C→A 5.9% L252L (CTG→CTT PP_4594 ← cystathionine gamma‑synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,212,9560CA 5.9% 218.8 / 3.7 68L252L (CTG→CTTPP_4594cystathionine gamma‑synthase
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base A (2/2);  total (34/34)
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.11e-01

TTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGCCACCACCACACCGCCGGTGAGGTCGGAGTGGCCACCGATGTACTTGCTGG  >  NC_002947/5212821‑5213104
                                                                                                                                       |                                                                                                                                                    
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                                                                                                                            >  5:6282/1‑148 (MQ=255)
              gctgcGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTc                                                                                                                            >  6:134344/1‑148 (MQ=255)
                  cGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGaa                                                                                                                         <  5:204073/147‑1 (MQ=255)
                      tGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGccc                                                                                                                   <  4:23228/149‑1 (MQ=255)
                        gCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAgg                                                                                                                                          <  2:265798/124‑1 (MQ=255)
                        gCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAgg                                                                                                                                          >  1:265798/1‑124 (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                                                                                            <  6:219206/149‑1 (MQ=255)
                                             cacCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGc                                                                                            <  3:206280/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                                                      <  8:281369/142‑1 (MQ=255)
                                                                                          cAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGgcc                                                      >  7:281369/1‑142 (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)
                                                                                                                                       cAGACCGCGCAGGGCCAGGTAGCAGTCGAACGTCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCCACGCATCGTTGGCCGCCACCACCACACCGCCGGTGAGGGCGGAGTGGCCACCGATGTACTTGCTgg  >  4:222224/1‑149 (MQ=255)
                                                                                                                                       cAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGCCACCACCACACCGCCGGTGAGGTCGGAGTGGCCACCGATGTACTTGCt    <  5:359266/147‑1 (MQ=255)
                                                                                                                                       cAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGCCACCACCACACCGCCGGTGAGGTCGGAGTGGCCACCGATGTACTTGCTg   <  6:6282/148‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                    
TTTGCACAGTTCATGCTGCGGGTGGCTTTCCAGCCCCGGGTAATACACCTGCTCGACCTGCGCGTGACCTTCCAGGAAGCGCGCCACCTGCAGGGCATTGGCGCACTGGCGCTCCATGCGCACATCCAGGGTCTTCAGACCGCGCAGGGCCAGGTAGCAGTCGAACGGCCCCTGCACCGCACCAATCGCCATGCTGATACGGCGCAGGCGTGCCAGCAACGCATCGTTGGCCGCCACCACCACACCGCCGGTGAGGTCGGAGTGGCCACCGATGTACTTGCTGG  >  NC_002947/5212821‑5213104

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: