Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,912,734 T→C 6.4% Q4R (CAA→CGA)  gudD ← (D)‑glucarate dehydratase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,912,7340TC 6.4% 244.8 / 4.1 89Q4R (CAA→CGA) gudD(D)‑glucarate dehydratase 1
Reads supporting (aligned to +/- strand):  ref base T (39/44);  new base C (3/3);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

CCTACGCCAGTGTGACCAGAATTATCTTTGATAATCACAATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATG  >  CP009273/2912587‑2912878
                                                                                                                                                   |                                                                                                                                                
ccTACGCCAGTGTGACCAGAATTATCTTTGATAATCACAATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATGACCTGCATTTCAGTAACAACAGGCGTCGTAAATTg                                                                                                                                                 >  1:1124562/1‑149 (MQ=255)
  tACGCCAGTGTGACCAGAATTATCTTTGATAATCACAATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACACCATGCGGCGTAACTTgag                                                                                                                                               >  2:100986/1‑149 (MQ=255)
  tACGCCAGTGTGACCAGAATTATCTTTGATAATCACAATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATCTCAGTAACAACACGCACCGTAACTCgag                                                                                                                                               >  2:893633/1‑149 (MQ=255)
   aCGCCAGTGTGACCAGAATTATCTTTGATAATCACAATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGTCCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTgaga                                                                                                                                              <  1:683965/149‑1 (MQ=255)
    cGCCAGTGTGACCAGAATTATCTTTGATAATCACAATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGaa                                                                                                                                             >  1:783877/1‑149 (MQ=255)
    cGCCAGTGTGACCAGAATTATCTTGGATAATCACAATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGaa                                                                                                                                             <  1:490292/149‑1 (MQ=255)
       cAGTGTGACCAGAATTATCTTTGATAATCACAATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTgaacaccc                                                                                                                                          >  2:643289/1‑143 (MQ=255)
          tgtgACCAGAATTATCTTTGATAATCACAATATTACGCGTCAAGAACGGTGCGTGTGCACCACTCAGCTTCATCCGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCCGTAACAACAGGCGTCGTAAATTGCGAACGCata                                                                                                                                       >  1:306692/1‑149 (MQ=255)
              aCCAGAATTATCTTTGATAATCACAATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTa                                                                                                                                   >  1:1194811/1‑149 (MQ=255)
               ccAGAATTATCTTTGATAATCACAATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTAt                                                                                                                                  <  1:211097/149‑1 (MQ=255)
                    aTTATCTTGGAAAATCCCAATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCt                                                                                                                             <  2:160875/149‑1 (MQ=255)
                     ttATCTTTGATAAACACAATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGGACTCATATTTATGTCCtt                                                                                                                            >  2:859628/1‑149 (MQ=255)
                         cTTTGATAATCAAAATATTACGCGTAAAGAACTGTGCGTGTGCACCCCTCAGATTCATCAGCATACTGTCATGTTCCGTCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTc                                                                                                                        <  1:793935/149‑1 (MQ=255)
                          tttGATAATCACAACATTACGCGTAAAGGACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCa                                                                                                                       <  1:933393/149‑1 (MQ=255)
                            tGATATTCACAATTTTACGCGTAGGGAGGGGGGCGTGTGCACCACTTAGATTCATGAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGa                                                                                                                     <  2:1122198/149‑1 (MQ=255)
                             gATAATCACAATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGaa                                                                                                                    >  1:1080571/1‑149 (MQ=255)
                              aTAATCACAATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAAt                                                                                                                   <  1:89551/149‑1 (MQ=255)
                                  tCACAATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAAt                                                                                                               <  1:1116464/149‑1 (MQ=255)
                                   cacaATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATg                                                                                                              >  1:1162818/1‑149 (MQ=255)
                                   cacaATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTAATATTTATGTCCTTTTTCCTAATctttt                                                                                                              >  2:808889/1‑145 (MQ=255)
                                     caATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATCGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGAc                                                                                                            >  1:1116470/1‑149 (MQ=255)
                                     caATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATCGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGAc                                                                                                            >  1:501807/1‑149 (MQ=255)
                                       atatTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACgg                                                                                                          >  1:517502/1‑149 (MQ=255)
                                       atatTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGGAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACgg                                                                                                          <  1:1032566/149‑1 (MQ=255)
                                              gcgTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCAGGACCCGCCAGAGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATCTTTATGTCCTTTTTCAGAATCAATGACGGCCGaaaa                                                                                                   >  1:294885/1‑149 (MQ=255)
                                                gTAAACAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACg                                                                                                 <  1:121191/149‑1 (MQ=255)
                                                  aaaGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACggg                                                                                               <  2:1174237/149‑1 (MQ=255)
                                                   aaGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCACCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTAAGAATAAATGAAGAACGAAAACGGGa                                                                                              >  2:426540/1‑149 (MQ=255)
                                                        cGGGGCGTCTGCACCCCTCAGATTCATCAGCATACTGCCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGATCTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGttt                                                                                         <  2:1124562/149‑1 (MQ=255)
                                                         ggTTCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCAGTTCCGTAACAACAGGCGTCGTAAATTGAGGACTAATATTTAGATCCGTTTTGAGAATCAATGACGGCCGAAACCGGGCCGtata                                                                                        >  2:1101616/1‑149 (MQ=255)
                                                         ggTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGGAAACGGGACGGTTa                                                                                        >  1:507728/1‑149 (MQ=255)
                                                             cGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCAGCCCCCGCGTGAGGCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGt                                                                                    <  1:35512/149‑1 (MQ=255)
                                                               tgtgCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGCCa                                                                                  >  1:978647/1‑149 (MQ=255)
                                                               tgtgCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTGAAGAAGCAACAAGCGTCGTAAATTGAGACCGCAAATTGAGGTCCTTTTTCAGAAACAATGCCGGCCGGCGACGGGTCGTTTGCGGtga                                                                                  >  1:234859/1‑147 (MQ=255)
                                                                 tgCACCGCTCAGATTCATCAGCATACTGTGATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCaaa                                                                                <  2:56178/149‑1 (MQ=255)
                                                                 tgCACCCCTCAGGTTCATCAGCATCCTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCaaa                                                                                <  2:671724/149‑1 (MQ=255)
                                                                 tgCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCaaa                                                                                <  1:827258/149‑1 (MQ=255)
                                                                 tgCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCaaa                                                                                <  1:291466/149‑1 (MQ=255)
                                                                 tgCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATCGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCaaa                                                                                <  2:501807/149‑1 (MQ=255)
                                                                 tgCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATCGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCaaa                                                                                <  2:1116470/149‑1 (MQ=255)
                                                                    accccccAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGt                                                                             <  2:72662/143‑1 (MQ=255)
                                                                            gATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCgg                                                                     <  2:936023/149‑1 (MQ=255)
                                                                            gATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCAGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCgg                                                                     <  2:693133/149‑1 (MQ=255)
                                                                               tcatcaGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGAGACGTTTACGGTCGAAGGTACACCCAGGGa                                                                  >  2:1113331/1‑149 (MQ=255)
                                                                               tcatcaGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACACATATTTATGTCCTTTTTCAGAATCAAGGGCAGCCGAAAACGGGCCGTTCTCGGTCACAGGTCCAGCGGGGGa                                                                  >  2:36554/1‑149 (MQ=255)
                                                                                catcaGCGTACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCATAATCAATGACGGCCGAAAACGGGACGGTTACGGTCAAAGGTCCGGCCggggct                                                                 >  2:608987/1‑147 (MQ=255)
                                                                                catcaGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGAt                                                                 <  1:399075/149‑1 (MQ=255)
                                                                                catcaGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAAACAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGCTCCAGCCGGGGAt                                                                 <  1:849610/149‑1 (MQ=255)
                                                                                 ttctgcattttGTCATGACCCGCTCACGGGATAACCTGCATTTCAGTAACAACGGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATc                                                                <  1:92594/139‑1 (MQ=255)
                                                                                 atcaGCATACTGTCCTGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCCTCGCAGATTGTGAACTCATATTGATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATc                                                                <  2:443548/149‑1 (MQ=255)
                                                                                 atcaGCATACTGTCATGACCGGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATc                                                                <  2:840705/149‑1 (MQ=255)
                                                                                 atcaGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATc                                                                <  2:753632/149‑1 (MQ=255)
                                                                                 atcaGCATACTGCCATGACCCGCCCCCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATCGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAACCGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATc                                                                <  2:841288/149‑1 (MQ=255)
                                                                                 ataacgttACTGTCATTACCCGCCACCGGGATAACATCCATTTCAGTAACGACAGGCGGCGGAAATTGAGAACTCATATTTATGTTCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATc                                                                <  2:956300/142‑1 (MQ=255)
                                                                                  tcaGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGGATCAATGACGGCCGAAGACGGGACCTTTACAGTCAAAGGTACCGGCGGGGacag                                                               >  2:424093/1‑146 (MQ=255)
                                                                                     gCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGt                                                            >  1:541850/1‑149 (MQ=255)
                                                                                      cATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTGCAGAAACGATGACGGCCGAAAACGGGACGGTTACGGTCAAATGTCCAGCCGGGGGTCCGGtc                                                           >  2:32853/1‑148 (MQ=255)
                                                                                      cATACTGTCATGACCCGCCACCGAGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTACAGCCGGGGATCAGGTa                                                           >  1:301533/1‑149 (MQ=255)
                                                                                        tACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATTTACTTTTTCAGAATCAATGACGGCCGAAAGCGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACt                                                         <  1:1160317/149‑1 (MQ=255)
                                                                                        tACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGCAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGACTCGGGTAcc                                                         >  1:858749/1‑148 (MQ=255)
                                                                                               aTGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCCTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCCCCCGGGGATCAGGTACTGCATCgg                                                  >  2:511418/1‑149 (MQ=255)
                                                                                                   cccGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTCTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCt                                              >  1:1192276/1‑149 (MQ=255)
                                                                                                    ccGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTg                                             >  2:741449/1‑149 (MQ=255)
                                                                                                     cGCCACCGTGATATCCTCTATTTCAGTAACAACAGTCGTCTTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATCACGGCCGAAACCTGGACTTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGc                                            <  2:623745/149‑1 (MQ=255)
                                                                                                            gggATCACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTCCTGCATCGGACCTGCCTCGTTa                                     >  2:830942/1‑149 (MQ=255)
                                                                                                            gggATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTa                                     <  1:741449/149‑1 (MQ=255)
                                                                                                             ggATAACCTGCATTTCAGTAACAACAGGCGTCGTAATTTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGCTCAGGTACTGCATCGGAGCTGCGTCGTTAc                                    >  1:728522/1‑149 (MQ=255)
                                                                                                               aTAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGAGCAGGTACTGCATGGGACCTGCGTCGTTAcgc                                  >  1:10939/1‑149 (MQ=255)
                                                                                                               aTAACCTGCATTTCAGTAACAACAGGCGTCGCAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTAcgc                                  <  1:831658/149‑1 (MQ=255)
                                                                                                               aTAACCTGCATTTCAGGAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTAcgc                                  <  2:302527/149‑1 (MQ=255)
                                                                                                                    cTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTAGTTACGCGCAcc                             >  1:1078991/1‑149 (MQ=255)
                                                                                                                      gCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGc                           <  1:1049056/149‑1 (MQ=255)
                                                                                                                       cATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGcc                          <  2:557675/149‑1 (MQ=255)
                                                                                                                        aTTTCATTAACAACAGGCGTCGTAAATTAAGAACTCAGATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAATCGGGACGTTTACGGTTAAAGGTCCAGCCGGGGATCAGGTATTGCATCGGACCTGCGTCGTTACGCGCACCGCCa                         <  2:387999/149‑1 (MQ=255)
                                                                                                                            cAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCa                     <  2:1104063/149‑1 (MQ=255)
                                                                                                                             aGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAg                    <  2:1001090/149‑1 (MQ=255)
                                                                                                                               tAACAGCAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGAc                  >  2:526078/1‑149 (MQ=255)
                                                                                                                               tAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGAc                  >  2:472309/1‑149 (MQ=255)
                                                                                                                               tAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGAc                  >  2:378405/1‑149 (MQ=255)
                                                                                                                                    acaGGCGTCGTAAATTTAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTta             <  1:984266/149‑1 (MQ=255)
                                                                                                                                    acaGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTta             <  1:1199492/149‑1 (MQ=255)
                                                                                                                                       ggCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAg          <  2:836434/149‑1 (MQ=255)
                                                                                                                                       ggCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAg          <  1:720731/149‑1 (MQ=255)
                                                                                                                                          gtcgtAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCACTGAAGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACGTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATACGcc       >  2:226892/1‑149 (MQ=255)
                                                                                                                                          gtcgtAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGcc       >  2:1041375/1‑149 (MQ=255)
                                                                                                                                            cgtAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTc     <  2:783877/149‑1 (MQ=255)
                                                                                                                                             gtAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCa    >  1:400894/1‑149 (MQ=255)
                                                                                                                                              tAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCAt   <  1:175106/149‑1 (MQ=255)
                                                                                                                                               aaaTTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGCCGTTTATAGGCCTCCtg  >  1:60546/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCTACGCCAGTGTGACCAGAATTATCTTTGATAATCACAATATTACGCGTAAAGAACGGTGCGTGTGCACCACTCAGATTCATCAGCATACTGTCATGACCCGCCACCGGGATAACCTGCATTTCAGTAACAACAGGCGTCGTAAATTGAGAACTCATATTTATGTCCTTTTTCAGAATCAATGACGGCCGAAAACGGGACGTTTACGGTCAAAGGTCCAGCCGGGGATCAGGTACTGCATCGGACCTGCGTCGTTACGCGCACCGCCAGGCAGACGTTTATAGGCCTCATG  >  CP009273/2912587‑2912878

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: