Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,550,907 C→G 7.2% T92T (ACC→ACG pydA → NADP‑dependent dihydropyrimidine dehydrogenase subunit PreA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,550,9070CG 7.2% 214.3 / 6.4 69T92T (ACC→ACGpydANADP‑dependent dihydropyrimidine dehydrogenase subunit PreA
Reads supporting (aligned to +/- strand):  ref base C (35/29);  new base G (3/2);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

GAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCC  >  NC_002947/4550763‑4551043
                                                                                                                                                |                                                                                                                                          
gAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAgg                                                                                                                                        >  2:227614/1‑149 (MQ=255)
gAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGCTACAGGGGATCAACAATATCGAGTTGATCACCGCCCGCTCTATGGAAATCAACCTGCGTGAAATCACCCAgg                                                                                                                                        >  2:190310/1‑149 (MQ=255)
    aCCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTgaa                                                                                                                                   <  2:76794/149‑1 (MQ=255)
      ccTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCAGGAAATCAACCTGCGTGAAATCACCCAGGTgaag                                                                                                                                  >  4:280944/1‑148 (MQ=255)
          ggCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGa                                                                                                                             >  5:141894/1‑149 (MQ=255)
              aaGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTgg                                                                                                                         <  7:130414/149‑1 (MQ=255)
              aaGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCGCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTgg                                                                                                                         >  2:271674/1‑149 (MQ=255)
               aGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGc                                                                                                                        <  6:286890/149‑1 (MQ=255)
                gACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGcc                                                                                                                       >  8:297670/1‑149 (MQ=255)
                   ccTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGa                                                                                                                    >  2:227187/1‑149 (MQ=255)
                        ggcggTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACcgcg                                                                                                               <  4:38201/149‑1 (MQ=255)
                          cggTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCa                                                                                                             <  1:227614/149‑1 (MQ=255)
                               aaCGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACACGGGAACAACAATATCGAGTTGATCACCGACCGCTCACTGGAAATCAACCTGCCAGAAATCACCCAGGTGAAGAAGGACAGGCCCGACCGCGCACTCAt                                                                                                        >  2:299993/1‑149 (MQ=255)
                                  gtgtCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCgt                                                                                                     <  1:227187/149‑1 (MQ=255)
                                          gcgcTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGa                                                                                             <  1:202813/149‑1 (MQ=255)
                                          gcgcTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGa                                                                                             <  7:85321/149‑1 (MQ=255)
                                            gcTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATg                                                                                           >  7:282892/1‑149 (MQ=255)
                                                cTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGt                                                                                         <  6:262783/147‑1 (MQ=255)
                                                cTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGc                                                                                       <  6:141894/149‑1 (MQ=255)
                                                      gcACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGCAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGc                                                                                  <  1:258605/149‑1 (MQ=255)
                                                      gcACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCa                                                                                 >  1:64519/1‑149 (MQ=255)
                                                        aCTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATa                                                                               <  2:279726/149‑1 (MQ=255)
                                                            tGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATaga                                                                             <  2:64519/147‑1 (MQ=255)
                                                               ccccAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAAt                                                                        >  2:208248/1‑149 (MQ=255)
                                                                 ccAACCGTCTGGTACAGGGGATCAACAATATCGAGTTTATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCg                                                                     <  1:190310/149‑1 (MQ=255)
                                                                     ccGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGa                                                                 <  7:297670/149‑1 (MQ=255)
                                                                       gTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAg                                                               <  6:168573/149‑1 (MQ=255)
                                                                       gTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAg                                                               >  5:168573/1‑149 (MQ=255)
                                                                       gTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAg                                                               <  1:208248/149‑1 (MQ=255)
                                                                        tCAGGTACAGGGGATCAACAATGTCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGt                                                              >  5:316237/1‑149 (MQ=255)
                                                                               cAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCt                                                       >  3:106185/1‑149 (MQ=255)
                                                                               cAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCt                                                       >  3:311640/1‑149 (MQ=255)
                                                                                  gggATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCAc                                                                                                            >  5:6728/1‑94 (MQ=255)
                                                                                  gggATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCAc                                                                                                            <  6:6728/94‑1 (MQ=255)
                                                                                         acaaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATagag                                                                            <  4:139506/119‑1 (MQ=255)
                                                                                         acaaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATagag                                                                            >  3:139506/1‑119 (MQ=255)
                                                                                         acaaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtggtg                                             >  1:116303/1‑149 (MQ=255)
                                                                                         acaaTATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtggtg                                             >  6:141806/1‑149 (MQ=255)
                                                                                               tCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGc                                        >  7:205333/1‑148 (MQ=255)
                                                                                               tCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGc                                        >  5:162645/1‑148 (MQ=255)
                                                                                               tCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGc                                        >  4:201525/1‑148 (MQ=255)
                                                                                                cGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCa                                      >  6:179171/1‑149 (MQ=255)
                                                                                                        tCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGcc                              >  1:253704/1‑149 (MQ=255)
                                                                                                        tCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGcc                              >  8:176212/1‑149 (MQ=255)
                                                                                                            cgaccgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAATAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATg                          <  4:32209/149‑1 (MQ=255)
                                                                                                                cgCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCAt                      <  8:282892/149‑1 (MQ=255)
                                                                                                                 gCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACggg                                  >  5:324319/1‑136 (MQ=255)
                                                                                                                 gCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATc                     >  1:317078/1‑149 (MQ=255)
                                                                                                                 gCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGGGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACggg                                  <  6:324319/136‑1 (MQ=255)
                                                                                                                   tctctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGAGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATc                     >  4:123969/1‑147 (MQ=255)
                                                                                                                     tctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGa                   <  1:299993/147‑1 (MQ=255)
                                                                                                                     tctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGc                 >  5:181722/1‑149 (MQ=255)
                                                                                                                     tctGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGc                 >  2:102814/1‑149 (MQ=255)
                                                                                                                       tGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGc                 <  7:176212/147‑1 (MQ=255)
                                                                                                                       tGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCt                >  6:274632/1‑148 (MQ=255)
                                                                                                                       tGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCt                >  7:125412/1‑148 (MQ=255)
                                                                                                                          ttcttcaCCTGGGTGATTTCACGCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAAc            >  4:236150/7‑149 (MQ=255)
                                                                                                                            aTCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCCCTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACt           <  1:271674/149‑1 (MQ=255)
                                                                                                                               caCCTGGGTGATTTCACGCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCgg       >  5:37094/2‑149 (MQ=255)
                                                                                                                               caCCTGGGTGATTTCACGCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCgg       >  4:91447/2‑149 (MQ=255)
                                                                                                                                aCCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGt      <  6:181722/149‑1 (MQ=255)
                                                                                                                                 ccTGGGTGATTTCACGCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGtt     <  6:37094/149‑1 (MQ=255)
                                                                                                                                 ccTGGGTGATTTCACGCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGtt     <  3:236150/149‑1 (MQ=255)
                                                                                                                                  cTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtggt                                              >  6:205838/1‑107 (MQ=255)
                                                                                                                                  cTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCtggt                                              <  5:205838/107‑1 (MQ=255)
                                                                                                                                   tGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTgg                                                                  >  1:132257/1‑86 (MQ=255)
                                                                                                                                   tGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTgg                                                                  <  2:132257/86‑1 (MQ=255)
                                                                                                                                   tGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTAATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGc   >  7:315885/1‑149 (MQ=255)
                                                                                                                                     cGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGcc  <  2:317078/148‑1 (MQ=255)
                                                                                                                                                ccAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATgg                         <  3:82913/114‑1 (MQ=255)
                                                                                                                                                ccAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATgg                         >  4:82913/1‑114 (MQ=255)
                                                                                                                                                ccAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTc         >  8:154169/1‑130 (MQ=255)
                                                                                                                                                ccAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTc         <  7:154169/130‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                          
GAAGACCCTGGGCGAAGACCCTGCGGCGGTCAACGTGTCGTCGCGCTACTCGGCGCACTATGGCCCCAACCGTCAGGTACAGGGGATCAACAATATCGAGTTGATCACCGACCGCTCTCTGGAAATCAACCTGCGTGAAATCACCCAGGTGAAGAAGGACTGGCCCGACCGCGCACTGATCGTGTCGCTGATGGTGCCGTGCATAGAGGAATCGTGGAAGTTCATCCTGCCGCTGGTGGAAGCCACGGGGGCCGATGGCATCGAGCTGAACTTCGGTTGCC  >  NC_002947/4550763‑4551043

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: