Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,909,174 T→C 8.5% Y40C (TAC→TGC)  PP_1709 ← fumarylacetoacetate hydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,909,1740TC 8.5% 157.6 / 5.2 59Y40C (TAC→TGC) PP_1709fumarylacetoacetate hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (26/28);  new base C (3/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

GGCGCCGCCTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATCCTCCTTACTGGAGTGGTCAGCG  >  NC_002947/1909029‑1909323
                                                                                                                                                  |                                                                                                                                                      
ggCGCCGCCTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTACGGCGCGGCGCCGCAGGa                                                                                                                                                    <  1:859182/150‑1 (MQ=255)
     cgccTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCg                                                                                                                                               <  1:800982/150‑1 (MQ=255)
       ccTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGc                                                                                                                                             >  1:655141/1‑150 (MQ=255)
         tttGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTACGGCGCGGCGCCGCAGGa                                                                                                                                                    <  2:48792/141‑1 (MQ=255)
         tttGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTACGGCGCGGCGCCGCAGGa                                                                                                                                                    >  1:48792/1‑141 (MQ=255)
             gCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACggg                                                                                                                                       <  2:57364/150‑1 (MQ=255)
                 tGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGt                                                                                                                                   >  1:777892/1‑150 (MQ=255)
                     cgcgACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCg                                                                                                                               <  2:655141/150‑1 (MQ=255)
                     cgcgACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCg                                                                                                                               >  2:499183/1‑150 (MQ=255)
                      gcgACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCgg                                                                                                                              >  1:849552/1‑150 (MQ=255)
                            aGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAgg                                                                                                                        <  1:676469/150‑1 (MQ=255)
                              cTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCa                                                                                                                      <  1:1098806/150‑1 (MQ=255)
                               tCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCCCGTTGCTTGCCGGTTTCATGAAAAACC‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGccg                                                                                                                     >  2:202294/1‑148 (MQ=255)
                                          gCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGcc                                                                                                         >  1:1364193/1‑150 (MQ=255)
                                                cTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCgg                                                                                                   <  2:945922/150‑1 (MQ=255)
                                                 tcgtcgGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAg                                                                                                                         <  1:872122/128‑1 (MQ=255)
                                                 tcgtcgGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCggg                                                                                                  >  1:1038015/1‑150 (MQ=255)
                                                 tcgtcgGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTACGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAg                                                                                                                         >  2:872122/1‑128 (MQ=255)
                                                    tcgGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAg                                                                                               <  2:494965/150‑1 (MQ=255)
                                                     cgGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGTATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAgc                                                                                              <  1:175642/150‑1 (MQ=255)
                                                                      ggCAGTTCGCCCAGTGCCTCGACCACTTTGCTTGCCGGTTTCATGAAAAACACCAAGC‑CCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACGCCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCgg                                                                             >  1:886386/1‑150 (MQ=255)
                                                                       gCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGa                                                                            <  2:777892/150‑1 (MQ=255)
                                                                       gCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGCTTAATGAAACCAA‑CGCGCTCCTGCGGCGCCGCGCCGTAGGATTACGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGa                                                                            >  1:1096857/1‑150 (MQ=255)
                                                                              gCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAg                                                                     <  2:29989/150‑1 (MQ=255)
                                                                                ccAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCa                                                                                                                      <  1:1094457/100‑1 (MQ=255)
                                                                                ccAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCa                                                                                                                      >  2:1094457/1‑100 (MQ=255)
                                                                                       cTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGAt                                                            >  2:1271716/1‑150 (MQ=255)
                                                                                       cTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGAt                                                            >  2:347640/1‑150 (MQ=255)
                                                                                        tCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATg                                                           <  2:1096857/150‑1 (MQ=255)
                                                                                        tCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATg                                                           <  2:849552/150‑1 (MQ=255)
                                                                                               cGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCcggcg                                                    <  1:1164827/150‑1 (MQ=255)
                                                                                                 ttgcttgcCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCggg                                                  <  1:1271716/150‑1 (MQ=255)
                                                                                                 ttgcttgcCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCggg                                                  >  1:1215953/1‑150 (MQ=255)
                                                                                                     ttgcCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGGGGCCCATCCGGGATGCCCTCCTGG‑GCAGCTCAGCGGTGCCGCGGATCGGCAGGCTGACGATAGCCGGGGGGGTGa                                              >  1:525498/1‑150 (MQ=255)
                                                                                                     ttgcCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCg                                                                         >  1:1018280/1‑123 (MQ=255)
                                                                                                     ttgcCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCg                                                                         <  2:1018280/123‑1 (MQ=255)
                                                                                                      tgcCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGaa                                             <  1:1200653/150‑1 (MQ=255)
                                                                                                         cGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACAcc                                                                                                              >  2:167251/1‑82 (MQ=255)
                                                                                                         cGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACAcc                                                                                                              <  1:167251/82‑1 (MQ=255)
                                                                                                             ttCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTaa                                      <  2:188869/150‑1 (MQ=255)
                                                                                                                 tGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATg                                                           >  1:1413540/1‑125 (MQ=255)
                                                                                                                 tGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATg                                                           <  2:1413540/125‑1 (MQ=255)
                                                                                                                 tGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGtga                                  >  1:1172274/1‑148 (MQ=255)
                                                                                                                 tGAAAAACA‑CAGGCTCCTGCGGCGCCACGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCa                                  >  1:831764/1‑150 (MQ=255)
                                                                                                                   aaaaaCA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATg                                >  2:232724/1‑150 (MQ=255)
                                                                                                                       aca‑cAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTg                            <  1:499183/150‑1 (MQ=255)
                                                                                                                         gccggaatCCTACGGCGCGGCGCCGCAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTgg                           >  2:1184896/8‑150 (MQ=255)
                                                                                                                              gCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGAtcct                      >  1:297358/1‑150 (MQ=255)
                                                                                                                              aatCCTACGGCGCGGCGCCGCAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGAtcct                      >  1:1293549/3‑150 (MQ=255)
                                                                                                                                 ccTGCGGCGCCGCGCCGTAGGATTCCGGCACCGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATcctcct                   <  2:337599/150‑1 (MQ=255)
                                                                                                                                 ccTACGGCGCGGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATcctcct                   >  1:360228/1‑150 (MQ=255)
                                                                                                                                   tGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACa                                           >  2:278747/1‑124 (MQ=255)
                                                                                                                                   tGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACa                                           <  1:278747/124‑1 (MQ=255)
                                                                                                                                         gccgcgccgTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATCCTCCTTACTGGAg           <  2:297358/150‑1 (MQ=255)
                                                                                                                                         gccgcgccgTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCCCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATCCTCCTTACTGGAg           <  1:455395/150‑1 (MQ=255)
                                                                                                                                                cgTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATCCTCCTTACTGGAGTGGTCAg    <  1:1228106/150‑1 (MQ=255)
                                                                                                                                                 gTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGc                                                         <  1:943622/96‑1 (MQ=255)
                                                                                                                                                 gTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGc                                                         >  2:943622/1‑96 (MQ=255)
                                                                                                                                                 gTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATCCTCCTTACTGGAGTGGTCAGc   >  1:40666/1‑150 (MQ=255)
                                                                                                                                                  tAGGATTCCGGCCACGGGTAGTTGCGGGCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATCCTCCTTACTGGAGTGGTCAGCg  >  2:603645/1‑150 (MQ=255)
                                                                                                                                                  tAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATCCTCCTTACt               >  1:670914/1‑137 (MQ=255)
                                                                                                                                                  tAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATCCTCCTTACt               <  2:670914/137‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GGCGCCGCCTTTGGCGATGGCCGCGACCAGCTCGATTTCATGGCAGAACTCGTCGGTCTGCGGCGGGAAGGGCAGTTCGCCCAGTGCCTCGACCACGTTGCTTGCCGGTTTCATGAAAAACA‑CAGGCTCCTGCGGCGCCGCGCCGTAGGATTCCGGCCACGGGTAGTTGCGGCCGAGGCAGA‑ACACCCGGCCGACCGGGAAGCGCGCATCGCTGCCGCGGATCGGCAGGCTGACGATGGCCGGCGGGTTGAACAGGTAAGTCATGGGTGGATCCTCCTTACTGGAGTGGTCAGCG  >  NC_002947/1909029‑1909323

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