Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,429,129 C→T 5.8% I65I (ATC→ATT) ‡ 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,1290CT 5.8% 384.3 / 7.9 121I65I (ATC→ATT) ‡fklBFKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base C (60/54);  new base T (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 = 7.01e-01

Only 100 of 122 total aligned reads displayed.

GGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTGAAGG  >  NC_000913/4428980‑4429275
                                                                                                                                                      |                                                                                                                                                   
      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)
                      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:613703/150‑1 (MQ=255)
                                                        gcCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAgg                                                                                            >  2:899083/1‑150 (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)
                                                                 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:1641175/1‑150 (MQ=255)
                                                                                 gCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGcc                                                                   >  1:1664486/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                                                         <  2:1042164/150‑1 (MQ=255)
                                                                                           tGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAaagaag                                                         <  1:246885/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)
                                                                                                  cAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGaa                                                  >  2:427288/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:326447/1‑150 (MQ=255)
                                                                                                                                  aTCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGa                  <  2:521344/150‑1 (MQ=255)
                                                                                                                                   tCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGAt                 <  2:773501/150‑1 (MQ=255)
                                                                                                                                          gcTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCa           >  1:1454549/1‑150 (MQ=255)
                                                                                                                                           cTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAg          <  2:1731038/150‑1 (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)
                                                                                                                                                   aaTCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTGAAgg  <  2:167621/150‑1 (MQ=255)
                                                                                                                                                   aaTCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTGCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCTATTCCGCGTGATCAACCAGGGTGAAgg  <  1:899083/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GGCATTGGTTTGCAGGTAGGGCAACAACTGAGTGAATCTGGCCTGGAAGGGCTGCTGCCAGAAGCACTGGTTGCAGGTATTGCCGATGCGCTGGAAGGCAAACATCCGGCTGTTCCGGTTGATGTGGTGCATCGCGCGCTGCGTGAAATCCACGAGCGCGCCGATGCCGTTCGTCGTCAGCGTTTCCAGGCGATGGCTGCTGAAGGTGTGAAATACCTGGAAGAAAACGCCAAAAAAGAAGGTGTGAATAGCACCGAATCTGGCCTGCAATTCCGCGTGATCAACCAGGGTGAAGG  >  NC_000913/4428980‑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: