Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,576,325 C→A 77.9% E171D (GAG→GAT nudK ← GDP‑mannose pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,576,3250CA77.9% 567.2 / 207.5 288E171D (GAG→GATnudKGDP‑mannose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (43/21);  new base A (121/103);  total (164/124)
Fisher's exact test for biased strand distribution p-value = 6.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 295 total aligned reads displayed.

GTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCT  >  CP009273/2576178‑2576472
                                                                                                                                                   |                                                                                                                                                   
gTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTAtctc                                                                                                                                                   >  1:1107434/1‑150 (MQ=255)
gTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTAtatc                                                                                                                                                   <  2:3142452/150‑1 (MQ=255)
 tATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCg                                                                                                                                                  >  2:4473499/1‑150 (MQ=255)
   tCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGcc                                                                                                                                                >  1:50317/1‑150 (MQ=255)
    cTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCg                                                                                                                                               >  1:2715972/1‑150 (MQ=255)
        gTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGttt                                                                                                                                           >  2:3695094/1‑150 (MQ=255)
          tGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTATTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTg                                                                                                                                         >  1:1007542/1‑150 (MQ=255)
           gcgcAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGa                                                                                                                                        <  2:640506/150‑1 (MQ=255)
            cgcAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGat                                                                                                                                       >  1:5105848/1‑150 (MQ=255)
             gcAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGatc                                                                                                                                      <  1:4500554/150‑1 (MQ=255)
              cAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGatca                                                                                                                                     <  2:3763481/150‑1 (MQ=255)
                 gTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCAtct                                                                                                                                  <  1:2071736/150‑1 (MQ=255)
                 gTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCAtct                                                                                                                                  <  1:2470595/150‑1 (MQ=255)
                    gcgcGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCa                                                                                                                               >  2:1994194/1‑150 (MQ=255)
                         gCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGcc                                                                                                                          >  1:3855840/1‑150 (MQ=255)
                             aaTCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGc                                                                                                                      <  2:4915981/150‑1 (MQ=255)
                                   aTTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACg                                                                                                                <  1:317234/150‑1 (MQ=255)
                                   aTTTGTCGGATTGAGCCAATAATAAACAACCAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACg                                                                                                                >  2:1175033/1‑150 (MQ=255)
                                    tttGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGACCgg                                                                                                               >  2:5246368/1‑150 (MQ=255)
                                      tGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCa                                                                                                             >  1:3600111/1‑150 (MQ=255)
                                        tCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGc                                                                                                           <  1:2139499/150‑1 (MQ=255)
                                         cGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCt                                                                                                          <  2:5052268/150‑1 (MQ=255)
                                            aTTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGa                                                                                                       >  1:5472913/1‑150 (MQ=255)
                                             ttGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAg                                                                                                      >  1:5164497/1‑150 (MQ=255)
                                                aGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCAc                                                                                                   >  1:4046356/1‑150 (MQ=255)
                                                   cAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTc                                                                                                >  2:2011512/1‑150 (MQ=255)
                                                     ataataAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCaa                                                                                              <  1:4473499/150‑1 (MQ=255)
                                                     ataataAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCaa                                                                                              >  1:722976/1‑150 (MQ=255)
                                                       aataaACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTTTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAata                                                                                            >  2:4228165/1‑150 (MQ=255)
                                                       aataaACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAata                                                                                            >  1:5174687/1‑150 (MQ=255)
                                                       aataaACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAata                                                                                            >  2:2708467/1‑150 (MQ=255)
                                                        ataaACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAatat                                                                                           >  1:5603700/1‑150 (MQ=255)
                                                           aacaacGAAATATTATTGTTATATTTCAGTCCATTAAATGTGACGTTTGCAAAGAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatctt                                                                                        <  1:3247855/150‑1 (MQ=255)
                                                              aacGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcat                                                                                     >  1:4760547/1‑150 (MQ=255)
                                                              aacGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcat                                                                                     >  1:4091894/1‑150 (MQ=255)
                                                              aacGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcat                                                                                     >  2:2075173/1‑150 (MQ=255)
                                                               acGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcatc                                                                                    >  1:4409049/1‑150 (MQ=255)
                                                                 gAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcatctt                                                                                  >  1:4113641/1‑150 (MQ=255)
                                                                   aaTATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCg                                                                                >  2:916176/1‑150 (MQ=255)
                                                                   aaTATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCg                                                                                >  2:75199/1‑150 (MQ=255)
                                                                     tattatTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGac                                                                              <  1:258491/150‑1 (MQ=255)
                                                                       ttattGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGacac                                                                            <  2:334754/150‑1 (MQ=255)
                                                                       ttattGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGacac                                                                            >  1:272340/1‑150 (MQ=255)
                                                                         attGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCg                                                                          <  2:1452591/150‑1 (MQ=255)
                                                                           tGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGcc                                                                        >  2:1101822/1‑150 (MQ=255)
                                                                             ttATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGTCACCGCCAc                                                                      >  2:5510057/1‑150 (MQ=255)
                                                                             ttATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCAc                                                                      <  2:2945231/150‑1 (MQ=255)
                                                                              tATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCAcc                                                                     >  2:3994863/1‑150 (MQ=255)
                                                                              tATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCAcc                                                                     <  2:1738269/150‑1 (MQ=255)
                                                                              tATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCAcc                                                                     <  1:2363001/150‑1 (MQ=255)
                                                                               aTTTTTCAGTCCATGAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCg                                                                    <  2:4886937/150‑1 (MQ=255)
                                                                                    tCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGtt                                                               >  1:5328920/1‑150 (MQ=255)
                                                                                      aGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAg                                                             >  2:984347/1‑150 (MQ=255)
                                                                                         ccATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAgcgc                                                          >  1:1006737/1‑150 (MQ=255)
                                                                                         ccATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAgcgc                                                          >  1:2924961/1‑150 (MQ=255)
                                                                                         ccATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAgcgc                                                          >  1:1184434/1‑150 (MQ=255)
                                                                                          cATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAgcgcg                                                         >  1:2797143/1‑150 (MQ=255)
                                                                                            ttAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCt                                                       <  1:133176/150‑1 (MQ=255)
                                                                                             tAATTGTGAAGTTTGCAAATAGTTAAGCAATAACTCCGTCTTACCGTCTCGTATATCGCCGGTTTTGATCATCTCCAGCTCCTGGCTGATCGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTg                                                      <  2:4815301/150‑1 (MQ=255)
                                                                                             tAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTg                                                      <  2:1471429/150‑1 (MQ=255)
                                                                                              agaTGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGa                                                     >  1:4755566/3‑150 (MQ=255)
                                                                                              aaaTGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGa                                                     >  1:2000174/1‑150 (MQ=255)
                                                                                                aTGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGAtt                                                   <  2:4699494/150‑1 (MQ=255)
                                                                                                   tgACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGAGTGTc                                                >  1:3437085/1‑150 (MQ=255)
                                                                                                    gACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCa                                               >  2:3257645/1‑150 (MQ=255)
                                                                                                           gCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTAtt                                        <  2:2469982/150‑1 (MQ=255)
                                                                                                           gCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTAtt                                        >  2:5659777/1‑150 (MQ=255)
                                                                                                            cAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTc                                       <  1:555773/150‑1 (MQ=255)
                                                                                                            cAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTc                                       <  1:4887544/150‑1 (MQ=255)
                                                                                                            cAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTc                                       <  2:1604419/150‑1 (MQ=255)
                                                                                                             aaaTAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCg                                      >  2:357346/1‑150 (MQ=255)
                                                                                                             aaaTAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCg                                      >  1:86321/1‑150 (MQ=255)
                                                                                                             aaaTAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCg                                      >  1:616019/1‑150 (MQ=255)
                                                                                                               aTAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGc                                    <  1:3535475/150‑1 (MQ=255)
                                                                                                                 aGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGa                                  >  2:5571387/1‑150 (MQ=255)
                                                                                                                    tAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaa                               >  2:2655442/1‑150 (MQ=255)
                                                                                                                    tAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaa                               <  1:419220/150‑1 (MQ=255)
                                                                                                                    tAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaa                               >  1:2735710/1‑150 (MQ=255)
                                                                                                                    tAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaa                               >  2:2158094/1‑150 (MQ=255)
                                                                                                                     aaGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaaa                              >  1:1570875/1‑150 (MQ=255)
                                                                                                                     aaGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaaa                              <  1:3695094/150‑1 (MQ=255)
                                                                                                                       gCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaaaaa                            <  2:4936733/150‑1 (MQ=255)
                                                                                                                         aaTAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAg                          >  1:306813/1‑150 (MQ=255)
                                                                                                                         aaTAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAg                          >  1:4691351/1‑150 (MQ=255)
                                                                                                                          aTAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGt                         >  1:5272926/1‑150 (MQ=255)
                                                                                                                          aTAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGt                         >  1:2455022/1‑150 (MQ=255)
                                                                                                                           tAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTg                        <  1:1994194/150‑1 (MQ=255)
                                                                                                                                  gTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGc                 <  1:1704848/150‑1 (MQ=255)
                                                                                                                                   tCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCt                >  1:2772080/1‑150 (MQ=255)
                                                                                                                                     ttACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCa              >  1:2118217/1‑150 (MQ=255)
                                                                                                                                     ttACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCa              >  2:2589359/1‑150 (MQ=255)
                                                                                                                                        ccGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTc           >  2:4828917/1‑150 (MQ=255)
                                                                                                                                          gTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTcac         <  2:5164497/150‑1 (MQ=255)
                                                                                                                                          gTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTcac         >  2:4150029/1‑150 (MQ=255)
                                                                                                                                          gTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGATCGAGCAATTCAATATCTTAATCTTAGACACCTCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTcac         <  2:4432743/150‑1 (MQ=255)
                                                                                                                                           tCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTcaca        <  2:4074839/150‑1 (MQ=255)
                                                                                                                                            cACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCacac       <  1:829793/150‑1 (MQ=255)
                                                                                                                                            cACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCacac       >  1:4286116/1‑150 (MQ=255)
                                                                                                                                             aCGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACAcc      <  2:1032134/150‑1 (MQ=255)
                                                                                                                                                tataTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGcc   >  1:2659128/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCCT  >  CP009273/2576178‑2576472

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: