Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,476,685 C→T 5.8% L245L (CTG→CTA)  fkpA ← 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_0009133,476,6850CT 5.8% 280.6 / 5.6 86L245L (CTG→CTA) fkpAFKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base C (44/37);  new base T (2/3);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

TTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGT  >  NC_000913/3476538‑3476833
                                                                                                                                                   |                                                                                                                                                    
ttCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCcagca                                                                                                                                                     <  3:714157/149‑1 (MQ=255)
  ccAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCagcagc                                                                                                                                                   <  3:416543/149‑1 (MQ=255)
   cAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGct                                                                                                                                                  >  3:442108/1‑149 (MQ=255)
         tAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACg                                                                                                                                             <  1:288460/148‑1 (MQ=255)
          aaTTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGt                                                                                                                                            >  1:62219/1‑148 (MQ=255)
             tataTCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCaaa                                                                                                                                        <  4:113512/149‑1 (MQ=255)
             tataTCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCaaa                                                                                                                                        <  4:510006/149‑1 (MQ=255)
             tataTCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCaa                                                                                                                                         <  1:346134/148‑1 (MQ=255)
                    ggCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGTAGCTCTACGTCAAACACCAgg                                                                                                                                 <  4:238959/149‑1 (MQ=255)
                    ggCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAg                                                                                                                                  <  2:886944/148‑1 (MQ=255)
                       cTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTa                                                                                                                              <  3:217493/149‑1 (MQ=255)
                             taaAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCATATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTc                                                                                                                        >  3:622846/1‑149 (MQ=255)
                              aaaaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGaa                                                                                                                           >  4:435820/1‑145 (MQ=255)
                              aaaaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGaa                                                                                                                           <  3:435820/145‑1 (MQ=255)
                              aaaaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTc                                                                                                                        <  1:957469/148‑1 (MQ=255)
                              aaaaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCg                                                                                                                       >  3:153956/1‑149 (MQ=255)
                               taagaaCCGCCGCCTGACCAGGCGGCGGTTTTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCg                                                                                                                       >  2:716435/5‑148 (MQ=255)
                                aaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCgg                                                                                                                      <  2:400966/148‑1 (MQ=255)
                                aaaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGt                                                                                                                     >  3:739393/1‑149 (MQ=255)
                                 aaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGt                                                                                                                     <  2:226160/148‑1 (MQ=255)
                                 aaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTg                                                                                                                    <  3:392510/149‑1 (MQ=255)
                                 aaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTg                                                                                                                    <  4:739393/149‑1 (MQ=255)
                                 aaaaCCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTCTTAGCATAATCTGCGGCTTTCGCGTCTGCTTCCGGCTTTGCATTAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTGCGTCAAACACCAGGAGAGAATTCGGTg                                                                                                                    >  4:209130/1‑149 (MQ=255)
                                   aaCCGCCGCCTGACCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGctc                                                                                                                                                 >  1:917132/1‑118 (MQ=255)
                                   aaCCGCCGCCTGACCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGctc                                                                                                                                                 <  2:917132/118‑1 (MQ=255)
                                     ccgccgCCTGACCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGa                                                                                                                 >  2:23122/1‑148 (MQ=255)
                                                tCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCagcagc                                                                                                                                                   <  1:917811/103‑1 (MQ=255)
                                                tCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCagcagc                                                                                                                                                   >  2:917811/1‑103 (MQ=255)
                                                    gccgcggTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCg                                                                                                  <  1:647070/145‑1 (MQ=255)
                                                         ggTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGc                                                                                            <  4:202699/149‑1 (MQ=255)
                                                            tCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGTCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCg                                                                                          >  2:95916/1‑148 (MQ=255)
                                                             cTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTa                                                                                        >  4:745866/1‑149 (MQ=255)
                                                              ttAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTaa                                                                                       >  3:614841/1‑149 (MQ=255)
                                                              ttAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTaa                                                                                       >  3:524651/1‑149 (MQ=255)
                                                                   gCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCg                                                                                                                       <  3:273568/112‑1 (MQ=255)
                                                                   gCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCg                                                                                                                       >  4:273568/1‑112 (MQ=255)
                                                                   gCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAg                                                                                  >  4:162723/1‑149 (MQ=255)
                                                                     ttATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGt                                                                                 >  1:630255/1‑148 (MQ=255)
                                                                      tATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGtt                                                                                >  2:595281/1‑148 (MQ=255)
                                                                      tATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTc                                                                               <  4:268409/149‑1 (MQ=255)
                                                                       aTTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCt                                                                              <  4:324382/149‑1 (MQ=255)
                                                                          ttttAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCtggt                                                                           <  3:360747/149‑1 (MQ=255)
                                                                          ttttAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCtgg                                                                            <  2:253403/148‑1 (MQ=255)
                                                                               gCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAAt                                                                      >  3:817550/1‑149 (MQ=255)
                                                                                cAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATa                                                                     >  3:597266/1‑149 (MQ=255)
                                                                                 aGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATaa                                                                    <  4:442108/149‑1 (MQ=255)
                                                                                        gCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGtt                                                              >  2:465353/1‑148 (MQ=255)
                                                                                          ggCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTg                                                            >  1:447832/1‑148 (MQ=255)
                                                                                           gCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCGCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGAt                                                          <  3:84462/149‑1 (MQ=255)
                                                                                           gCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGa                                                           <  1:465353/148‑1 (MQ=255)
                                                                                           gCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGAt                                                          <  3:545571/149‑1 (MQ=255)
                                                                                            cTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGAt                                                          <  2:447832/148‑1 (MQ=255)
                                                                                             tttCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATc                                                                                                               >  3:499260/1‑94 (MQ=255)
                                                                                             tttCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATc                                                                                                               <  4:499260/94‑1 (MQ=255)
                                                                                             tttCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCt                                                        >  3:666840/1‑149 (MQ=255)
                                                                                              ttCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAg                                                                >  4:366894/1‑140 (MQ=255)
                                                                                              ttCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAg                                                                <  3:366894/140‑1 (MQ=255)
                                                                                                cGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTAc                                                     <  4:597266/149‑1 (MQ=255)
                                                                                                         ttCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAg                                                                                                                                  >  2:222926/1‑63 (MQ=255)
                                                                                                         ttCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAg                                                                                                                                  <  1:222926/63‑1 (MQ=255)
                                                                                                         ttCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCt                                            >  4:105506/1‑149 (MQ=255)
                                                                                                         ttCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCt                                            >  3:406948/1‑149 (MQ=255)
                                                                                                            cGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTg                                          <  2:941656/148‑1 (MQ=255)
                                                                                                             ggCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGTGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGa                                         <  2:920452/148‑1 (MQ=255)
                                                                                                             ggCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGa                                         >  1:465898/1‑148 (MQ=255)
                                                                                                                tttGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGtt                                     >  4:147470/1‑149 (MQ=255)
                                                                                                                tttGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGt                                      <  1:792739/148‑1 (MQ=255)
                                                                                                                    cATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTa                                                                                        >  1:595413/1‑94 (MQ=255)
                                                                                                                    cATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTa                                                                                        <  2:595413/94‑1 (MQ=255)
                                                                                                                    cATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGttcttc                                 <  4:524651/149‑1 (MQ=255)
                                                                                                                     aTCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCa                                >  4:888328/1‑149 (MQ=255)
                                                                                                                     aTCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCa                                >  3:269483/1‑149 (MQ=255)
                                                                                                                      tCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATaa                                                                    <  1:443408/112‑1 (MQ=255)
                                                                                                                      tCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATaa                                                                    >  2:443408/1‑112 (MQ=255)
                                                                                                                      tCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTg                                          <  2:569739/138‑1 (MQ=255)
                                                                                                                      tCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTg                                          >  1:569739/1‑138 (MQ=255)
                                                                                                                      tCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAg                               >  3:32477/1‑149 (MQ=255)
                                                                                                                       cAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAg                               >  2:338166/1‑148 (MQ=255)
                                                                                                                          ccTTCGGCGCTGGTTTCACATCCAGTAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGAt                                                          >  4:215118/1‑118 (MQ=255)
                                                                                                                          ccTTCGGCGCTGGTTTCACATCCAGTAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGAt                                                          <  3:215118/118‑1 (MQ=255)
                                                                                                                           cTTCGGCGCTGGTTTCACATCCAGTAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGa                                                           >  2:839197/1‑116 (MQ=255)
                                                                                                                           cTTCGGCGCTGGTTTCACATCCAGTAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGa                                                           <  1:839197/116‑1 (MQ=255)
                                                                                                                                   cTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGATTTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCGGGTGGAATAACCAGTTTGCTCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCaa                  >  4:687095/1‑149 (MQ=255)
                                                                                                                                    tGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCaa                  >  2:80099/1‑148 (MQ=255)
                                                                                                                                      gTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAAccc               >  4:189144/1‑149 (MQ=255)
                                                                                                                                       nnncacaTCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAAccc               >  2:122989/4‑148 (MQ=255)
                                                                                                                                            caTCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGAt          <  2:66952/148‑1 (MQ=255)
                                                                                                                                              tCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATaa        >  2:419209/1‑148 (MQ=255)
                                                                                                                                              tCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATaa        >  1:662014/1‑148 (MQ=255)
                                                                                                                                              tCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTGTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAac       >  3:853507/1‑149 (MQ=255)
                                                                                                                                                   cagcTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACAccat  >  3:571435/1‑147 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTCCAGCACTAATTATATCCGGCCTGTAATAAAAAAACCGCCGCCTGGTCAGGCGGCGGTTCTTAATGCTTATTTTTTAGCAGAATCTGCGGCTTTCGCATCAGCTTCCGGCTTTGCATCAGCCTTCGGCGCTGGTTTCACATCCAGCAGCTCTACGTCAAACACCAGGGTAGAATTCGGTGGGATCCCCGGAACACCCGCTTTGCCGTAAGCCAGTTCTGGTGGAATAACCAGTTTGATCTTACCGCCTTTCTTGATGTTCTTCAGACCTTCTGTCCAACCCGGGATAACACCGT  >  NC_000913/3476538‑3476833

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑