Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,995,769 T→A 5.4% L689Q (CTG→CAG)  xdhA → xanthine dehydrogenase, molybdenum binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,995,7690TA 5.4% 1036.6 / 35.8 300L689Q (CTG→CAG) xdhAxanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base T (152/130);  new base A (4/13);  total (157/143)
Fisher's exact test for biased strand distribution p-value = 2.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 309 total aligned reads displayed.

GCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCGCTGACGC  >  CP009273/2995620‑2995918
                                                                                                                                                     |                                                                                                                                                     
gctgGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACt                                                                                                                                                       <  1:1665296/150‑1 (MQ=255)
     cAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGaaa                                                                                                                                                  >  2:4992053/1‑150 (MQ=255)
        aaGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAgcgc                                                                                                                                              >  1:3755758/1‑150 (MQ=255)
           gTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGt                                                                                                                                            >  2:3019051/1‑150 (MQ=255)
             cAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTtc                                                                                                                                          <  1:2247309/150‑1 (MQ=255)
              aGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTtcg                                                                                                                                         >  1:135313/1‑150 (MQ=255)
                gTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTtcgtc                                                                                                                                       >  1:122037/1‑150 (MQ=255)
                 tACACGGCGGAAGGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTtcgtcg                                                                                                                                      >  2:212156/1‑150 (MQ=255)
                  acacGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGa                                                                                                                                     >  2:4574576/1‑150 (MQ=255)
                      ggcggAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACAGGAAAGCGCGTTCGTCGAAATc                                                                                                                                 <  2:51464/150‑1 (MQ=255)
                      ggcggAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCAACAACTGGAAAGCGCGTTCGTCGAAATc                                                                                                                                 >  2:1686170/1‑150 (MQ=255)
                         ggaatgggaatgGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAAt                                                                                                                              <  2:2256577/150‑1 (MQ=255)
                         ggaatgggaatgGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAAt                                                                                                                              <  2:1823390/150‑1 (MQ=255)
                          gaatgggaatggGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATg                                                                                                                             >  1:40358/1‑150 (MQ=255)
                          gaatgggaatggGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATg                                                                                                                             >  1:3955377/1‑150 (MQ=255)
                           aatgggaatgggCATTGGCTGGGCGCTATTTGAAGAGATGAACATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGa                                                                                                                            <  2:2852761/150‑1 (MQ=255)
                                     ggCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATcc                                                                                                                  <  2:862714/150‑1 (MQ=255)
                                     ggCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATcc                                                                                                                  >  2:1188165/1‑150 (MQ=255)
                                         ttGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCAt                                                                                                              <  2:1389663/150‑1 (MQ=255)
                                         ttGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACAGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCAt                                                                                                              <  2:4720755/150‑1 (MQ=255)
                                           ggCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATAc                                                                                                            >  2:3780245/1‑150 (MQ=255)
                                           ggCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATAc                                                                                                            <  1:1299442/150‑1 (MQ=255)
                                            gCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACg                                                                                                           <  1:662457/150‑1 (MQ=255)
                                                  cgcTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATa                                                                                                     >  1:4048933/1‑150 (MQ=255)
                                                  cgcTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATa                                                                                                     <  1:2466350/150‑1 (MQ=255)
                                                  cgcTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGAGTACAAAATGCCGACCATGCCGCATCTGCCACAACTGGAAACCGCGTTCGTCGAAATCAAACTGACTCAATCCGCATAAGAACATa                                                                                                     >  2:2704358/1‑150 (MQ=255)
                                                   gcTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATaa                                                                                                    <  2:2450756/150‑1 (MQ=255)
                                                     tATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGt                                                                                                  <  2:809222/150‑1 (MQ=255)
                                                       tttGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCa                                                                                                <  1:1271559/150‑1 (MQ=255)
                                                         tGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACt                                                                                              >  2:1006920/1‑150 (MQ=255)
                                                          gAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTg                                                                                             <  2:3585247/150‑1 (MQ=255)
                                                            agagATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTggg                                                                                           <  2:602157/150‑1 (MQ=255)
                                                            agagATGATCATCGATGCGAACAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTggg                                                                                           <  2:4010239/150‑1 (MQ=255)
                                                              agaTGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTgggcg                                                                                         <  2:3521854/150‑3 (MQ=255)
                                                              agaTGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACAGAAGTCACTGGGTg                                                                                         >  2:4176964/1‑150 (MQ=255)
                                                               gatgatCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGa                                                                                        <  2:5213792/150‑1 (MQ=255)
                                                               gatgatCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGa                                                                                        >  1:4196101/1‑150 (MQ=255)
                                                               gatgatCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGa                                                                                        >  2:1272876/1‑150 (MQ=255)
                                                               gatgatCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGa                                                                                        >  2:3889156/1‑150 (MQ=255)
                                                                atgatCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAg                                                                                       >  1:2833298/1‑150 (MQ=255)
                                                                  gatCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGcc                                                                                     <  1:3282194/150‑1 (MQ=255)
                                                                   atcatcGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGcca                                                                                    >  1:2393566/1‑150 (MQ=255)
                                                                     catcGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGccacc                                                                                  >  2:2222760/1‑150 (MQ=255)
                                                                         gATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCaata                                                                              >  1:5161536/1‑150 (MQ=255)
                                                                          aTGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGTGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCaataa                                                                             <  2:4857631/150‑1 (MQ=255)
                                                                          aTGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCaataa                                                                             >  1:3100004/1‑150 (MQ=255)
                                                                            gCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAAtt                                                                           >  1:5130780/1‑150 (MQ=255)
                                                                               aaaaGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCt                                                                        >  1:2486161/1‑150 (MQ=255)
                                                                               aaaaGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCt                                                                        <  1:1450541/150‑1 (MQ=255)
                                                                                 aaGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGt                                                                      <  1:2260469/150‑1 (MQ=255)
                                                                                 aaGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACAGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGt                                                                      <  2:5031045/150‑1 (MQ=255)
                                                                                    cggcgTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGc                                                                   <  1:2304361/150‑1 (MQ=255)
                                                                                        gTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACAGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCgct                                                               <  2:3330363/150‑1 (MQ=255)
                                                                                         tGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCgctg                                                              >  2:2055923/1‑150 (MQ=255)
                                                                                           gTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCgctgct                                                            >  1:1049928/1‑150 (MQ=255)
                                                                                           gTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCgctgct                                                            >  2:4499418/1‑150 (MQ=255)
                                                                                               gTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTc                                                        >  1:5053116/1‑150 (MQ=255)
                                                                                                 aaCCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGt                                                      <  2:908444/150‑1 (MQ=255)
                                                                                                 aaCCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGt                                                      >  1:4742127/1‑150 (MQ=255)
                                                                                                 aaCCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGt                                                      >  1:3374112/1‑150 (MQ=255)
                                                                                                  aCCCCAATCTGCTGGATTACAAAATGCCGTCCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGATCCACCAATAATTCCTGTTGCCGCTGCTATTCGTa                                                     <  2:4082929/150‑1 (MQ=255)
                                                                                                    cccAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAAc                                                   <  2:1846034/150‑1 (MQ=255)
                                                                                                     ccAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTTCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGTTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAAcg                                                  >  2:4956384/1‑150 (MQ=255)
                                                                                                        aTCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCgg                                               <  1:1418906/150‑1 (MQ=255)
                                                                                                            gctgGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGaa                                           <  2:2391013/150‑1 (MQ=255)
                                                                                                              tgGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGa                                         >  2:599403/1‑150 (MQ=255)
                                                                                                               ggATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGAt                                        >  1:3014056/1‑150 (MQ=255)
                                                                                                                gATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATg                                       >  1:9026/1‑150 (MQ=255)
                                                                                                                  ttACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGc                                     <  1:1989927/150‑1 (MQ=255)
                                                                                                                    aCAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTa                                   <  1:783266/150‑1 (MQ=255)
                                                                                                                      aaaaTGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTAcc                                 >  1:2119686/1‑150 (MQ=255)
                                                                                                                      aaaaTGCCGACCATGCAGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTAcc                                 <  2:3143130/150‑1 (MQ=255)
                                                                                                                        aaTGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCgg                               >  1:1928685/1‑150 (MQ=255)
                                                                                                                         aTGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGgt                              <  1:3804606/150‑1 (MQ=255)
                                                                                                                         aTGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGgt                              >  2:4675597/1‑150 (MQ=255)
                                                                                                                          tGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGgtg                             >  1:3722568/1‑150 (MQ=255)
                                                                                                                               aCCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGca                        >  2:1609809/1‑150 (MQ=255)
                                                                                                                               aCCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGca                        >  1:5093537/1‑150 (MQ=255)
                                                                                                                                ccATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTCTTCGTAACGCGGTGAAGATGGCTACCGGTGTTGcaa                       >  2:3106916/1‑150 (MQ=255)
                                                                                                                                  aTGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGcaatc                     <  2:1561553/150‑1 (MQ=255)
                                                                                                                                  aTGCCGGATCTGCCACAACTGGAAAGCGAGTTAGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGTTGAAGATGGCTACCGGTGTTGcaatt                     <  2:311847/150‑2 (MQ=255)
                                                                                                                                   tGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGcaatca                    <  2:3290909/150‑1 (MQ=255)
                                                                                                                                       ggATCTGCCACAACAGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATac                <  1:3494717/150‑1 (MQ=255)
                                                                                                                                        taTCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCTATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATaca               <  1:2700034/149‑1 (MQ=255)
                                                                                                                                        gATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATaca               >  1:3642090/1‑150 (MQ=255)
                                                                                                                                        gATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATaca               <  1:204995/150‑1 (MQ=255)
                                                                                                                                         aTCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATacac              <  1:5024622/150‑1 (MQ=255)
                                                                                                                                           cTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTg            >  2:2967769/1‑150 (MQ=255)
                                                                                                                                             gCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGcc          <  2:3490372/150‑1 (MQ=255)
                                                                                                                                              ccACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCg         <  2:4272792/150‑1 (MQ=255)
                                                                                                                                              ccACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCg         <  1:4600859/150‑1 (MQ=255)
                                                                                                                                              ccACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCg         <  1:1600080/150‑1 (MQ=255)
                                                                                                                                                acaacaGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCGCt       >  2:3137802/1‑150 (MQ=255)
                                                                                                                                                 caACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTTTTGACGCTGCTATTCGTAACGCGGTGAAGATGGCTACAGGTGTTGCAATCAATACACTGCCGCTg      >  2:3672537/1‑150 (MQ=255)
                                                                                                                                                  aaCTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCGCTGa     >  2:4641448/1‑150 (MQ=255)
                                                                                                                                                  aaCTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCGCTGa     >  1:4338765/1‑150 (MQ=255)
                                                                                                                                                   aCTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCGCTGAc    <  2:3421745/150‑1 (MQ=255)
                                                                                                                                                   aCTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCGCTGAc    <  2:1919928/150‑1 (MQ=255)
                                                                                                                                                    cTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCGCTGACg   >  1:1743381/1‑150 (MQ=255)
                                                                                                                                                     tGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCGCTGACGc  >  1:4379574/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTGGCAGAAGGTCAGGTACACGGCGGAATGGGAATGGGCATTGGCTGGGCGCTATTTGAAGAGATGATCATCGATGCTAAAAGCGGCGTGGTCCGTAACCCCAATCTGCTGGATTACAAAATGCCGACCATGCCGGATCTGCCACAACTGGAAAGCGCGTTCGTCGAAATCAATGAGCCGCAATCCGCATACGGACATAAGTCACTGGGTGAGCCACCAATAATTCCTGTTGCCGCTGCTATTCGTAACGCGGTGAAGATGGCTACCGGTGTTGCAATCAATACACTGCCGCTGACGC  >  CP009273/2995620‑2995918

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