Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,429,133:1 +C 5.8% coding (199/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,1331.C 5.8% 501.5 / 23.5 121coding (199/621 nt)fklBFKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base . (60/54);  new base C (2/5);  total (62/59)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 121 total aligned reads displayed.

GGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTGAAGG  >  NC_000913/4428986‑4429275
                                                                                                                                                     |                                                                                                                                              
ggTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAATTGCTCCTGCGCGAAGAACTGGTTGCAGGTATTGCGGATGCTCTGGAAGGCAAAAATCAGGCTGATAAGGTTGATGTGGTGCATC‑GCGCTCTGCGAGAAATCCACG‑Ag                                                                                                                                              >  1:238891/1‑150 (MQ=255)
ggTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑Ag                                                                                                                                              >  1:939050/1‑150 (MQ=255)
ggTGTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCGGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑Ag                                                                                                                                              <  1:976236/150‑1 (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:210951/150‑1 (MQ=38)
           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)
                 acaacTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTtcg                                                                                                                             >  2:789893/1‑150 (MQ=255)
                 acaacTGAGGGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTGCCGGTTGATGTGGTGCATC‑GCGCGCTTCGTGAAATCCACG‑AGCGCGCCGATGCCGTtcg                                                                                                                             <  1:1274358/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:183299/1‑150 (MQ=255)
                      tGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCgtcgtc                                                                                                                        >  1:426520/1‑150 (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‑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)
                                    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)
                                             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)
                                                  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)
                                                     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)
                                                                         ttGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACg                                                                     <  1:1590548/150‑1 (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                                                         <  2:1042164/150‑1 (MQ=255)
                                                                                     tGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaagaag                                                         <  1:246885/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)
                                                                                            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)
                                                                                                     ggCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGa                                         >  2:669456/1‑150 (MQ=255)
                                                                                                         gTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCTCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATCGCACCGAATCt                                     >  2:1639460/1‑150 (MQ=255)
                                                                                                         gTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCt                                     <  1:1347366/150‑1 (MQ=255)
                                                                                                         gTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCt                                     <  1:1483105/150‑1 (MQ=255)
                                                                                                          gtCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCAGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTg                                    >  1:632118/2‑150 (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)
                                                                                                                       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)
                                                                                                                            aTC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGa                  >  2:326447/1‑150 (MQ=255)
                                                                                                                             tC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGAt                 <  2:773501/150‑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)
                                                                                                                                       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)
                                                                                                                                             aaTCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTGAAgg  <  2:167621/150‑1 (MQ=255)
                                                                                                                                             aaTCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTGCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCTATTCCGCGTGATCAACCAGGGTGAAgg  <  1:899083/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
GGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATC‑GCGCGCTGCGTGAAATCCACG‑AGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTGAAGG  >  NC_000913/4428986‑4429275

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