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_0009131,695,3100TC99.5% 629.7 / ‑5.4 197Q252R (CAG→CGG) hdhANAD‑dependent 7alpha‑hydroxysteroid dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (80/116);  total (80/117)
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.90e-01

Only 100 of 202 total aligned reads displayed.

CCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCT  >  NC_000913/1695220‑1695399
                                                                                          |                                                                                         
ccAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGt                                                                                         <  1:3424694/93‑1 (MQ=255)
 cAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTa                                                                                        >  1:613876/1‑93 (MQ=255)
  aGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTAc                                                                                       <  2:1828661/93‑1 (MQ=255)
  aGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTAc                                                                                       >  1:2621953/1‑93 (MQ=255)
   gTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTAcc                                                                                      >  1:3063349/1‑93 (MQ=255)
     tATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTAcccc                                                                                    <  1:5456555/93‑1 (MQ=255)
      aTTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCcca                                                                                   >  2:4725152/1‑93 (MQ=255)
        ttAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCccacc                                                                                 <  1:6461979/93‑1 (MQ=255)
        ttAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCccacc                                                                                 <  2:5387991/93‑1 (MQ=255)
          aGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCcaccac                                                                               >  2:2738734/1‑93 (MQ=255)
             agagACGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCgg                                                                            <  1:6250768/93‑1 (MQ=255)
             agagACGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCgg                                                                            >  2:4244556/1‑93 (MQ=255)
              gagaCGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCCCCggg                                                                           >  2:827303/1‑92 (MQ=255)
              gagaCGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGga                                                                           >  1:1412178/1‑93 (MQ=255)
              gagaCGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGga                                                                           <  1:1769153/93‑1 (MQ=255)
              gagaCGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTACGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGga                                                                           >  2:916970/1‑93 (MQ=255)
                gaCGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGgaga                                                                         >  1:2540904/1‑93 (MQ=255)
                  cGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACg                                                                       <  2:2905079/93‑1 (MQ=255)
                   gTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACgg                                                                      <  1:1875235/93‑1 (MQ=255)
                    ttGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGt                                                                     <  1:6289725/93‑1 (MQ=255)
                      gatgaGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTga                                                                   >  1:2937076/1‑93 (MQ=255)
                      gatgaGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTga                                                                   <  2:5451936/93‑1 (MQ=255)
                        tgaGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTgaga                                                                 >  2:2374086/1‑93 (MQ=255)
                        tgaGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTgaga                                                                 >  1:5177281/1‑93 (MQ=255)
                         gcgTTTATCCCGGAAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGaa                                                                <  2:2051641/91‑1 (MQ=255)
                         gaGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGaa                                                                <  1:382269/93‑1 (MQ=255)
                           gTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAAtt                                                              <  2:5221744/93‑1 (MQ=255)
                              tATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGt                                                           >  1:5880094/1‑93 (MQ=255)
                               aTCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTc                                                          <  1:4028054/93‑1 (MQ=255)
                                tCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTcc                                                         >  2:3643489/1‑93 (MQ=255)
                                tCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTcc                                                         >  2:4084038/1‑93 (MQ=255)
                                tCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTcc                                                         >  1:5608524/1‑93 (MQ=255)
                                tCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTcc                                                         >  2:3456508/1‑93 (MQ=255)
                                  ccGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGc                                                       <  1:1602906/93‑1 (MQ=255)
                                  ccGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGc                                                       <  2:3890448/93‑1 (MQ=255)
                                   cGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCt                                                      <  1:3759340/93‑1 (MQ=255)
                                    ggTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCtt                                                     >  1:5046680/1‑93 (MQ=255)
                                     gTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTa                                                    >  2:5015273/1‑93 (MQ=255)
                                     gTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTa                                                    <  1:2978942/93‑1 (MQ=255)
                                        aacaacACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTAccc                                                 <  1:3193468/93‑1 (MQ=255)
                                        aacaacACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTAccc                                                 <  2:312904/93‑1 (MQ=255)
                                        aacaacACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTAccc                                                 <  2:3046748/93‑1 (MQ=255)
                                         acaacaCGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCa                                                >  1:2936338/1‑93 (MQ=255)
                                         acaacaCGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCa                                                <  1:3835503/93‑1 (MQ=255)
                                           aacaCGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGc                                              <  1:2496697/93‑1 (MQ=255)
                                           aacaCGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGc                                              <  2:103712/93‑1 (MQ=255)
                                           aacaCGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGc                                              <  2:3238020/93‑1 (MQ=255)
                                           aacaCGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGc                                              >  1:4196378/1‑93 (MQ=255)
                                            acacGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCt                                             >  1:1019456/1‑93 (MQ=255)
                                             cacgtacgGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTc                                            <  1:3730762/93‑1 (MQ=255)
                                             cacgcacgGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTc                                            <  2:3729869/93‑1 (MQ=255)
                                              acgcacgGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCg                                           <  1:3362076/93‑1 (MQ=255)
                                                   cgGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCgcagca                                      >  1:4598246/1‑93 (MQ=255)
                                                    ggTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGcagcag                                     <  2:6628442/93‑1 (MQ=255)
                                                      ttgttTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAggg                                   >  2:420032/1‑92 (MQ=255)
                                                      ttgttTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGc                                   >  1:5794881/1‑93 (MQ=255)
                                                        gttTACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGa                                 <  2:5570838/93‑1 (MQ=255)
                                                           tACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCa                              <  1:4964186/93‑1 (MQ=255)
                                                           tACCGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCa                              >  2:5708740/1‑93 (MQ=255)
                                                             ccGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCaaa                            >  2:6246964/1‑93 (MQ=255)
                                                              cGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTGCCCAGCTCGCAGCAGGCGAGCAAAg                           <  2:4397541/93‑1 (MQ=255)
                                                              cGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAg                           <  1:3163471/93‑1 (MQ=255)
                                                              cGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAg                           <  2:2496835/93‑1 (MQ=255)
                                                              cGGTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAg                           <  1:435359/93‑1 (MQ=255)
                                                                gTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGa                         <  1:3650060/93‑1 (MQ=255)
                                                                gTCCGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGa                         <  2:334124/93‑1 (MQ=255)
                                                                   cGTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACa                      <  2:864362/93‑1 (MQ=255)
                                                                    gTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAg                     <  2:4344850/93‑1 (MQ=255)
                                                                    gTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAg                     <  2:5977732/93‑1 (MQ=255)
                                                                    gTTAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAg                     <  1:1179136/93‑1 (MQ=255)
                                                                     ttAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAgc                    <  1:5963834/93‑1 (MQ=255)
                                                                     ttAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAgc                    >  2:5626628/1‑93 (MQ=255)
                                                                      tAGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAgcg                   >  2:1874031/1‑93 (MQ=255)
                                                                       aGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAgcgc                  <  2:5816836/93‑1 (MQ=255)
                                                                       aGTGTATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAgcgc                  <  1:811148/93‑1 (MQ=255)
                                                                          gtATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGc               <  1:1378048/93‑1 (MQ=255)
                                                                          gtATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGc               <  1:3884124/93‑1 (MQ=255)
                                                                           tATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCg              <  2:3500504/93‑1 (MQ=255)
                                                                           tATTAATTGAGCTCCCGTACCCCACCACCGGAGACGGGGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGGGCTGCg              >  1:6507717/1‑93 (MQ=255)
                                                                            aTTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGt             <  2:1953738/93‑1 (MQ=255)
                                                                            aTTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGt             <  2:511210/93‑1 (MQ=255)
                                                                            aTTAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGAGCGCAAAGGAAGAGCGCTGCGt             >  1:1170550/1‑93 (MQ=255)
                                                                            aTTAATGGAGATCCCGTACCCAACCACCGGAGACGGTGAGAATTTGTCCGCTTGCCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGt             <  1:2878352/93‑1 (MQ=255)
                                                                             ttAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGtt            <  1:4857041/93‑1 (MQ=255)
                                                                             ttAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGtt            >  1:1530134/1‑93 (MQ=255)
                                                                             ttAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGtt            >  2:5738444/1‑93 (MQ=255)
                                                                             ttAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGtt            >  1:3943492/1‑93 (MQ=255)
                                                                             ttAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGtt            >  2:1284403/1‑93 (MQ=255)
                                                                             ttAATTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGtt            >  2:4460367/1‑93 (MQ=255)
                                                                               aaTTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAg          <  2:5919492/93‑1 (MQ=255)
                                                                               aaTTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAg          >  1:2096899/1‑93 (MQ=255)
                                                                               aaTTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAg          <  2:3396217/93‑1 (MQ=255)
                                                                                aTTGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGc         >  2:6146289/1‑93 (MQ=255)
                                                                                  tGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCaa       <  2:5503981/93‑1 (MQ=255)
                                                                                  tGAGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCaa       <  1:2238820/93‑1 (MQ=255)
                                                                                    aGCTCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAata     <  2:1696447/93‑1 (MQ=255)
                                                                                       tCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGGGAGCAAAGGAACAGCGCTGCGTTAGCAATATCt  >  2:2858751/1‑93 (MQ=255)
                                                                                       tCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCt  <  1:4665000/93‑1 (MQ=255)
                                                                                       tCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCt  <  1:5264183/93‑1 (MQ=255)
                                                                                       tCCCGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCt  >  2:5083251/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCAGTTATTTAGCAGAGACGTTGATGAGTTTATCCCGGTAAACAACACGCACGGTTGTTTACCGGTCCGTTAGTGTATTAATTGAGCTCCTGTACCCCACCACCGGAGACGGTGAGAATTTGTCCGCTTACCCAGCTCGCAGCAGGCGAGCAAAGGAACAGCGCTGCGTTAGCAATATCT  >  NC_000913/1695220‑1695399

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