Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,429,125 A→G 5.8% 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 5.8% 360.2 / 9.0 121E64G (GAA→GGA) fklBFKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base A (61/52);  new base G (2/5);  total (63/58)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

Only 100 of 122 total aligned reads displayed.

TACGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTGAA  >  NC_000913/4428977‑4429273
                                                                                                                                                     |                                                                                                                                                     
tACGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGaa                                                                                                                                                      >  2:496225/1‑150 (MQ=255)
   ggCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATc                                                                                                                                                   >  1:872029/1‑150 (MQ=255)
         ggTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAATTGCTCCTGCGCGAAGAACTGGTTGCAGGTATTGCGGATGCTCTGGAAGGCAAAAATCAGGCTGATAAGGTTGATGTGGTGCATC‑GCGCTCTGCGAGAAATCCACG‑Ag                                                                                                                                            >  1:238891/1‑150 (MQ=255)
         ggTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑Ag                                                                                                                                            >  1:939050/1‑150 (MQ=255)
           gttGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑Agcg                                                                                                                                          <  1:1298192/149‑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)
                   tAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATg                                                                                                                                  >  1:96521/1‑150 (MQ=255)
                   gaGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCTGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATg                                                                                                                                  <  2:147649/149‑1 (MQ=255)
                    aGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGGCGCGCT‑CGTGGATTTCACGCAGCGCG‑CGATGc                                                                                                                                 <  1:1509801/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)
                               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)
                                  gTGAATCTGGCCTGGAAGTGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGc                                                                                                                   >  2:215912/1‑150 (MQ=255)
                                   tGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCg                                                                                                                  <  2:685726/150‑1 (MQ=255)
                                     aaTCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGAGtt                                                                                                                >  2:253906/1‑150 (MQ=255)
                                      aTCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGttt                                                                                                               >  2:892125/1‑150 (MQ=255)
                                        cTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTcc                                                                                                             >  2:246885/1‑150 (MQ=255)
                                           gCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAgg                                                                                                          >  2:1686936/1‑150 (MQ=255)
                                           gCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGGTCGTCGTCAGCGTTTCCAgg                                                                                                          >  1:167621/1‑150 (MQ=255)
                                             cTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCg                                                                                                        <  2:804448/150‑1 (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:773501/1‑150 (MQ=255)
                                                    ggCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgctg                                                                                                 >  1:61162/1‑150 (MQ=255)
                                                     gctgctGCCAGAAGCACTGGGTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgctgc                                                                                                <  1:1726774/150‑1 (MQ=255)
                                                      ctgctgCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgctgct                                                                                               <  1:328087/150‑1 (MQ=255)
                                                      ctgctgCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGATTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgctgct                                                                                               <  2:387379/150‑1 (MQ=255)
                                                        gctgcCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGa                                                                                             >  1:1731038/1‑150 (MQ=255)
                                                         ctgcCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGaa                                                                                            >  1:521344/1‑150 (MQ=255)
                                                          tgcCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAg                                                                                           >  1:1689896/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)
                                                              aGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGgtgt                                                                                       >  1:1222075/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)
                                                                                    gCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGcc                                                                 >  1:1664486/1‑150 (MQ=255)
                                                                                    gCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGcc                                                                 >  1:1641175/1‑150 (MQ=255)
                                                                                       gATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCaaa                                                              >  2:1195070/1‑150 (MQ=255)
                                                                                         tGCGCTGGAATGCAAACATCCGGCTGTTCCGGTTGATGTTGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCaaaaa                                                            <  2:387590/150‑1 (MQ=255)
                                                                                            gcTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaaga                                                         >  2:748435/1‑150 (MQ=255)
                                                                                             cTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGAGCATC‑GCTCGCTGCGTGAAATCCAAG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaagaa                                                        >  1:62604/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                                                  >  1:1673940/1‑150 (MQ=255)
                                                                                                   ggCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATCCCGTTCGTCGTCAGCGTTTCCATGCGATGGCTGCTGAAGGGGTGAACTACCTCCAAGCAAACCCCAACCACGCAGGtgtg                                                  >  1:301213/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)
                                                                                                              ggCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGa                                       >  2:669456/1‑150 (MQ=255)
                                                                                                                  gTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCTCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATCGCACCGAATCt                                   >  2:1639460/1‑150 (MQ=255)
                                                                                                                  gTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCt                                   <  1:1483105/150‑1 (MQ=255)
                                                                                                                   gtCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCAGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTg                                  >  1:632118/2‑150 (MQ=255)
                                                                                                                      cGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCCCG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGAGGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGcc                               <  2:130917/150‑1 (MQ=255)
                                                                                                                      cGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGcc                               <  1:496225/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)
                                                                                                                            aTGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAAAACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTCGCCTGCAAt                         >  2:467315/1‑150 (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)
                                                                                                                                      tC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGAt               <  2:773501/150‑1 (MQ=255)
                                                                                                                                         gcgcgcTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCa             >  2:421023/1‑150 (MQ=255)
                                                                                                                                             gcTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCa         >  1:1454549/1‑150 (MQ=255)
                                                                                                                                               ggCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAgg       <  1:215912/149‑1 (MQ=255)
                                                                                                                                                gagTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCTGGCGATGGCTGAGGCAGGTGTGAAATACCTGGAAGACAACGCCAAAAAAGCAGGTGTGAATAGGAACGGATCTGGCCGGCCAGGCCGCGTGCGCATCCAggg      >  2:836438/3‑150 (MQ=255)
                                                                                                                                                 ggTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTTAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAAGGt     <  1:116619/149‑1 (MQ=255)
                                                                                                                                                  gTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTg    >  1:814802/1‑150 (MQ=255)
                                                                                                                                                  gTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGAGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTg    >  2:1116009/1‑150 (MQ=255)
                                                                                                                                                  gTGAAATCCACG‑AGCGCGCCGAGGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTg    <  1:1165511/150‑1 (MQ=255)
                                                                                                                                                   tGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTGa   >  1:422290/1‑150 (MQ=255)
                                                                                                                                                    gAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTGaa  >  2:40152/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTGAA  >  NC_000913/4428977‑4429273

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: ‑