Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,242,090 T→C 5.1% T379A (ACC→GCC)  hmgB ← fumarylacetoacetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,242,0900TC 5.1% 235.3 / 2.2 79T379A (ACC→GCC) hmgBfumarylacetoacetase
Reads supporting (aligned to +/- strand):  ref base T (40/35);  new base C (2/2);  total (42/37)
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.90e-01

TGCGATGACCGTGCCGCGGCACTCACCGAAACCGATGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCGCCACGGTCCAGTACATGCTGCGGGTGTTGCTC  >  NC_002947/5241941‑5242238
                                                                                                                                                     |                                                                                                                                                     
tGCGCTGACAGTGACGCGGCACCCACCGCAACCGATGGCGGCCACGCCGTCACGTGTGCCGCGGGCACGAAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGt                                                                                                                                                       <  1:400822/150‑1 (MQ=255)
 gCGATGACCGTGCCGCGGCACTCACCGAAACCGATGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTg                                                                                                                                                      <  1:633199/150‑1 (MQ=255)
  cGATGACCGTGCCGCGGCACTCACCGAAACCGATGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGa                                                                                                                                                     >  1:782036/1‑150 (MQ=255)
      gACCGTGCCGCGGCACTCACCGAAACCGATGCGGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCCTCCAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGAtctc                                                                                                                                                 <  2:117406/150‑1 (MQ=255)
            gCCGCGGCACTCACCGAAACCGATGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCcagcag                                                                                                                                           <  1:488006/150‑1 (MQ=255)
            gCCGCGGCACTCACCGAAACCGATGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCcagcag                                                                                                                                           <  2:1197991/150‑1 (MQ=255)
                 ggCACTCACCGAACCCGATGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCCGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTAc                                                                                                                                      <  1:1798033/150‑1 (MQ=255)
                     cTCACCGAAACTGATGCTGGTCACGCCGTCACGTGTGCAGAGGGCACGCCGGATGATTTCGTCGCCGTCTTCGACGAAGTTGTGCACCTCGCCGCGCGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGaa                                                                                                                                  <  2:527724/150‑1 (MQ=255)
                      tCACCGAAACCGATGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCTCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAAc                                                                                                                                 >  2:633913/1‑150 (MQ=255)
                      tCACCGAAACCGATGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAAc                                                                                                                                 <  1:340354/150‑1 (MQ=255)
                      tCACCGAAACCGATGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAAc                                                                                                                                 >  1:1314852/1‑150 (MQ=255)
                      tCACCGAAACCGATGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAAc                                                                                                                                 >  2:1740403/1‑150 (MQ=255)
                        aCCGAAACCGATGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGCTTTCGTCGCCGTCTTCGAGGCATTTGCGCACCTCGCCGCTGGCCAGCTCCACAGGCAGTTTGCCCCCTTCGGTGACCTCCAGCAGGCTACAGACCGa                                                                                                                               >  2:957620/1‑150 (MQ=255)
                        aCCGAAACCGATGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGa                                                                                                                               <  1:1680304/150‑1 (MQ=255)
                        aCCGAAACCGATGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGa                                                                                                                               <  1:277917/150‑1 (MQ=255)
                        aCCGAAACCGATGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGa                                                                                                                               >  1:581053/1‑150 (MQ=255)
                            aaaCCGATGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATccc                                                                                                                           <  2:1517039/150‑1 (MQ=255)
                             aaCCGATGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCg                                                                                                                          <  2:1808819/150‑1 (MQ=255)
                                cGATGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGCGATCTCCAGCAGGCTACCGAACGATCCCGGCg                                                                                                                       >  1:267370/1‑150 (MQ=255)
                                  aTGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGCGAtc                                                                                                                                                   <  2:437745/120‑1 (MQ=255)
                                  aTGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGCGAtc                                                                                                                                                   >  1:437745/1‑120 (MQ=255)
                                    ggtGGCCACGCCGTCACGTGTTCAGCGGGTACGCAGGATGATTTCGTCGCCGTCTTCGATGAATTTGCGCACCTCGCCGTGGGCCAGCTCGGCTGGCCGTTTGCCGCCTGCGGTGAGCTCCAGCAGGCTACCGAATGATCCCGGCGTGgc                                                                                                                   <  2:970732/148‑1 (MQ=255)
                                    gCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGgc                                                                                                                   >  2:15180/1‑150 (MQ=255)
                                     cTGGCCACTCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGgcg                                                                                                                  >  2:484773/1‑150 (MQ=255)
                                              ccGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGc                                                                                                         <  1:1058778/150‑1 (MQ=255)
                                                    cGTGTGCAGCGGGCACGTAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGccc                                                                                                   <  1:316206/150‑1 (MQ=255)
                                                           aGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGcccgacccga                                                                                            >  1:980514/1‑150 (MQ=255)
                                                            gCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGaa                                                                                           >  1:1300892/1‑150 (MQ=255)
                                                                        gatgatTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGt                                                                                                       <  1:285018/126‑1 (MQ=255)
                                                                        gatgatTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGt                                                                                                       >  2:285018/1‑126 (MQ=255)
                                                                          tgatTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGt                                                                                                                      <  1:1586105/109‑1 (MQ=255)
                                                                          tgatTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGt                                                                                                                      >  2:1586105/1‑109 (MQ=255)
                                                                          tgatTGCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCt                                                                             <  1:260912/150‑1 (MQ=255)
                                                                                     ccGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAg                                                                 >  2:255456/1‑150 (MQ=255)
                                                                                          ttCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGtt                                                            >  2:915113/1‑150 (MQ=255)
                                                                                          ttCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGtt                                                            <  1:1328421/150‑1 (MQ=255)
                                                                                          ttCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGtt                                                            >  2:1328421/1‑150 (MQ=255)
                                                                                             gAGGAATTTGCGCACCTCGCCGCTGGGCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATATCCAGCAGGCCACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCAAACAGGTCACCTCGCTGCAGCTGACAGACGTTTac                                                         >  2:341132/1‑150 (MQ=255)
                                                                                             gAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGac                                                         >  1:801875/1‑150 (MQ=255)
                                                                                              aGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGaca                                                        >  1:852550/1‑150 (MQ=255)
                                                                                                   tttGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTAGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACtgtg                                                   <  2:1314852/150‑1 (MQ=255)
                                                                                                        gcACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTgg                                                                               >  1:1746233/1‑118 (MQ=255)
                                                                                                        gcACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTgg                                                                               <  2:1746233/118‑1 (MQ=255)
                                                                                                          aCCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGcc                                            >  2:1925538/1‑150 (MQ=255)
                                                                                                            cTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACa                                                                                         <  2:1587193/104‑1 (MQ=255)
                                                                                                            cTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACa                                                                                         >  1:1587193/1‑104 (MQ=255)
                                                                                                            cTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGccac                                          >  1:324172/1‑150 (MQ=255)
                                                                                                            cTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCCCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGccac                                          >  2:719034/1‑150 (MQ=255)
                                                                                                             tCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGccacc                                         >  1:393401/1‑150 (MQ=255)
                                                                                                               gccgcTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGCGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAg                                       <  2:267370/150‑1 (MQ=255)
                                                                                                                  gcTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTg                                                                     <  1:1697265/117‑1 (MQ=255)
                                                                                                                  gcTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTg                                                                     >  2:1697265/1‑117 (MQ=255)
                                                                                                                     ggCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTgcgc                                 <  1:484773/150‑1 (MQ=255)
                                                                                                                        cAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCGCCAc                              >  2:372729/1‑150 (MQ=255)
                                                                                                                        cAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCGCCAc                              >  2:372421/1‑150 (MQ=255)
                                                                                                                                 cGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCCAGCTGCAGccaggt                                                           >  2:1430846/1‑108 (MQ=255)
                                                                                                                                 cGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCCAGCTGCAGccaggt                                                           <  1:1430846/112‑5 (MQ=255)
                                                                                                                                 cGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGc                                                                   >  1:279589/1‑104 (MQ=255)
                                                                                                                                 cGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGc                                                                   <  2:279589/104‑1 (MQ=255)
                                                                                                                                 cGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCGCCACGGTCCAGTa                     <  2:801875/150‑1 (MQ=255)
                                                                                                                                  ggCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCGCCACGGTCCAGTAc                    >  2:1368979/1‑150 (MQ=255)
                                                                                                                                   gCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCTCCACGGTCCAGTACa                   >  1:697476/1‑150 (MQ=255)
                                                                                                                                   gCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCGCCACGGTCCAGTACa                   >  2:1552661/1‑150 (MQ=255)
                                                                                                                                          gccgccTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCGCCACGGTCCAGTACATGCTGCg            >  1:374617/1‑150 (MQ=255)
                                                                                                                                          gccgccTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCGCCACGGTCCAGTACATGCTGCg            >  1:1025349/1‑150 (MQ=255)
                                                                                                                                           ccgccTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCGCCACGGTCCAGTACATGCTGCgg           <  1:633913/150‑1 (MQ=255)
                                                                                                                                           ccgccTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCGCCACGGTCCAGTACATGCTGCgg           <  1:915113/150‑1 (MQ=255)
                                                                                                                                           ccgccTTCGGTGATCTCCAGCAGGCTACCGAACGAGCCCGGCGTGGCGCCCGACATCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCGCCAAGGTCCAGTACATGCTGCgg           <  1:957620/150‑1 (MQ=255)
                                                                                                                                            cgccTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGc                                                                   >  1:1477896/1‑93 (MQ=255)
                                                                                                                                            cgccTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGc                                                                   <  2:1477896/93‑1 (MQ=255)
                                                                                                                                             gccTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCGCCACGGTCCAGTACATGCTGCGGgt         <  1:1740403/150‑1 (MQ=255)
                                                                                                                                              ccTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCATGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCCGCTGCGCCACGGTCCAGTACATGCTGCggacg        >  1:334118/1‑147 (MQ=255)
                                                                                                                                               cTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCGCCACGGTCCAGTACATGCTGCGGgtgt       <  1:15180/150‑1 (MQ=255)
                                                                                                                                               cTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCGCCACGGTCCAGTACATGCTGCGGgtgt       <  1:1552661/150‑1 (MQ=255)
                                                                                                                                                ttCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTgc                                   >  1:1837259/1‑121 (MQ=255)
                                                                                                                                                ttCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTgc                                   <  2:1837259/121‑1 (MQ=255)
                                                                                                                                                ttCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCGCCACGGTCCAGTACATGCTGCGGGTGtt      >  1:664962/1‑150 (MQ=255)
                                                                                                                                                ttCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCGCCACGGTCCAGTACATGCTGCGGGTGtt      <  2:852550/150‑1 (MQ=255)
                                                                                                                                                    gTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCGCCACGGTCCAGTGCATGCTGCGGGTGTTGCTc  >  2:1277638/1‑150 (MQ=255)
                                                                                                                                                    gTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCGCCACGGTCCAGTACATGCTGCGGGTGTTGCTc  >  1:1484353/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCGATGACCGTGCCGCGGCACTCACCGAAACCGATGCTGGCCACGCCGTCACGTGTGCAGCGGGCACGCAGGATGATTTCGTCGCCGTCTTCGAGGAATTTGCGCACCTCGCCGCTGGCCAGCTCGACCGGCAGTTTGCCGCCTTCGGTGATCTCCAGCAGGCTACCGAACGATCCCGGCGTGGCGCCCGACAGCGTGCCCGACCCGAACAGGTCACCTGGCTGCAGCTGGCAGCCGTTGACACTGTGGTGCGCCACCAGCTGCGCCACGGTCCAGTACATGCTGCGGGTGTTGCTC  >  NC_002947/5241941‑5242238

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: