Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 325,740 G→T 57.8% P189P (CCG→CCT betT → choline transporter of high affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273325,7400GT57.8% 38.1 / inf 325P189P (CCG→CCTbetTcholine transporter of high affinity
Reads supporting (aligned to +/- strand):  ref base G (68/69);  new base T (87/101);  total (155/170)
Fisher's exact test for biased strand distribution p-value = 5.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 330 total aligned reads displayed.

GTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGC  >  CP009273/325591‑325889
                                                                                                                                                     |                                                                                                                                                     
gTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCg                                                                                                                                                       >  2:5049373/1‑150 (MQ=255)
    ggACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAg                                                                                                                                                   <  2:1570670/150‑1 (MQ=255)
    ggACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAg                                                                                                                                                   <  2:2812419/150‑1 (MQ=255)
     gACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAAAGGGCCTATAgg                                                                                                                                                  >  2:1811938/1‑150 (MQ=255)
        gCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCa                                                                                                                                               <  2:1923948/150‑1 (MQ=255)
        gCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCa                                                                                                                                               <  1:1074588/150‑1 (MQ=255)
        gCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCa                                                                                                                                               <  2:3430415/150‑1 (MQ=255)
          tGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACt                                                                                                                                             <  1:2281139/150‑1 (MQ=255)
            tttCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCa                                                                                                                                           <  2:3515286/150‑1 (MQ=255)
            tttCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCa                                                                                                                                           <  1:2640511/150‑1 (MQ=255)
               cACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTATCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACTGATTAACGGGCCGATAGGTCACTCAGTg                                                                                                                                        >  2:4057422/1‑150 (MQ=255)
               cACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTATCGGGCCTATAGGTCACTCAGTg                                                                                                                                        <  2:560314/150‑1 (MQ=255)
                    cGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGatat                                                                                                                                   <  1:3060210/150‑1 (MQ=255)
                    cGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATAGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGatat                                                                                                                                   <  1:4280500/150‑1 (MQ=255)
                     ggCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATAtt                                                                                                                                  <  1:4434557/150‑1 (MQ=255)
                         tAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAg                                                                                                                              <  1:1969447/150‑1 (MQ=255)
                         tAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATGTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAg                                                                                                                              <  1:284739/150‑1 (MQ=255)
                          aaCCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGc                                                                                                                             <  1:626623/150‑1 (MQ=255)
                              ggCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTg                                                                                                                         <  1:198786/150‑1 (MQ=255)
                                cTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGAt                                                                                                                       >  2:320477/1‑150 (MQ=255)
                                cTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGAt                                                                                                                       <  2:584551/150‑1 (MQ=255)
                                 tGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATc                                                                                                                      <  2:137585/150‑1 (MQ=255)
                                  ggTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCg                                                                                                                     <  2:2338067/150‑1 (MQ=255)
                                   gTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCgg                                                                                                                    <  2:4763002/150‑1 (MQ=255)
                                   gTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCgg                                                                                                                    <  1:715986/150‑1 (MQ=255)
                                     cGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCa                                                                                                                  >  1:1862313/1‑150 (MQ=255)
                                      gATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATCTTGCAGCGGTGATCGGCAc                                                                                                                 >  1:1104733/1‑150 (MQ=255)
                                       aTGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACt                                                                                                                >  1:4870175/1‑150 (MQ=255)
                                         gTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTAt                                                                                                              <  1:949114/150‑1 (MQ=255)
                                          tATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATc                                                                                                             <  2:4898486/150‑1 (MQ=255)
                                             gcgcTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTc                                                                                                          <  2:2681517/150‑1 (MQ=255)
                                                cggATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGt                                                                                                       <  2:2952803/148‑1 (MQ=255)
                                                  gATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTAt                                                                                                     >  1:3430153/1‑150 (MQ=255)
                                                   aTGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTAtt                                                                                                    <  2:2715023/150‑1 (MQ=255)
                                                      ggCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGcc                                                                                                 <  1:4024413/150‑1 (MQ=255)
                                                         aTGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACt                                                                                              <  1:3162908/150‑1 (MQ=255)
                                                             cgcTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGc                                                                                          >  2:4336971/1‑150 (MQ=255)
                                                              gcTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCt                                                                                         <  1:483764/150‑1 (MQ=255)
                                                              gcTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCt                                                                                         >  2:166681/1‑150 (MQ=255)
                                                               cTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTc                                                                                        <  2:3355857/150‑1 (MQ=255)
                                                               cTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTc                                                                                        <  1:4953156/150‑1 (MQ=255)
                                                               cTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTc                                                                                        <  1:4484755/150‑1 (MQ=255)
                                                                   gATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTa                                                                                    >  1:5074772/1‑150 (MQ=255)
                                                                     tACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATc                                                                                  <  2:1067849/150‑1 (MQ=255)
                                                                       cTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCgg                                                                                >  1:861817/1‑150 (MQ=255)
                                                                           aGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGtgtg                                                                            >  1:1284359/1‑150 (MQ=255)
                                                                           aGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGtgtg                                                                            <  1:1337172/150‑1 (MQ=255)
                                                                               aTCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTCTTGCCACTACGCTCGGTATCGGTGTGGTGc                                                                        >  2:685015/1‑150 (MQ=255)
                                                                                  gTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGc                                                                     >  1:2223265/1‑150 (MQ=255)
                                                                                      taATTTGCCGCTCACCATCCGCTCGGCTCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTaa                                                                 <  2:3282328/150‑1 (MQ=255)
                                                                                       aaTTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAAc                                                                <  2:4279529/150‑1 (MQ=255)
                                                                                       aaTTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAAc                                                                <  1:2804914/150‑1 (MQ=255)
                                                                                         tttGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTa                                                              >  1:5119408/1‑150 (MQ=255)
                                                                                            tccGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTTTGGTGCAGCTTAACTATgg                                                           <  2:3881427/149‑1 (MQ=255)
                                                                                                cTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTg                                                       <  1:4754973/150‑1 (MQ=255)
                                                                                                cTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTg                                                       <  2:2207841/150‑1 (MQ=255)
                                                                                                 tCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGa                                                      >  2:1717446/1‑150 (MQ=255)
                                                                                                 tCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGa                                                      >  1:790205/1‑150 (MQ=255)
                                                                                                 tCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGa                                                      >  2:58024/1‑150 (MQ=255)
                                                                                                  cACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAg                                                     <  2:953247/150‑1 (MQ=255)
                                                                                                   aCCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGc                                                    >  1:613627/1‑150 (MQ=255)
                                                                                                    ccATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCAGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCg                                                   >  1:3261372/1‑150 (MQ=255)
                                                                                                    ccATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCg                                                   >  1:344276/1‑150 (MQ=255)
                                                                                                     cATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGt                                                  >  2:3778195/1‑150 (MQ=255)
                                                                                                        ccGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACt                                               <  1:3477252/150‑1 (MQ=255)
                                                                                                        ccGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACt                                               <  1:4340034/150‑1 (MQ=255)
                                                                                                        ccGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACt                                               <  1:762706/150‑1 (MQ=255)
                                                                                                         cGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTg                                              <  2:4580725/150‑1 (MQ=255)
                                                                                                             cGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTg                                          <  2:4468012/150‑1 (MQ=255)
                                                                                                              ggCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGa                                         <  1:1664615/150‑1 (MQ=255)
                                                                                                                cgcTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGata                                       <  2:4691259/150‑1 (MQ=255)
                                                                                                                  cTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATAtt                                     <  1:1159367/150‑1 (MQ=255)
                                                                                                                  cTGTACCCGATCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATAtt                                     <  1:2403252/150‑1 (MQ=255)
                                                                                                                      aCCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCg                                 >  2:1451828/1‑150 (MQ=255)
                                                                                                                          gATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTc                             <  2:1999442/150‑1 (MQ=255)
                                                                                                                           aTCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCg                            <  1:2769253/150‑1 (MQ=255)
                                                                                                                           aTCTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCg                            <  2:2419558/150‑1 (MQ=255)
                                                                                                                             cTTCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGAt                          <  2:3604203/150‑1 (MQ=255)
                                                                                                                             cTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGAt                          <  1:1480709/150‑1 (MQ=255)
                                                                                                                              ttCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATg                         >  2:5120393/1‑150 (MQ=255)
                                                                                                                               tCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTg                                          >  1:1195382/1‑132 (MQ=255)
                                                                                                                               tCGGTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTg                                          <  2:1195382/132‑1 (MQ=255)
                                                                                                                                 ggTAAACGGATTAACGTGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCTGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGAAGCTTAACTATGGCTAGAGCGTAATTTTTGATATTACCGATTCGAtttcg                      >  1:1266466/1‑146 (MQ=255)
                                                                                                                                 ggTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATggcg                      <  1:1593637/150‑1 (MQ=255)
                                                                                                                                 ggTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATggcg                      <  2:1977270/150‑1 (MQ=255)
                                                                                                                                  gTAAACGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATggcgg                     >  2:660572/1‑150 (MQ=255)
                                                                                                                                   tAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATggcggc                    >  2:3552081/1‑150 (MQ=255)
                                                                                                                                    aaaCGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGgcggcg                   <  1:4936741/150‑1 (MQ=255)
                                                                                                                                    aaaCGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGgcggcg                   >  1:3144160/1‑150 (MQ=255)
                                                                                                                                     aaCGGATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGa                  >  1:4173779/1‑150 (MQ=255)
                                                                                                                                       cGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGaaa                <  2:3990661/150‑1 (MQ=255)
                                                                                                                                         gATTAACGGGCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGc              >  1:556494/1‑150 (MQ=255)
                                                                                                                                         gATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGc              >  1:3921186/1‑150 (MQ=255)
                                                                                                                                         gATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGc              >  2:5216624/1‑150 (MQ=255)
                                                                                                                                          aTTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCg             >  1:1905640/1‑150 (MQ=255)
                                                                                                                                            tAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGc           >  2:1123609/1‑150 (MQ=255)
                                                                                                                                               cGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACt        >  2:2228393/1‑150 (MQ=255)
                                                                                                                                                gggCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTg       <  1:3754016/150‑1 (MQ=255)
                                                                                                                                                  gCCTATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGAt     >  1:5203573/1‑150 (MQ=255)
                                                                                                                                                     gATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGc  >  2:2105277/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCTGGACGCTGTTTCACTACGGCTTAACCGGCTGGTCGATGTATGCGCTGATGGGCATGGCGCTCGGATACTTTAGCTATCGTTATAATTTGCCGCTCACCATCCGCTCGGCGCTGTACCCGATCTTCGGTAAACGGATTAACGGGCCGATAGGTCACTCAGTGGATATTGCAGCGGTGATCGGCACTATCTTCGGTATTGCCACTACGCTCGGTATCGGTGTGGTGCAGCTTAACTATGGCTTGAGCGTACTGTTTGATATTCCCGATTCGATGGCGGCGAAAGCGGCACTGATCGC  >  CP009273/325591‑325889

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: