Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,757,759 (G)6→5 coding (1118/1392 nt) selA ← selenocysteine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,757,7540G.99.6% 1198.9 / ‑3.5 255coding (1123/1392 nt)selAselenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base . (119/135);  total (119/136)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01

Only 100 of 294 total aligned reads displayed.

CGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAGGGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTT  >  NC_000913/3757662‑3757843
                                                                                            |                                                                                         
cGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAg                                                                                           <  2:97080/93‑1 (MQ=255)
cGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAg                                                                                           <  2:4008187/93‑1 (MQ=255)
cGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAg                                                                                           <  2:1508143/93‑1 (MQ=255)
 ggCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAgg                                                                                          <  1:4104009/93‑1 (MQ=255)
 ggCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAgg                                                                                          <  1:1633418/93‑1 (MQ=255)
  gCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAggg                                                                                         <  1:5248918/93‑1 (MQ=255)
   cAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAgggg                                                                                        <  1:6462478/93‑1 (MQ=255)
    aGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAggggg                                                                                       >  2:2450018/1‑93 (MQ=255)
     ggCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTATCGCAAACTCCGCGCCGTAATGTGCGGCACgggggc                                                                                      >  1:576146/1‑92 (MQ=255)
     ggCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAgggggc                                                                                      <  1:5012213/93‑2 (MQ=255)
     ggCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAgggggc                                                                                      >  2:3843869/1‑92 (MQ=255)
      gCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAgggggcc                                                                                     <  1:1073674/93‑3 (MQ=255)
       cGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAgggggcct                                                                                    <  2:2373760/93‑4 (MQ=255)
       cGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAgggggcct                                                                                    <  2:3655640/93‑4 (MQ=255)
         aTCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGt                                                                                 >  2:2787633/1‑93 (MQ=255)
          tCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTa                                                                                >  2:4281093/1‑93 (MQ=255)
           cAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTaa                                                                               >  1:2320590/1‑93 (MQ=255)
            aaCCGGCAGCGAACCATTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTaaa                                                                              <  1:709320/93‑1 (MQ=255)
            aaCCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTaaa                                                                              <  2:3874113/93‑1 (MQ=255)
            aaCCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTaaa                                                                              <  2:2576875/93‑1 (MQ=255)
             aCCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAAc                                                                             >  2:1990674/1‑93 (MQ=255)
                 gCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTtg                                                                         <  1:2033775/93‑1 (MQ=255)
                  cAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTtgt                                                                        <  1:3046517/93‑1 (MQ=255)
                    gCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGc                                                                      <  2:5868836/93‑1 (MQ=255)
                     cGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGcc                                                                     <  1:5401521/93‑1 (MQ=255)
                      gAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAGGGGGGCCTGTAAACGTTGTGcc                                                                     <  1:566953/93‑1 (MQ=255)
                      gAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCt                                                                    <  2:4524967/93‑1 (MQ=255)
                      gAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCt                                                                    >  1:3164133/1‑93 (MQ=255)
                       aaCCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTg                                                                   >  1:1096383/1‑93 (MQ=255)
                       aaCCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTg                                                                   <  1:1906756/93‑1 (MQ=255)
                       aaCCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTg                                                                   <  1:471622/93‑1 (MQ=255)
                       aaCCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTg                                                                   <  2:622198/93‑1 (MQ=255)
                       aaCCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTg                                                                   >  2:3004954/1‑93 (MQ=255)
                       aaCCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTg                                                                   >  2:5482824/1‑93 (MQ=255)
                       aaCCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTg                                                                   >  2:4422571/1‑93 (MQ=255)
                       aaCCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGACGTAACGCGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTg                                                                   >  1:1319055/1‑93 (MQ=255)
                         ccACTGTCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGa                                                                 <  2:5984092/93‑1 (MQ=255)
                         ccACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGa                                                                 <  2:3452196/93‑1 (MQ=255)
                         ccACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGa                                                                 <  1:3504940/93‑1 (MQ=255)
                           aCTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGAtt                                                               >  1:42669/1‑93 (MQ=255)
                           aCTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGAtt                                                               >  2:5078462/1‑93 (MQ=255)
                                cAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAAt                                                          <  1:3695603/93‑1 (MQ=255)
                                 aaTCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATg                                                         <  1:2710834/93‑1 (MQ=255)
                                 aaTCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATg                                                         >  2:5555958/1‑93 (MQ=255)
                                 aaTCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATg                                                         <  1:2179267/93‑1 (MQ=255)
                                  aTCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGa                                                        <  2:1701995/93‑1 (MQ=255)
                                    cTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGAcc                                                      >  1:2476131/1‑93 (MQ=255)
                                     tGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACct                                                     <  1:4637749/93‑1 (MQ=255)
                                         aaaGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTctgc                                                 >  1:3628709/1‑93 (MQ=255)
                                          aaGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTctgcg                                                >  1:787632/1‑93 (MQ=255)
                                          aaGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTctgcg                                                >  1:3759555/1‑93 (MQ=255)
                                          aaGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTctgcg                                                >  1:5490308/1‑93 (MQ=255)
                                          aaGACATGGCATAACCTGTACCGCAAACTACGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTctgcg                                                >  1:4807722/1‑93 (MQ=255)
                                            gACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTctgcgct                                              >  1:1465416/1‑93 (MQ=255)
                                            gACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTctgcgct                                              <  1:4239522/93‑1 (MQ=255)
                                            gACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTctgcgct                                              >  2:1393469/1‑93 (MQ=255)
                                             aCATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTctgcgctg                                             >  2:2743141/1‑93 (MQ=255)
                                               aTGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTctgcgctgcg                                           <  2:1890655/93‑1 (MQ=255)
                                                 ggCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGTGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCggg                                         <  1:3570209/93‑1 (MQ=255)
                                                 ggCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCggg                                         >  1:665932/1‑93 (MQ=255)
                                                 ggCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCggg                                         <  2:100416/93‑1 (MQ=255)
                                                  gCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGt                                        <  1:93193/93‑1 (MQ=255)
                                                   cATAACCTGTACGGCAACGTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGAATTTGAATGACCTCTGCGCTGCGGGTa                                       <  1:3320273/93‑1 (MQ=255)
                                                   cATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTa                                       <  2:3188886/93‑1 (MQ=255)
                                                    aTAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCGTCGGGTaa                                      >  2:352681/1‑93 (MQ=255)
                                                     tAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAg                                     <  2:2292446/93‑1 (MQ=255)
                                                     tAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAg                                     <  1:4296808/93‑1 (MQ=255)
                                                     tAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGATGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGGAAg                                     >  1:6346582/1‑93 (MQ=255)
                                                      aaCCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGc                                    <  2:5753247/93‑1 (MQ=255)
                                                      aaCCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGc                                    >  2:5326995/1‑93 (MQ=255)
                                                       aCCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCa                                   >  2:3778934/1‑93 (MQ=255)
                                                        ccTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAg                                  <  1:1510254/93‑1 (MQ=255)
                                                        ccTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAg                                  >  1:3439609/1‑93 (MQ=255)
                                                          tGTACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGgc                                <  2:6148690/93‑1 (MQ=255)
                                                            tACCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGgcgc                              <  1:1434221/93‑1 (MQ=255)
                                                             aCCGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCa                             <  2:700135/93‑1 (MQ=255)
                                                              ccGCAACCTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGT‑AACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAgg                           <  1:6217755/93‑1 (MQ=255)
                                                              ccGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAg                            <  2:2688433/93‑1 (MQ=255)
                                                              ccGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAg                            <  2:2601587/93‑1 (MQ=255)
                                                               cGCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAgg                           <  2:2492279/93‑1 (MQ=255)
                                                                gCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAggg                          <  2:4664400/93‑1 (MQ=255)
                                                                gCAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAggg                          <  1:2561320/93‑1 (MQ=255)
                                                                 cAAACTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGt                         <  1:5989327/93‑1 (MQ=255)
                                                                  aaccTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTc                        <  1:5645919/90‑1 (MQ=255)
                                                                  aaaCTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTc                        <  1:4194482/93‑1 (MQ=255)
                                                                   aaCTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGAGCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCg                       >  2:484862/1‑93 (MQ=255)
                                                                     cTCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGt                     >  2:1321047/1‑93 (MQ=255)
                                                                      tCCGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTa                    >  1:2925080/1‑93 (MQ=255)
                                                                       ccGCGCCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTaa                   <  1:5589147/93‑1 (MQ=255)
                                                                          cgcCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAAttt                <  2:3347565/93‑1 (MQ=255)
                                                                          cgcCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAAttt                <  1:4413193/93‑1 (MQ=255)
                                                                           gcCGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAAtttt               >  1:2743418/1‑93 (MQ=255)
                                                                             cGTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTc             >  1:1263818/1‑93 (MQ=255)
                                                                              gTAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCa            >  1:2407244/1‑93 (MQ=255)
                                                                               tAATGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCAc           <  1:3589614/93‑1 (MQ=255)
                                                                                 aTGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTc         <  2:527599/93‑1 (MQ=255)
                                                                                 aTGTGCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTc         >  2:565095/1‑93 (MQ=255)
                                                                                  tgtgCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGCAATTTTTCACTCa        >  1:1803201/1‑93 (MQ=255)
                                                                                   gtgCGGCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCag       <  2:119081/93‑1 (MQ=255)
                                                                                       ggCAA‑GGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCt   <  1:3932764/93‑1 (MQ=255)
                                                                                            |                                                                                         
CGGCAGGCGATCAACCGGCAGCGAACCACTGCCAATCTGCGAAAGACATGGCATAACCTGTACCGCAAACTCCGCGCCGTAATGTGCGGCAAGGGGGGCCTGTAAACGTTGTGCCTGGATTTGAATGACCTCTGCGCTGCGGGTAAGCAGGCGCAGGGTCGGTAATTTTTCACTCAGAGCTT  >  NC_000913/3757662‑3757843

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