Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,995,609 C→A 7.6% L636I (CTT→ATT)  xdhA → xanthine dehydrogenase, molybdenum binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,995,6090CA 7.6% 1074.7 / 40.3 289L636I (CTT→ATT) xdhAxanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base C (139/124);  new base A (10/16);  total (149/140)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 290 total aligned reads displayed.

CGGCGGGCAGCTCTCTGCTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGG  >  CP009273/2995461‑2995756
                                                                                                                                                    |                                                                                                                                                   
  gcggGCAGCTCTCTGCTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCtt                                                                                                                                                   <  1:1322343/149‑1 (MQ=255)
  gcggGCAGCTCTCTGCTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCtt                                                                                                                                                   <  1:2456609/149‑1 (MQ=255)
     ggCAGCTCTCTGCTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAAt                                                                                                                                                <  2:3437606/149‑1 (MQ=255)
          ctctctGCTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCAct                                                                                                                                           <  1:2801344/149‑1 (MQ=255)
            ctctGCTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTTTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCAACAACCCCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCAcagc                                                                                                                                         >  1:1192175/1‑149 (MQ=255)
             tctGCTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGCTGCAGGGCCTAGTCTTCATCAActgct                                                                                                                                        >  2:2721722/1‑149 (MQ=255)
             tctGCTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCActgct                                                                                                                                        >  2:2721844/1‑149 (MQ=255)
             tctGCTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCActgct                                                                                                                                        <  2:2807760/149‑1 (MQ=255)
              ctgctgAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCActgctg                                                                                                                                       <  2:683789/149‑1 (MQ=255)
               tgctgAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTgg                                                                                                                                      <  2:1357792/149‑1 (MQ=255)
                gctgAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGc                                                                                                                                     <  1:2854061/149‑1 (MQ=255)
                  tgAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAg                                                                                                                                   >  2:1072897/1‑149 (MQ=255)
                  tgAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTATTAATACACTGCTGGCAg                                                                                                                                   >  2:1694641/1‑149 (MQ=255)
                      aGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAgg                                                                                                                               <  2:5098965/149‑1 (MQ=255)
                      aGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAgg                                                                                                                               >  2:2410063/1‑149 (MQ=255)
                             tCAAAACCACCACTAACCCACCGGCGTTTGGCTGTAACTTTGTTGATCTGACGGTCGATATTTCGCTGTGCCAAGTCACCATCAACCCCATCCACAACGTTCATGATTCAGGGCATATTATTAATACACAGCTGGCAGAAGGTCAGGTa                                                                                                                        >  1:4533883/1‑149 (MQ=255)
                                  accaccACTAACCCACCGGCGTTTGGCTGTACCTTTGTTTATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCAAAAACCGCATCCTCAACGTTCATGATTCCGGGCATATTCTTAATACACTGCTGGCAGAAGGTCAGGTACAcgg                                                                                                                   >  1:4692716/1‑149 (MQ=255)
                                  accaccACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcgg                                                                                                                   >  2:4191340/1‑149 (MQ=255)
                                   ccaccaCTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggc                                                                                                                  <  1:3827565/149‑1 (MQ=255)
                                    caccacTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcg                                                                                                                 <  1:2992150/149‑1 (MQ=255)
                                    caccacTAACCCACCGGCGTGTGGCTGTACCTTTGTTGATCTGACGGGCGATATTGCGCTGTGCCAAGTCACCATCAACCGCATCCACAACGTTCATGATTCAGGGCCTATTCTTAATCCACGGCTGGCCGAAGGTCAGGTACAcggcg                                                                                                                 >  1:585154/1‑149 (MQ=255)
                                     accacTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg                                                                                                                <  1:4778273/149‑1 (MQ=255)
                                     accacTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACAcggcgg                                                                                                                <  1:3061521/149‑1 (MQ=255)
                                        acTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGCCGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCggaat                                                                                                             >  2:2685185/1‑149 (MQ=255)
                                        acTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCggaat                                                                                                             <  1:2221913/149‑1 (MQ=255)
                                        acTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCggaat                                                                                                             <  1:1863486/149‑1 (MQ=255)
                                        acTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTATTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCggaat                                                                                                             <  2:2933170/149‑1 (MQ=255)
                                                aCCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTTCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATACACTGCTGGCAGAAGGTCAGGTACACGGCGgaatgggaatgg                                                                                                     >  1:4325797/1‑149 (MQ=255)
                                                  cGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTATTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGc                                                                                                   >  1:3571781/1‑149 (MQ=255)
                                                   ggCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGTTTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGTCTGGGAATGGGcc                                                                                                  >  2:5128640/1‑148 (MQ=255)
                                                   ggCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCa                                                                                                  >  1:356745/1‑149 (MQ=255)
                                                   ggCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCa                                                                                                  >  2:394872/1‑149 (MQ=255)
                                                    gCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCAt                                                                                                 <  1:1092162/149‑1 (MQ=255)
                                                    gCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTATTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCAt                                                                                                 <  1:1771526/149‑1 (MQ=255)
                                                       tttGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTgg                                                                                              <  1:2037681/149‑1 (MQ=255)
                                                         tGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCTGAATGGGAATGGGCATTGGCt                                                                                            <  1:4009993/149‑1 (MQ=255)
                                                         tGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCt                                                                                            >  1:2915718/1‑149 (MQ=255)
                                                          ggCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTg                                                                                           <  2:4426074/149‑1 (MQ=255)
                                                          ggCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTg                                                                                           >  2:1204825/1‑149 (MQ=255)
                                                          ggCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTg                                                                                           <  1:4857529/149‑1 (MQ=255)
                                                           gCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTgg                                                                                          >  1:3912093/1‑149 (MQ=255)
                                                           gCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTgg                                                                                          >  2:3027276/1‑149 (MQ=255)
                                                           gCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTgg                                                                                          >  2:2558894/1‑149 (MQ=255)
                                                             tGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGgc                                                                                        <  2:763840/149‑1 (MQ=255)
                                                              gTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGgcg                                                                                       >  1:2506361/1‑149 (MQ=255)
                                                               tACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGgcgc                                                                                      <  1:3202501/149‑1 (MQ=255)
                                                                aCCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCt                                                                                     <  1:3760690/149‑1 (MQ=255)
                                                                  cTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTAt                                                                                   <  2:3716441/149‑1 (MQ=255)
                                                                   tttGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTAtt                                                                                  >  2:1768070/1‑149 (MQ=255)
                                                                     tgttgATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTg                                                                                <  2:199463/149‑1 (MQ=255)
                                                                      gttgATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGa                                                                               >  2:4799666/1‑149 (MQ=255)
                                                                      gttgATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGa                                                                               >  1:982225/1‑149 (MQ=255)
                                                                            cTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGACGCTATTTGAAGAgat                                                                         <  1:1524374/149‑1 (MQ=255)
                                                                            cTGACGGTCGATAGTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAgat                                                                         <  2:4120735/149‑1 (MQ=255)
                                                                             tGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAgatg                                                                        <  2:3498872/149‑1 (MQ=255)
                                                                             tGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTATTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAgatg                                                                        <  1:3686312/149‑1 (MQ=255)
                                                                              gACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAgatga                                                                       >  1:3042125/1‑149 (MQ=255)
                                                                              gACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAgatga                                                                       >  2:1573218/1‑149 (MQ=255)
                                                                                 ggTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGatca                                                                    >  2:1864480/1‑149 (MQ=255)
                                                                                  gTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGatcat                                                                   >  1:4174181/1‑149 (MQ=255)
                                                                                  gTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGatcat                                                                   >  1:3663116/1‑149 (MQ=255)
                                                                                     gATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGa                                                                <  1:1641938/149‑1 (MQ=255)
                                                                                     gATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGAGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGa                                                                >  1:4468584/1‑149 (MQ=255)
                                                                                      atatTGCGCTGTGCTAAGTCAACATCTACCGCATCCTCTACGTTCATGATTCAGTGCATATTCTTAATCCTCTGCTGGGAGAAGGTCAGGTACACGGCGTAATGGGAATGGGCATTTGCTGGGCGCTATTTGAAGAGATGATCATCGAt                                                               <  2:4669878/149‑1 (MQ=255)
                                                                                      atatTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGAt                                                               >  2:283625/1‑149 (MQ=255)
                                                                                        atTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGc                                                             >  1:5001549/1‑149 (MQ=255)
                                                                                        atTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGc                                                             >  1:1453591/1‑149 (MQ=255)
                                                                                        atTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTATTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGc                                                             <  1:4566398/149‑1 (MQ=255)
                                                                                         ttGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCt                                                            <  1:2546317/149‑1 (MQ=255)
                                                                                         ttGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCt                                                            >  1:1369457/1‑149 (MQ=255)
                                                                                           gcgcTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTaa                                                          >  2:5031517/1‑149 (MQ=255)
                                                                                             gcTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTaaaa                                                        <  2:4490751/149‑1 (MQ=255)
                                                                                               tgtgCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAgc                                                      >  2:5011657/1‑149 (MQ=255)
                                                                                                    aaaGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCGTCGATGCTAAAAGCGGCGt                                                 >  2:3831935/1‑149 (MQ=255)
                                                                                                           ccATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGt                                          >  2:2194940/1‑149 (MQ=255)
                                                                                                                 aCCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAAcccc                                    >  2:2358651/1‑149 (MQ=255)
                                                                                                                 aCCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAAcccc                                    >  2:577357/1‑149 (MQ=255)
                                                                                                                   cGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCaa                                  <  1:1044832/149‑1 (MQ=255)
                                                                                                                   cGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCaa                                  <  1:1695708/149‑1 (MQ=255)
                                                                                                                    gCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAAt                                 >  1:2409647/1‑149 (MQ=255)
                                                                                                                     cATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATc                                <  2:4010866/149‑1 (MQ=255)
                                                                                                                     cATCCTCAACGTTCATGATTCAGGGCATATTATTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATc                                <  2:2654261/149‑1 (MQ=255)
                                                                                                                        ccTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATctgc                             >  1:4439987/1‑149 (MQ=255)
                                                                                                                             aCGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGAt                        >  2:1953320/1‑149 (MQ=255)
                                                                                                                              cGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGAtt                       <  2:3203269/149‑1 (MQ=255)
                                                                                                                               gTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTa                      >  1:4158984/1‑149 (MQ=255)
                                                                                                                                ttCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTAc                     >  2:3891838/1‑149 (MQ=255)
                                                                                                                                   aTGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACaaa                  >  1:2057898/1‑149 (MQ=255)
                                                                                                                                    tGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACaaaa                 <  2:1131524/149‑1 (MQ=255)
                                                                                                                                    tGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACaaaa                 >  2:4315177/1‑149 (MQ=255)
                                                                                                                                          aGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGa           >  2:108666/1‑149 (MQ=255)
                                                                                                                                          aGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGa           <  2:2992083/149‑1 (MQ=255)
                                                                                                                                           gggCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGAc          >  1:750050/1‑149 (MQ=255)
                                                                                                                                           gggCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGAc          >  1:3290247/1‑149 (MQ=255)
                                                                                                                                             gCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCa        >  2:1292194/1‑149 (MQ=255)
                                                                                                                                                tatTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGc     <  2:4721794/149‑1 (MQ=255)
                                                                                                                                                tatTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGc     <  2:1420412/149‑1 (MQ=255)
                                                                                                                                                  ttCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCg   >  1:4347727/1‑149 (MQ=255)
                                                                                                                                                  ttCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCg   >  1:5138323/1‑149 (MQ=255)
                                                                                                                                                   tCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCgg  >  2:4064831/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGCGGGCAGCTCTCTGCTGAAAGCTCCATCAAAACCACCACTAACCCACCGGCGTTTGGCTGTACCTTTGTTGATCTGACGGTCGATATTGCGCTGTGCAAAGTCACCATCAACCGCATCCTCAACGTTCATGATTCAGGGCATATTCTTAATCCACTGCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGG  >  CP009273/2995461‑2995756

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