Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,429,125 A→G 8.6% E64G (GAA→GGA)  fklB → FKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,429,1250AG 8.6% 319.7 / 17.7 117E64G (GAA→GGA) fklBFKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base A (56/48);  new base G (3/7);  total (62/55)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

Only 100 of 119 total aligned reads displayed.

TACGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTG  >  NC_000913/4428977‑4429271
                                                                                                                                                     |                                                                                                                                                   
tACGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGaa                                                                                                                                                    >  2:365652/1‑150 (MQ=255)
 aCGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGaaa                                                                                                                                                   <  1:1657971/150‑1 (MQ=255)
  cGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAAt                                                                                                                                                  >  2:1241451/1‑150 (MQ=255)
         ggTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGTTGCATC‑GCGCGCTGCGTGAAATCCACG‑Ag                                                                                                                                          >  1:356110/1‑150 (MQ=255)
         ggTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCATG‑Ag                                                                                                                                          <  2:1158986/150‑1 (MQ=255)
         ggTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑Ag                                                                                                                                          <  2:1070994/150‑1 (MQ=255)
          gTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGATGCCAGAAGCAAAGGTTGCAGGTATTGCAGATGAGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGGTGCGTGAAATCCACG‑Agc                                                                                                                                         >  2:303390/1‑150 (MQ=255)
           tttGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑Agcg                                                                                                                                        <  1:397886/150‑1 (MQ=255)
           tttCCAGGGAGGACAACAACTGAGTGAATCTGGCCTTGAAGGGCTGCTTCCAGAAGCACTGGTAGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATCTGCTTCATC‑GCGCGCTGCGTGAAATCCACG‑Agcg                                                                                                                                        <  1:296461/150‑1 (MQ=255)
           gttGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑Agcg                                                                                                                                        <  1:399422/149‑1 (MQ=255)
              gCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑Agcgcgc                                                                                                                                     >  2:1079787/1‑150 (MQ=255)
                 ggTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGa                                                                                                                                  >  1:462141/1‑150 (MQ=255)
                  gTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGGCGCGCT‑CGTGGATTTCACGCAgcgcgcgat                                                                                                                                  <  1:583766/150‑4 (MQ=38)
                  gTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGGCGCGCT‑CGTGGATTTCACGCAgcgcgcgat                                                                                                                                  <  2:319180/150‑4 (MQ=38)
                      ggCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCg                                                                                                                             >  1:530391/1‑150 (MQ=255)
                       gCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGAt                                                                                                                                 >  2:1692119/1‑145 (MQ=255)
                       gCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGAt                                                                                                                                 <  1:1692119/145‑1 (MQ=255)
                       gCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGt                                                                                                                            >  1:1364270/1‑150 (MQ=255)
                        caacaaCTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGtt                                                                                                                           >  1:768040/1‑150 (MQ=255)
                           caacTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑CGCGCGCCGATGCCGTtcgt                                                                                                                        >  1:1705443/1‑150 (MQ=255)
                             acTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTtcgtcg                                                                                                                      >  2:982710/1‑150 (MQ=255)
                                   tGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCg                                                                                                                <  1:399464/150‑1 (MQ=255)
                                        cTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTcc                                                                                                           >  1:505154/1‑150 (MQ=255)
                                        cTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTcc                                                                                                           <  2:792287/150‑1 (MQ=255)
                                         tGGCCTGGAAGGGCTGCTGCCAGACGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCCAACATCCGGCTGTTCCGGTTGATGTGGTGCCTC‑GCGCGCTGCGGGCAATCCACG‑AGCGCGCCGCTGCCGTTCGGCGTCCGCGTTTCCa                                                                                                          >  1:145551/1‑150 (MQ=255)
                                           gCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAgg                                                                                                        <  1:391987/150‑1 (MQ=255)
                                            ccTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCATCGTCAGTGTTTCCAGGc                                                                                                       >  1:631631/1‑150 (MQ=255)
                                            ccTGCAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGCTGCCGTTCGTCGTCAGCGTTTCCAGGc                                                                                                       >  1:857541/1‑150 (MQ=255)
                                                gacgggCGGCTGGCAGAAGCTCTGGTTGCAGGTTTTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATg                                                                                                   <  2:130843/147‑1 (MQ=255)
                                                  cgggCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgc                                                                                                 <  2:1560993/149‑1 (MQ=255)
                                                  aGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgc                                                                                                 >  2:439076/1‑150 (MQ=255)
                                                    ggCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgctg                                                                                               <  2:462141/150‑1 (MQ=255)
                                                    ggCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgctg                                                                                               >  2:1089415/1‑150 (MQ=255)
                                                      ctgctgCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgctgct                                                                                             <  1:739944/150‑1 (MQ=255)
                                                       tgctgcCAGAAGCACTGGTTGCGGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGctgctg                                                                                            <  2:642872/150‑1 (MQ=255)
                                                        gctgcCAGTAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGGTCCGGTGGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGa                                                                                           <  1:1389852/150‑1 (MQ=255)
                                                        gctgcCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGa                                                                                           >  1:1658351/1‑150 (MQ=255)
                                                           gcCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAgg                                                                                        >  1:1280331/1‑150 (MQ=255)
                                                           gcCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAgg                                                                                        >  2:752380/1‑150 (MQ=255)
                                                             cAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGgtg                                                                                      <  1:281625/150‑1 (MQ=255)
                                                             cAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGgtg                                                                                      <  1:1241451/150‑1 (MQ=255)
                                                             cAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGgtg                                                                                      >  2:281625/1‑150 (MQ=255)
                                                               gAAGCACTGGGTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGtgtg                                                                                    <  1:900180/150‑1 (MQ=255)
                                                                   cACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAAt                                                                                >  1:104040/1‑150 (MQ=255)
                                                                    aCTGGTTGCATGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGACGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATa                                                                               <  2:284101/150‑1 (MQ=255)
                                                                     cTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGAGAAATAc                                                                              >  1:63210/1‑150 (MQ=255)
                                                                        gTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTg                                                                           <  2:768040/150‑1 (MQ=255)
                                                                        gTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTg                                                                           <  2:1439278/150‑1 (MQ=255)
                                                                          tGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTTGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGga                                                                         <  1:1320957/150‑1 (MQ=255)
                                                                          tGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGga                                                                         >  2:1432243/1‑150 (MQ=255)
                                                                          tGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGga                                                                         >  2:1679108/1‑150 (MQ=255)
                                                                           gCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGgaa                                                                        <  1:1311683/150‑1 (MQ=255)
                                                                           gCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGgaa                                                                        >  1:1358360/1‑150 (MQ=255)
                                                                            cAGGGATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGgaag                                                                       <  1:929610/150‑1 (MQ=255)
                                                                                tATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGaaaa                                                                   >  1:780283/1‑150 (MQ=255)
                                                                                tATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGGGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGaaaa                                                                   <  2:881370/150‑1 (MQ=255)
                                                                                 aTTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATAACTGGAAGAAAAc                                                                  >  2:750149/1‑150 (MQ=255)
                                                                                  ttGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACg                                                                 <  2:1650514/150‑1 (MQ=255)
                                                                                     ccGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑CGCGCGCCGCTGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCa                                                              >  1:952994/1‑150 (MQ=255)
                                                                                     ccGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCa                                                              >  1:179375/1‑150 (MQ=255)
                                                                                       gATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCCGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCaaa                                                            >  1:1096516/1‑150 (MQ=255)
                                                                                           cgcTGGATGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaag                                                        >  2:1082300/1‑150 (MQ=255)
                                                                                            gcTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaaga                                                       >  1:795422/1‑150 (MQ=255)
                                                                                             cTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaagaa                                                      <  1:1079787/150‑1 (MQ=255)
                                                                                              tGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaagaag                                                     <  2:1406003/150‑1 (MQ=255)
                                                                                               ggAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAgg                                                    >  1:527845/1‑150 (MQ=255)
                                                                                               ggAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTCCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAgg                                                    <  1:911636/150‑1 (MQ=255)
                                                                                               ggAAGGCAAACATCCGGCTGTTCCGGGTGATGTGGTGCATC‑GCGCGCTGCGTGCAATCCACG‑CGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGCAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAgg                                                    >  1:1483878/1‑150 (MQ=255)
                                                                                                   ggCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGtgtg                                                <  1:732007/150‑1 (MQ=255)
                                                                                                     cAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCTATGGCTGCTGAAGGTGTGAAATACCTgg                                                                          <  1:553971/122‑1 (MQ=255)
                                                                                                     cAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTgg                                                                          >  2:553971/1‑122 (MQ=255)
                                                                                                     cAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCACCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGaa                                              >  1:1206646/1‑150 (MQ=255)
                                                                                                       aaCATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATa                                            >  1:1264924/1‑150 (MQ=255)
                                                                                                        aCATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAg                                           <  1:752380/150‑1 (MQ=255)
                                                                                                         cATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCTAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGc                                          >  2:815227/1‑150 (MQ=255)
                                                                                                         cATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGc                                          <  2:63210/150‑1 (MQ=255)
                                                                                                             cGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCg                                      >  2:452715/1‑150 (MQ=255)
                                                                                                                     ccGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGc                              <  1:1467361/150‑1 (MQ=255)
                                                                                                                     ccGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGc                              >  1:352701/1‑150 (MQ=255)
                                                                                                                         ttGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGc                          >  2:328260/1‑150 (MQ=255)
                                                                                                                         gtGATGTGGTGCATC‑GCGCGCTGCGTGAATTCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGc                          <  1:303390/149‑1 (MQ=255)
                                                                                                                            aTGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAACCGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAAt                       >  1:1282602/1‑150 (MQ=255)
                                                                                                                             tgtgGTGCATC‑GCGCGCTGCGTGCAATCCACG‑AGCGCGCGGATGCCGTTCGTCGTCAGCGTGTCCAGGGGATGGCTGCTGAAGGTGTGATATACCTGGAAGATAACGCCAAAAAAGAAGGTGTGTATAGCACCGAATCTGGGCTGCACtt                      >  1:223792/1‑150 (MQ=255)
                                                                                                                               tggtgCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTcc                    >  1:138041/1‑150 (MQ=255)
                                                                                                                                ggtgCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCcg                   >  1:407383/1‑150 (MQ=255)
                                                                                                                                  tgCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCcgcg                 >  1:1592321/1‑150 (MQ=255)
                                                                                                                                   gCATCGGCGCGCT‑CGTGGATTTCACGCAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGtgtg                                                <  2:899392/119‑1 (MQ=37)
                                                                                                                                   gCATCGGCGCGCT‑CGTGGATTTCACGCAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAAGACCTGGAAGAAAACGCCAAAAAAGAAGGtgtg                                                >  1:905497/1‑119 (MQ=37)
                                                                                                                                   gCATCGGCGCGCT‑CGTGGATTTCACGCAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAACACCTGGAAGAAAACGCCAAAAAAGAAGGtgtg                                                <  2:905497/119‑1 (MQ=37)
                                                                                                                                       c‑gcgcgCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATc            <  2:505154/150‑1 (MQ=255)
                                                                                                                                         gcgcgcTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCa           <  2:179375/150‑1 (MQ=255)
                                                                                                                                         gcgcgcT‑CGTGGATTTCACGCAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCa           >  2:69410/1‑150 (MQ=255)
                                                                                                                                          cgcgcTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGATGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCaa          >  2:929339/1‑150 (MQ=255)
                                                                                                                                               tGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAgg     >  2:372921/1‑150 (MQ=255)
                                                                                                                                               ggCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAgg     <  1:676102/149‑1 (MQ=255)
                                                                                                                                                gCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAggg    <  2:626034/150‑1 (MQ=255)
                                                                                                                                                gCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAggg    <  1:1089415/150‑1 (MQ=255)
                                                                                                                                                gCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAggg    <  1:365652/150‑1 (MQ=255)
                                                                                                                                                 cGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGt   >  2:1543129/1‑150 (MQ=255)
                                                                                                                                                  gTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTg  >  1:1063269/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TACGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTG  >  NC_000913/4428977‑4429271

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