Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,695,310	T→C	Q252R (CAG→CGG)	hdhA ←	NAD‑dependent 7alpha‑hydroxysteroid dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,695,310	0	T	C	100.0%	276.3 / NA	91	Q252R (CAG→CGG)	hdhA	NAD‑dependent 7alpha‑hydroxysteroid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (39/52);  total (39/52)

CAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCT  >  NC_000913/1695221‑1695399                                                                                          |                                                                                          cAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTa                                                                                        <  2:2439687/93‑1 (MQ=255) cAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTa                                                                                        >  2:4507432/1‑93 (MQ=255)  aGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTAc                                                                                       <  1:65254/93‑1 (MQ=255)   gTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTAcc                                                                                      >  2:1726854/1‑93 (MQ=255)     tATTTAGCAGAGACGTTGATGAGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTAcccc                                                                                    >  1:2316052/1‑93 (MQ=255)     tATTTAGCAGAGACGTTGATGAGTTTATCCCGGCAAACAACACGCACGGTGGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTAcccc                                                                                    <  1:3297863/93‑1 (MQ=255)      aTTTAGCAGAGACGTTGATGAGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTAccccc                                                                                   >  2:3205721/1‑92 (MQ=255)      aTTTAGCAGAGACGTTGATGAGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCcca                                                                                   >  1:3933680/1‑93 (MQ=255)       tttAGCAGAGACGTTGATGAGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCccac                                                                                  >  1:436758/1‑93 (MQ=255)        ttAGCAGAGACGTTGATGAGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCccacc                                                                                 >  1:2983709/1‑93 (MQ=255)              gagaCGTTGATGAGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGga                                                                           >  2:4140393/1‑93 (MQ=255)              gagaCGTTGATGAGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAAATGAGCTCCCGTACCCCACCACCGga                                                                           >  1:1469552/1‑93 (MQ=255)                gaCGTTGATGAGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGgaga                                                                         >  1:2363445/1‑93 (MQ=255)                   gTTGATGAGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACgg                                                                      >  1:1227936/1‑93 (MQ=255)                   gTTGATGAGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACgg                                                                      <  2:409104/93‑1 (MQ=255)                   gTTGATGAGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACgg                                                                      <  1:2659691/93‑1 (MQ=255)                   gTTGATGAGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACgg                                                                      >  1:2171697/1‑93 (MQ=255)                    ttGATGAGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACggg                                                                     >  2:2589378/1‑92 (MQ=255)                     tgatgaGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTg                                                                    <  1:4213391/93‑1 (MQ=255)                       atgaGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTgag                                                                  >  2:2097689/1‑93 (MQ=255)                       atgaGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGCTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGCGACGGTgag                                                                  >  2:2097688/1‑93 (MQ=255)                        tgaGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTgaga                                                                 >  1:1617587/1‑93 (MQ=255)                         gaGTTTATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGaa                                                                >  1:1544516/1‑93 (MQ=255)                           gTTTATCCCGGCAAACAACACGCACGGGTGTTTACCGGTCCGTTAGTGTATTTATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGCAtt                                                              >  1:1862682/1‑93 (MQ=255)                             ttATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTg                                                            >  2:4091333/1‑93 (MQ=255)                              tATCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGt                                                           <  1:2586955/93‑1 (MQ=255)                                tCCCGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTcc                                                         <  1:3023869/93‑1 (MQ=255)                                   cGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCt                                                      <  2:1469552/93‑1 (MQ=255)                                   cGGCAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCt                                                      <  1:65072/93‑1 (MQ=255)                                    ggcaaaCAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCtt                                                     <  1:2009964/90‑1 (MQ=255)                                      caaaCAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTAc                                                   >  2:4327490/2‑93 (MQ=255)                                       aaaCAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTccc                                                  >  1:3470710/1‑93 (MQ=255)                                         acaacaCGCGCGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCa                                                <  2:2169272/93‑1 (MQ=255)                                          caacaCGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAg                                               >  2:1011859/1‑93 (MQ=255)                                          caacaCGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAg                                               <  2:424450/93‑1 (MQ=255)                                           aacaCGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGc                                              >  1:1824145/1‑93 (MQ=255)                                            acacGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCt                                             >  2:242295/1‑93 (MQ=255)                                            acacGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCt                                             <  1:976239/93‑1 (MQ=255)                                            acacGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCt                                             <  2:24633/93‑1 (MQ=255)                                             cacgcacgGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTc                                            >  2:2544879/1‑93 (MQ=255)                                             cacgcacgGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTc                                            >  2:4016898/1‑93 (MQ=255)                                                gcacgGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCgca                                         <  2:2824363/93‑1 (MQ=255)                                                 cacgGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCgcag                                        <  1:3601180/93‑1 (MQ=255)                                                    ggTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGcagcag                                     >  1:2795506/1‑93 (MQ=255)                                                     gttgttTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAgg                                    >  2:4636102/1‑93 (MQ=255)                                                      ttgttTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGc                                   <  1:659620/93‑1 (MQ=255)                                                       tgttTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCg                                  <  1:3985224/93‑1 (MQ=255)                                                        gttTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGa                                 >  2:1891221/1‑93 (MQ=255)                                                        gttTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGa                                 <  2:3220784/93‑1 (MQ=255)                                                        gttTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGa                                 >  1:3070803/1‑93 (MQ=255)                                                         tttACCGTTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAg                                <  2:126410/93‑1 (MQ=255)                                                         tttACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAg                                <  1:722866/93‑1 (MQ=255)                                                         tttACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAg                                <  2:762228/93‑1 (MQ=255)                                                          ttACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGc                               <  1:2853004/93‑1 (MQ=255)                                                            aCCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCaa                             <  2:1021820/93‑1 (MQ=255)                                                             ccGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCaaa                            <  2:3021991/93‑1 (MQ=255)                                                             ccGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCaaa                            <  1:2929133/93‑1 (MQ=255)                                                              cGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAg                           <  1:2343193/93‑1 (MQ=255)                                                              cGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAg                           <  1:3445941/93‑1 (MQ=255)                                                              cGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGACAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAg                           <  1:3385356/93‑1 (MQ=255)                                                               ggTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAgg                          >  2:4251414/1‑93 (MQ=255)                                                               ggTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAgg                          <  2:1232958/93‑1 (MQ=255)                                                                gTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGa                         <  1:2376555/93‑1 (MQ=255)                                                                gTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGa                         >  2:3320493/1‑93 (MQ=255)                                                                  ccGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAAc                       <  1:152790/93‑1 (MQ=255)                                                                   cGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACa                      <  1:3116992/93‑1 (MQ=255)                                                                   cGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACa                      <  2:3166496/93‑1 (MQ=255)                                                                   cGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACa                      >  1:4054633/1‑93 (MQ=255)                                                                    gTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAg                     <  2:4245152/93‑1 (MQ=255)                                                                    gTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAg                     <  2:498992/93‑1 (MQ=255)                                                                    gTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAg                     <  1:3803216/93‑1 (MQ=255)                                                                    gTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAg                     <  2:1038536/93‑1 (MQ=255)                                                                     ttAGTGTATTAATTGAGCTCCCGTACCCCACCGCCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAgc                    <  2:3057546/93‑1 (MQ=255)                                                                     ttAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAgc                    >  1:1395772/1‑93 (MQ=255)                                                                        gtgtATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCt                 >  2:4125134/1‑93 (MQ=255)                                                                         tgtATTAATGGACCTCCCGCCCCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTg                <  2:187171/93‑1 (MQ=255)                                                                          gtATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGc               >  1:1574102/1‑93 (MQ=255)                                                                           tATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCg              >  1:1751542/1‑93 (MQ=255)                                                                              tAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTa           <  2:3504627/93‑1 (MQ=255)                                                                              tAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTa           <  1:615564/93‑1 (MQ=255)                                                                               aaTTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCNGCGTTAg          <  1:2846547/93‑1 (MQ=255)                                                                                 ttGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCa        <  2:3812641/93‑1 (MQ=255)                                                                                 ttGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCa        <  1:4578429/93‑1 (MQ=255)                                                                                  tGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCaa       <  1:372730/93‑1 (MQ=255)                                                                                   gAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAat      <  2:851090/93‑1 (MQ=255)                                                                                   gAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAat      >  1:2091209/1‑93 (MQ=255)                                                                                   gAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAat      >  2:3817218/1‑93 (MQ=255)                                                                                   gAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAat      <  2:3797799/93‑1 (MQ=255)                                                                                    aGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAata     <  1:4656176/93‑1 (MQ=255)                                                                                     gCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAatat    <  1:4257585/93‑1 (MQ=255)                                                                                     gCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAatat    <  1:4577752/93‑1 (MQ=255)                                                                                      cTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATc   <  2:4225941/93‑1 (MQ=255)                                                                                      cTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATc   <  1:2393330/93‑1 (MQ=255)                                                                                       tCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGTGCTGCGTTAGCAATATCt  <  2:3019371/93‑1 (MQ=255)                                                                                          |                                                                                          CAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCT  >  NC_000913/1695221‑1695399

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