Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,429,120 Δ1 bp 6.0% coding (186/621 nt) 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,1200G. 6.0% 360.9 / 17.3 116coding (186/621 nt)fklBFKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base G (57/51);  new base . (2/5);  total (59/57)
Fisher's exact test for biased strand distribution p-value = 2.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

Only 100 of 117 total aligned reads displayed.

GCTACGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGG  >  NC_000913/4428975‑4429268
                                                                                                                                                  |                                                                                                                                                     
gcTACGGCATTGGTTTGCAGGTAGGGCAACAACTGCGGGAATCTGGCCTGGAAGGTCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTg                                                                                                                                                   <  1:1330916/150‑1 (MQ=255)
  tACGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGaa                                                                                                                                                 >  2:496225/1‑150 (MQ=255)
    tgccATTGGTTTGCAGGTAGGGCAACAAGTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCATAAGCACAGGTTGCAGGTATTGCCGATGCGCTGGAAGGCTAACATCCGGCTGTTCCGGATGATGTGGTGCAGC‑GCGCGCTGCGTGAAAt                                                                                                                                               <  1:825544/147‑1 (MQ=255)
     ggCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATc                                                                                                                                              >  1:872029/1‑150 (MQ=255)
           ggTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑Ag                                                                                                                                       >  1:939050/1‑150 (MQ=255)
           ggTGTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCGGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑Ag                                                                                                                                       <  1:976236/150‑1 (MQ=255)
                  aGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑Agcgcgcga                                                                                                                                <  1:735538/150‑3 (MQ=255)
                  aGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGGTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGGCGCGCT‑CGTGGATTTCACGCAgcgcgcg                                                                                                                                 <  2:907940/150‑2 (MQ=255)
                    gTAGGGCAACAACTGTGTGAATCTGGCCTGTACGGGCTGCTGCCAGAATCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGAt                                                                                                                              <  1:1131021/150‑1 (MQ=255)
                    gTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGAt                                                                                                                              <  1:625323/150‑1 (MQ=255)
                     gaGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCTGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATg                                                                                                                             <  2:147649/149‑1 (MQ=255)
                      aGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGGCGCGCT‑CGTGGATTTCACGCAGCGCG‑CGATGc                                                                                                                            <  1:210951/150‑1 (MQ=38)
                       gggCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGCTGcc                                                                                                                           >  1:196483/1‑150 (MQ=255)
                        ggCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCg                                                                                                                          >  1:1238832/1‑150 (MQ=255)
                           aacaacTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGGCGCGCT‑CGTGGATTTCACGCAGCGCG‑CGATGcaccac                                                                                                                       <  2:1698508/150‑6 (MQ=37)
                           aacaacTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGGTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGGCGCGCT‑CGTGGATTTCACGCAGCGCG‑CGATGcaccac                                                                                                                       <  2:1691953/150‑6 (MQ=25)
                            acaacTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTtcg                                                                                                                      >  2:789893/1‑150 (MQ=255)
                            acaacTGAGGGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTGCCGGTTGATGTGGTGCATC‑GCGCGCTTCGTGAAATCCACG‑AGCGCGCCGATGCCGTtcg                                                                                                                      <  1:1274358/150‑1 (MQ=255)
                              aacTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTGGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTtcgtc                                                                                                                    <  1:1357971/150‑1 (MQ=255)
                               acTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTtcgtcg                                                                                                                   <  1:490072/150‑1 (MQ=255)
                               acTGAGTGAAGCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTtcgtcg                                                                                                                   <  2:1483902/150‑1 (MQ=255)
                                 tGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCgtcgtc                                                                                                                 >  1:426520/1‑150 (MQ=255)
                                 tGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCgtcgtc                                                                                                                 >  1:183299/1‑150 (MQ=255)
                                  gAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCgtcgtcc                                                                                                                >  1:279445/1‑149 (MQ=255)
                                   aGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAg                                                                                                               <  1:1250793/150‑1 (MQ=255)
                                    gTGAATCTGGCCTGGAAGTGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGc                                                                                                              >  2:215912/1‑150 (MQ=255)
                                     tGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCg                                                                                                             <  2:685726/150‑1 (MQ=255)
                                        aTCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGttt                                                                                                          >  2:892125/1‑150 (MQ=255)
                                          cTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTcc                                                                                                        >  2:246885/1‑150 (MQ=255)
                                          cTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGGTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTcc                                                                                                        <  2:553527/150‑1 (MQ=255)
                                             gCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGGTCGTCGTCAGCGTTTCCAgg                                                                                                     >  1:167621/1‑150 (MQ=255)
                                               cTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCg                                                                                                   >  2:116619/1‑150 (MQ=255)
                                                  gAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGAGGCGCTGGAAGGCAAACATCCGGCAGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGTAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATg                                                                                                <  1:1503727/150‑1 (MQ=255)
                                                    aGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgc                                                                                              >  1:1276797/1‑150 (MQ=255)
                                                      ggCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAATGCAAACATACGGCTGTACCCGTTGATGAGGTGCATC‑CCGCCCTGCGTGAAATCCACG‑AGCGCGCCGATCCCGTACTTCGTCAGCGTGTCCAGGCGATGgcag                                                                                            >  2:672493/1‑148 (MQ=255)
                                                      ggCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgctg                                                                                            >  1:61162/1‑150 (MQ=255)
                                                       gctgctGCCAGAAGCACTGGGTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgctgc                                                                                           <  1:1726774/150‑1 (MQ=255)
                                                        ctgctgCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGATTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgctgct                                                                                          <  2:387379/150‑1 (MQ=255)
                                                          gctgcCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGa                                                                                        >  1:1731038/1‑150 (MQ=255)
                                                          gctgcCAGAAGCACTGGTTGCAGGTAGTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGa                                                                                        <  1:1156923/150‑1 (MQ=255)
                                                           ctgcCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGaa                                                                                       >  1:521344/1‑150 (MQ=255)
                                                             gcCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAgg                                                                                     >  2:899083/1‑150 (MQ=255)
                                                             gcCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAgg                                                                                     <  2:613703/150‑1 (MQ=255)
                                                               cAGAAGCTCTGGTGGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGgtg                                                                                   <  2:481774/150‑1 (MQ=255)
                                                               cAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCTCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGgtg                                                                                   >  2:1113529/1‑150 (MQ=255)
                                                                   aGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGaa                                                                               <  2:100690/150‑1 (MQ=255)
                                                                    gCACTGGTTGCAGGTTTTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGaaa                                                                              <  1:357602/150‑1 (MQ=255)
                                                                     cACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGAGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAAt                                                                             <  2:279445/150‑1 (MQ=255)
                                                                      aCTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCAATGGCTGCTGAAGGTGTGAAATa                                                                            >  1:370623/1‑150 (MQ=255)
                                                                         ggTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCt                                                                         >  1:1243405/1‑150 (MQ=255)
                                                                          gTTGCAGGTATTGCCGATGCGCTGGCAGGCAAACATCCGGCTGTTCCGGTGGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTg                                                                        <  1:266429/150‑1 (MQ=255)
                                                                           ttGCAGGTATTGCCGTTGCGCTGGAAGGCATACATCCGGCTGTGCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTgg                                                                       <  2:939050/150‑1 (MQ=255)
                                                                              cAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGgaag                                                                    >  1:1025160/1‑150 (MQ=255)
                                                                                ggTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAgg                                                                                                     >  1:1251258/1‑115 (MQ=255)
                                                                                ggTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAgg                                                                                                     <  2:1251258/115‑1 (MQ=255)
                                                                                ggTAGTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGgaagaa                                                                  <  1:56367/150‑1 (MQ=255)
                                                                                    ttGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACg                                                              <  1:1590548/150‑1 (MQ=255)
                                                                                      gCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGcc                                                            >  1:1641175/1‑150 (MQ=255)
                                                                                      gCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGcc                                                            >  1:1664486/1‑150 (MQ=255)
                                                                                         gATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCaaa                                                         >  2:1195070/1‑150 (MQ=255)
                                                                                           tGCGCTGGAATGCAAACATCCGGCTGTTCCGGTTGATGTTGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCaaaaa                                                       <  2:387590/150‑1 (MQ=255)
                                                                                             cgcTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaag                                                     >  1:24778/1‑150 (MQ=255)
                                                                                              gcTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaaga                                                    >  2:748435/1‑150 (MQ=255)
                                                                                                tGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaagaag                                                  <  1:246885/150‑1 (MQ=255)
                                                                                                tGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaagaag                                                  <  2:1042164/150‑1 (MQ=255)
                                                                                                 ggAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAgg                                                 <  1:708023/150‑1 (MQ=255)
                                                                                                     ggCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGtgtg                                             <  2:872029/150‑1 (MQ=255)
                                                                                                     ggCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGtgtg                                             >  1:1673940/1‑150 (MQ=255)
                                                                                                     ggCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATCCCGTTCGTCGTCAGCGTTTCCATGCGATGGCTGCTGAAGGGGTGAACTACCTCCAAGCAAACCCCAACCACGCAGGtgtg                                             >  1:301213/1‑150 (MQ=255)
                                                                                                       cAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGaa                                           <  1:1686936/150‑1 (MQ=255)
                                                                                                       cAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGaa                                           >  2:427288/1‑150 (MQ=255)
                                                                                                        aaaCATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAAt                                          >  1:618551/1‑150 (MQ=255)
                                                                                                          aCATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAg                                        >  1:427764/1‑150 (MQ=255)
                                                                                                               cGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGTTGGCTGCTGAAGGTGTGAATTACCGGGAAGAAAACGCCAAAACAGAAGGTGGGCTTAGCCCCg                                   >  1:1136541/1‑150 (MQ=255)
                                                                                                               cGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCg                                   <  1:1113529/150‑1 (MQ=255)
                                                                                                                ggaacagccGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGa                                  >  2:218932/8‑150 (MQ=255)
                                                                                                                ggCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGa                                  >  1:1129311/1‑150 (MQ=255)
                                                                                                                    gTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCTCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATCGCACCGAATCt                              >  2:1639460/1‑150 (MQ=255)
                                                                                                                    gTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCt                              <  1:1483105/150‑1 (MQ=255)
                                                                                                                    gTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCt                              <  1:1347366/150‑1 (MQ=255)
                                                                                                                     gtCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCAGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTg                             >  1:632118/2‑150 (MQ=255)
                                                                                                                        cGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCCCG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGAGGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGcc                          <  2:130917/150‑1 (MQ=255)
                                                                                                                          gTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGGCAGCGTTTCCAGGCGAGGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTg                        <  1:915435/150‑1 (MQ=255)
                                                                                                                           ttGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGc                       <  1:806836/150‑1 (MQ=255)
                                                                                                                            tGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑CGCGCGGCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCa                      >  1:695063/1‑150 (MQ=255)
                                                                                                                             gATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCaa                     >  2:486667/1‑150 (MQ=255)
                                                                                                                                  ggtgCATC‑GCGCGCTGAGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCATAAAAGAAGGTGTTAATAGCACCGATTCTGGCCTGCAATTCcg                <  1:368679/150‑1 (MQ=255)
                                                                                                                                   gtgCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCcgc               >  2:1239503/1‑150 (MQ=255)
                                                                                                                                     gCATC‑GCGCGCTGCGTGAAAACCACG‑AGCGCGCCGATGCCGTTCGGCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGAGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGt             <  1:599487/150‑1 (MQ=255)
                                                                                                                                      catcggcgcgcT‑CGTGGATTTCACGCAGCGCG‑CGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTg            >  1:246269/6‑150 (MQ=38)
                                                                                                                                       atcggcgcgcT‑CGTGGATTTCACGCAGCGCGAAGATGCCGTTCGTCGTCATAGTTTCCAGGCGATGGCTGCTGAAGGTGTGAATTACCTGGAAGAAAACGACAAAATAGAAGGTGTGAATATCAACGAATCTGTCCTGCAATTCCTCGTg            >  1:954250/5‑150 (MQ=255)
                                                                                                                                       aTC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGa           <  2:521344/150‑1 (MQ=255)
                                                                                                                                       aTC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGa           >  2:326447/1‑150 (MQ=255)
                                                                                                                                        tC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGAt          <  2:773501/150‑1 (MQ=255)
                                                                                                                                        gc‑gcgcgCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGAt          <  2:719151/149‑1 (MQ=255)
                                                                                                                                           gcgcgcTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCa        >  2:421023/1‑150 (MQ=255)
                                                                                                                                              cgcTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAAcc     <  1:700409/150‑1 (MQ=255)
                                                                                                                                               gcTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCa    >  1:1454549/1‑150 (MQ=255)
                                                                                                                                                cTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAg   <  2:1731038/150‑1 (MQ=255)
                                                                                                                                                 ggCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAgg  <  1:215912/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCTACGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGG  >  NC_000913/4428975‑4429268

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