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.

TACGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTGAA  >  NC_000913/4428977‑4429273
                                                                                                                                                     |                                                                                                                                                     
tACGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGaa                                                                                                                                                      >  2:496225/1‑150 (MQ=255)
   ggCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATc                                                                                                                                                   >  1:872029/1‑150 (MQ=255)
         ggTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAATTGCTCCTGCGCGAAGAACTGGTTGCAGGTATTGCGGATGCTCTGGAAGGCAAAAATCAGGCTGATAAGGTTGATGTGGTGCATCGCGCTCTGCGAGAAATCCACGAg                                                                                                                                            >  1:238891/1‑150 (MQ=255)
         ggTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAg                                                                                                                                            >  1:939050/1‑150 (MQ=255)
           gttGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAgcg                                                                                                                                          <  1:1298192/149‑1 (MQ=255)
                aGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAgcgcgcga                                                                                                                                     <  1:735538/150‑3 (MQ=255)
                aGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGGTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGGCGCGCTCGTGGATTTCACGCAgcgcgcg                                                                                                                                      <  2:907940/150‑2 (MQ=255)
                  gTAGGGCAACAACTGTGTGAATCTGGCCTGTACGGGCTGCTGCCAGAATCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGAt                                                                                                                                   <  1:1131021/150‑1 (MQ=255)
                  gTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGAt                                                                                                                                   <  1:625323/150‑1 (MQ=255)
                   tAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATg                                                                                                                                  >  1:96521/1‑150 (MQ=255)
                   gaGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCTGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATg                                                                                                                                  <  2:147649/149‑1 (MQ=255)
                    aGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGGCGCGCTCGTGGATTTCACGCAGCGCGCGATGc                                                                                                                                 <  1:1509801/150‑1 (MQ=38)
                     gggCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGCTGcc                                                                                                                                >  1:196483/1‑150 (MQ=255)
                      ggCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCg                                                                                                                               >  1:1238832/1‑150 (MQ=255)
                         aacaacTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGGCGCGCTCGTGGATTTCACGCAGCGCGCGATGcaccac                                                                                                                            <  2:1698508/150‑6 (MQ=37)
                         aacaacTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGGTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGGCGCGCTCGTGGATTTCACGCAGCGCGCGATGcaccac                                                                                                                            <  2:1691953/150‑6 (MQ=25)
                          acaacTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTtcg                                                                                                                           >  2:789893/1‑150 (MQ=255)
                          acaacTGAGGGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTGCCGGTTGATGTGGTGCATCGCGCGCTTCGTGAAATCCACGAGCGCGCCGATGCCGTtcg                                                                                                                           <  1:1274358/150‑1 (MQ=255)
                            aacTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTGGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTtcgtc                                                                                                                         <  1:1357971/150‑1 (MQ=255)
                               tGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCgtcgtc                                                                                                                      >  1:426520/1‑150 (MQ=255)
                               tGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCgtcgtc                                                                                                                      >  1:183299/1‑150 (MQ=255)
                                gAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCgtcgtcc                                                                                                                     >  1:279445/1‑149 (MQ=255)
                                  gTGAATCTGGCCTGGAAGTGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGc                                                                                                                   >  2:215912/1‑150 (MQ=255)
                                   tGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCg                                                                                                                  <  2:685726/150‑1 (MQ=255)
                                     aaTCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGAGtt                                                                                                                >  2:253906/1‑150 (MQ=255)
                                      aTCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGttt                                                                                                               >  2:892125/1‑150 (MQ=255)
                                        cTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTcc                                                                                                             >  2:246885/1‑150 (MQ=255)
                                           gCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAgg                                                                                                          >  2:1686936/1‑150 (MQ=255)
                                           gCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGGTCGTCGTCAGCGTTTCCAgg                                                                                                          >  1:167621/1‑150 (MQ=255)
                                             cTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCg                                                                                                        <  2:804448/150‑1 (MQ=255)
                                                gAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGAGGCGCTGGAAGGCAAACATCCGGCAGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGTAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATg                                                                                                     <  1:1503727/150‑1 (MQ=255)
                                                  aGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgc                                                                                                   >  1:1276797/1‑150 (MQ=255)
                                                    ggCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAATGCAAACATACGGCTGTACCCGTTGATGAGGTGCATCCCGCCCTGCGTGAAATCCACGAGCGCGCCGATCCCGTACTTCGTCAGCGTGTCCAGGCGATGgcag                                                                                                 >  2:672493/1‑148 (MQ=255)
                                                    ggCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgctg                                                                                                 >  1:773501/1‑150 (MQ=255)
                                                    ggCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgctg                                                                                                 >  1:61162/1‑150 (MQ=255)
                                                     gctgctGCCAGAAGCACTGGGTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgctgc                                                                                                <  1:1726774/150‑1 (MQ=255)
                                                      ctgctgCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgctgct                                                                                               <  1:328087/150‑1 (MQ=255)
                                                      ctgctgCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGATTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGgctgct                                                                                               <  2:387379/150‑1 (MQ=255)
                                                        gctgcCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGa                                                                                             >  1:1731038/1‑150 (MQ=255)
                                                         ctgcCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGaa                                                                                            >  1:521344/1‑150 (MQ=255)
                                                          tgcCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAg                                                                                           >  1:1689896/1‑150 (MQ=255)
                                                           gcCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAgg                                                                                          >  2:899083/1‑150 (MQ=255)
                                                           gcCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAgg                                                                                          <  2:613703/150‑1 (MQ=255)
                                                             cAGAAGCTCTGGTGGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGgtg                                                                                        <  2:481774/150‑1 (MQ=255)
                                                             cAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCTCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGgtg                                                                                        >  2:1113529/1‑150 (MQ=255)
                                                              aGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGgtgt                                                                                       >  1:1222075/1‑150 (MQ=255)
                                                                 aGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGaa                                                                                    <  2:100690/150‑1 (MQ=255)
                                                                  gCACTGGTTGCAGGTTTTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGaaa                                                                                   <  1:357602/150‑1 (MQ=255)
                                                                   cACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGAGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAAt                                                                                  <  2:279445/150‑1 (MQ=255)
                                                                    aCTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCAATGGCTGCTGAAGGTGTGAAATa                                                                                 >  1:370623/1‑150 (MQ=255)
                                                                       ggTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCt                                                                              >  1:1243405/1‑150 (MQ=255)
                                                                        gTTGCAGGTATTGCCGATGCGCTGGCAGGCAAACATCCGGCTGTTCCGGTGGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTg                                                                             <  1:266429/150‑1 (MQ=255)
                                                                         ttGCAGGTATTGCCGTTGCGCTGGAAGGCATACATCCGGCTGTGCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTgg                                                                            <  2:939050/150‑1 (MQ=255)
                                                                            cAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGgaag                                                                         >  1:1025160/1‑150 (MQ=255)
                                                                              ggTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAgg                                                                                                          >  1:1251258/1‑115 (MQ=255)
                                                                              ggTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAgg                                                                                                          <  2:1251258/115‑1 (MQ=255)
                                                                              ggTAGTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGgaagaa                                                                       <  1:56367/150‑1 (MQ=255)
                                                                                    gCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGcc                                                                 >  1:1664486/1‑150 (MQ=255)
                                                                                    gCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGcc                                                                 >  1:1641175/1‑150 (MQ=255)
                                                                                       gATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCaaa                                                              >  2:1195070/1‑150 (MQ=255)
                                                                                         tGCGCTGGAATGCAAACATCCGGCTGTTCCGGTTGATGTTGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCaaaaa                                                            <  2:387590/150‑1 (MQ=255)
                                                                                            gcTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaaga                                                         >  2:748435/1‑150 (MQ=255)
                                                                                             cTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGAGCATCGCTCGCTGCGTGAAATCCAAGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaagaa                                                        >  1:62604/1‑150 (MQ=255)
                                                                                              tGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaagaag                                                       <  1:246885/150‑1 (MQ=255)
                                                                                              tGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaagaag                                                       <  2:1042164/150‑1 (MQ=255)
                                                                                               ggAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAgg                                                      <  1:708023/150‑1 (MQ=255)
                                                                                                   ggCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGtgtg                                                  >  1:1673940/1‑150 (MQ=255)
                                                                                                   ggCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATCCCGTTCGTCGTCAGCGTTTCCATGCGATGGCTGCTGAAGGGGTGAACTACCTCCAAGCAAACCCCAACCACGCAGGtgtg                                                  >  1:301213/1‑150 (MQ=255)
                                                                                                        aCATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAg                                             >  1:427764/1‑150 (MQ=255)
                                                                                                             cGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGTTGGCTGCTGAAGGTGTGAATTACCGGGAAGAAAACGCCAAAACAGAAGGTGGGCTTAGCCCCg                                        >  1:1136541/1‑150 (MQ=255)
                                                                                                             cGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCg                                        <  1:1113529/150‑1 (MQ=255)
                                                                                                              ggaacagccGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGa                                       >  2:218932/8‑150 (MQ=255)
                                                                                                              ggCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGa                                       >  1:1129311/1‑150 (MQ=255)
                                                                                                              ggCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGa                                       >  2:669456/1‑150 (MQ=255)
                                                                                                                  gTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCTCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATCGCACCGAATCt                                   >  2:1639460/1‑150 (MQ=255)
                                                                                                                  gTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCt                                   <  1:1483105/150‑1 (MQ=255)
                                                                                                                   gtCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCAGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTg                                  >  1:632118/2‑150 (MQ=255)
                                                                                                                      cGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCCCGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGAGGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGcc                               <  2:130917/150‑1 (MQ=255)
                                                                                                                      cGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGcc                               <  1:496225/150‑1 (MQ=255)
                                                                                                                        gTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGGCAGCGTTTCCAGGCGAGGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTg                             <  1:915435/150‑1 (MQ=255)
                                                                                                                         ttGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGc                            <  1:806836/150‑1 (MQ=255)
                                                                                                                          tGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGCGCGCGGCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCa                           >  1:695063/1‑150 (MQ=255)
                                                                                                                           gATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCaa                          >  2:486667/1‑150 (MQ=255)
                                                                                                                            aTGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAAAACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTCGCCTGCAAt                         >  2:467315/1‑150 (MQ=255)
                                                                                                                                 gtgCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCcgc                    >  2:1239503/1‑150 (MQ=255)
                                                                                                                                   gCATCGCGCGCTGCGTGAAAACCACGAGCGCGCCGATGCCGTTCGGCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGAGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGt                  <  1:599487/150‑1 (MQ=255)
                                                                                                                                    catcggcgcgcTCGTGGATTTCACGCAGCGCGCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTg                 >  1:246269/6‑150 (MQ=38)
                                                                                                                                     atcggcgcgcTCGTGGATTTCACGCAGCGCGAAGATGCCGTTCGTCGTCATAGTTTCCAGGCGATGGCTGCTGAAGGTGTGAATTACCTGGAAGAAAACGACAAAATAGAAGGTGTGAATATCAACGAATCTGTCCTGCAATTCCTCGTg                 >  1:954250/5‑150 (MQ=255)
                                                                                                                                     aTCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGa                <  2:521344/150‑1 (MQ=255)
                                                                                                                                      tCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGAt               <  2:773501/150‑1 (MQ=255)
                                                                                                                                         gcgcgcTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCa             >  2:421023/1‑150 (MQ=255)
                                                                                                                                             gcTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCa         >  1:1454549/1‑150 (MQ=255)
                                                                                                                                               ggCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAgg       <  1:215912/149‑1 (MQ=255)
                                                                                                                                                gagTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCTGGCGATGGCTGAGGCAGGTGTGAAATACCTGGAAGACAACGCCAAAAAAGCAGGTGTGAATAGGAACGGATCTGGCCGGCCAGGCCGCGTGCGCATCCAggg      >  2:836438/3‑150 (MQ=255)
                                                                                                                                                 ggTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTTAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAAGGt     <  1:116619/149‑1 (MQ=255)
                                                                                                                                                  gTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTg    >  1:814802/1‑150 (MQ=255)
                                                                                                                                                  gTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGAGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTg    >  2:1116009/1‑150 (MQ=255)
                                                                                                                                                  gTGAAATCCACGAGCGCGCCGAGGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTg    <  1:1165511/150‑1 (MQ=255)
                                                                                                                                                   tGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTGa   >  1:422290/1‑150 (MQ=255)
                                                                                                                                                    gAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTGaa  >  2:40152/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACGGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTGAA  >  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: