Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,518,082 C→G 9.1% noncoding (35/76 nt) alaX ← tRNA‑Ala

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,518,0820CG 9.1% 299.1 / 18.4 88noncoding (35/76 nt)alaXtRNA‑Ala
Reads supporting (aligned to +/- strand):  ref base C (34/46);  new base G (3/5);  total (37/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

CAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGC  >  NC_000913/2517935‑2518228
                                                                                                                                                   |                                                                                                                                                  
cAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcc                                                                                                                                                   >  4:802943/1‑149 (MQ=255)
cAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcc                                                                                                                                                   <  3:692480/149‑1 (MQ=255)
 aaTGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcca                                                                                                                                                  >  4:670345/1‑149 (MQ=255)
 aaTGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcca                                                                                                                                                  <  4:131327/149‑1 (MQ=255)
        cTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGCGACTAAAAACCGCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGACCCGCTGACCTCTTGCATGCCATGCAAgc                                                                                                                                           <  2:638341/149‑1 (MQ=255)
        cTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAga                                                                                                                                           <  1:237427/149‑2 (MQ=255)
        cTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAga                                                                                                                                           <  4:713117/149‑2 (MQ=255)
         tttCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGGAGAGAAAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAgag                                                                                                                                          <  4:232374/149‑3 (MQ=255)
         tttCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAgag                                                                                                                                          <  4:622257/149‑3 (MQ=255)
         tttCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAgag                                                                                                                                          <  1:156043/149‑3 (MQ=255)
         tttCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAgag                                                                                                                                          >  3:622257/1‑147 (MQ=255)
         tttCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAgcg                                                                                                                                          <  2:420307/149‑1 (MQ=255)
         tttCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAgcg                                                                                                                                          >  2:521919/1‑149 (MQ=255)
             aCCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGctct                                                                                                                                      <  3:291491/149‑1 (MQ=255)
               cAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTcc                                                                                                                                    <  2:720296/149‑1 (MQ=255)
                    ggcACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCt                                                                                                                               <  3:116468/149‑1 (MQ=255)
                      cACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGa                                                                                                                             <  2:71017/149‑1 (MQ=255)
                      cACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGa                                                                                                                             <  1:806182/149‑1 (MQ=255)
                             ccGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAg                                                                                                                      >  4:229647/1‑149 (MQ=255)
                                 aCGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGcccc                                                                                                                  <  1:558412/149‑1 (MQ=255)
                                   gTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCAc                                                                                                                >  2:56691/1‑149 (MQ=255)
                                     aGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGa                                                                                                              >  2:354730/1‑149 (MQ=255)
                                         tGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCg                                                                                                          >  4:481128/1‑149 (MQ=255)
                                             ttACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTa                                                                                                      <  3:210461/149‑1 (MQ=255)
                                             ttACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTCCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTa                                                                                                      >  4:221677/1‑149 (MQ=255)
                                                 aaaaGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTa                                                                                                  <  3:858622/149‑1 (MQ=255)
                                                 aaaaGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTa                                                                                                  >  1:593193/1‑149 (MQ=255)
                                                   aaGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTAcc                                                                                                >  4:26915/1‑149 (MQ=255)
                                                     gCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCaa                                                                                              <  4:371147/149‑1 (MQ=255)
                                                          aaaaaGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGa                                                                                                                             <  3:468788/113‑1 (MQ=255)
                                                          aaaaaGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGa                                                                                                                             >  4:468788/1‑113 (MQ=255)
                                                          aaaaaGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTg                                                                                         >  3:662841/1‑149 (MQ=255)
                                                          aaaaaGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTg                                                                                         >  3:87130/1‑149 (MQ=255)
                                                          aaaaaGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTAGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGtttc                                                                                         >  2:379173/1‑148 (MQ=255)
                                                            aaaGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGcc                                                                                       <  3:229647/149‑2 (MQ=255)
                                                            aaaGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCt                                                                                       <  2:716553/149‑1 (MQ=255)
                                                             aaGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTg                                                                                      <  3:670345/149‑1 (MQ=255)
                                                             aaGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTg                                                                                      <  1:345085/149‑1 (MQ=255)
                                                              aGCAGAGACAAAAAGCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCGGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTgg                                                                                     <  2:793384/149‑1 (MQ=255)
                                                               gCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTgcg                                                                                    >  4:184769/1‑147 (MQ=255)
                                                                cAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGt                                                                                   <  3:802943/149‑1 (MQ=255)
                                                                cAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGt                                                                                   <  1:424032/149‑1 (MQ=255)
                                                                cAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGt                                                                                   <  3:221677/149‑1 (MQ=255)
                                                                      cAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCaaaa                                                                             <  3:317301/149‑1 (MQ=255)
                                                                       aaaaaaCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGCGATGCGAATATTGACCAAGTTTCCTGCGTgaacata                                                                            >  4:6774/1‑143 (MQ=18)
                                                                       aaaaaaCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAAt                                                                            >  3:199887/1‑149 (MQ=255)
                                                                       aaaaaaCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCACGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAAt                                                                            >  2:630882/1‑149 (MQ=255)
                                                                        aaaaaCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAAtt                                                                           <  4:662841/149‑1 (MQ=255)
                                                                         aaaaCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAAttt                                                                          <  1:638377/149‑1 (MQ=255)
                                                                          aaaCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtg                                                                         <  3:858094/149‑1 (MQ=255)
                                                                            aCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggt                                                                       <  3:348111/149‑1 (MQ=255)
                                                                            aCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggt                                                                       <  2:231068/149‑1 (MQ=255)
                                                                            aCCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggt                                                                       <  2:286505/149‑1 (MQ=255)
                                                                             cccccGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggtg                                                                      <  2:593193/149‑1 (MQ=255)
                                                                             cccccGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTtggtg                                                                      <  1:354730/149‑1 (MQ=255)
                                                                               cccGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGa                                                                    >  3:569354/1‑149 (MQ=255)
                                                                               cccGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGa                                                                    <  1:816946/149‑1 (MQ=255)
                                                                                ccacTTTGCAGGGAGGGTGGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGc                                                                                                                     <  3:679894/96‑1 (MQ=255)
                                                                                ccGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGGATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAg                                                                   >  4:603033/1‑149 (MQ=255)
                                                                                ccGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGc                                                                                                                     >  4:679894/1‑99 (MQ=255)
                                                                                ccGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAg                                                                   >  4:292362/1‑149 (MQ=255)
                                                                                ccGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAg                                                                   >  4:151172/1‑149 (MQ=255)
                                                                                  gCTTTGCAGCGTGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCt                                                                 >  3:775027/1‑149 (MQ=255)
                                                                                  gCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTa                                                                                                      >  1:312978/1‑112 (MQ=255)
                                                                                  gCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCt                                                                 >  3:497846/1‑149 (MQ=255)
                                                                                  gCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCACACGATGCGTTTa                                                                                                      <  2:312978/112‑1 (MQ=255)
                                                                                                gTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAAc                                                   <  2:683707/149‑1 (MQ=255)
                                                                                                    gAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCt                                               >  3:25118/1‑149 (MQ=255)
                                                                                                    gAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGAGGTTACTTACCAAGTTTGCTGGGTGCAAAATTTGGTTGAGCTAAGCGGGATCGAACCGCt                                               >  2:6684/1‑149 (MQ=255)
                                                                                                      aaTTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGa                                             <  2:875110/149‑1 (MQ=255)
                                                                                                       aTTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGc                                                <  3:534837/145‑1 (MQ=255)
                                                                                                       aTTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGc                                                >  4:534837/1‑145 (MQ=255)
                                                                                                       aTTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGAc                                            >  1:473870/1‑149 (MQ=255)
                                                                                                        tttGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGAcc                                           >  3:109201/1‑149 (MQ=255)
                                                                                                             tggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCCCTGACCTCTTg                                      >  3:88026/1‑149 (MQ=255)
                                                                                                              ggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGc                                     <  4:109201/149‑1 (MQ=255)
                                                                                                              ggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGc                                     <  4:710904/149‑1 (MQ=255)
                                                                                                                    aaGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgcca                               <  4:119554/149‑1 (MQ=255)
                                                                                                                      gCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGcatgccatg                             <  1:515778/149‑1 (MQ=255)
                                                                                                                           aTCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAga                        <  4:199887/149‑2 (MQ=255)
                                                                                                                           aTCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAga                        <  1:521919/149‑2 (MQ=255)
                                                                                                                           aTCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAgc                        <  3:124902/149‑1 (MQ=255)
                                                                                                                              gAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGct                     >  3:115674/1‑149 (MQ=255)
                                                                                                                               aaCCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCtt                                       >  4:661809/1‑130 (MQ=255)
                                                                                                                               aaCCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCtt                                       <  3:661809/130‑1 (MQ=255)
                                                                                                                               aaCCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGctc                    <  4:775027/149‑1 (MQ=255)
                                                                                                                                    cTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGCGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCa               >  4:193244/1‑149 (MQ=255)
                                                                                                                                         ctctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGa          <  1:700480/149‑1 (MQ=255)
                                                                                                                                         ctctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGa          <  1:856764/149‑1 (MQ=255)
                                                                                                                                          tctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAg         >  3:538613/1‑149 (MQ=255)
                                                                                                                                          gctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAg         >  4:788816/2‑149 (MQ=255)
                                                                                                                                          gctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAg         >  1:197402/2‑149 (MQ=255)
                                                                                                                                                atgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAg   >  2:596750/1‑149 (MQ=255)
                                                                                                                                                 tgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCTCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCACTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGc  <  3:603033/149‑1 (MQ=255)
                                                                                                                                                 tgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGc  >  3:727343/1‑149 (MQ=255)
                                                                                                                                                 tgccacgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGc  <  1:630882/149‑1 (MQ=255)
                                                                                                                                                   ccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCgg                                                           <  1:627769/90‑1 (MQ=255)
                                                                                                                                                   ccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCgg                                                           >  2:627769/1‑90 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAATGATGCTTTCACCAGGCGGCACTAACCCGGACGTAGACTGTATTACAAAAGCGGCAAAAAGCAGAGACAAAAAACCCCCGCTTTGCAGCGAGGGTTGGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGC  >  NC_000913/2517935‑2518228

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