Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,576,325 C→A 92.6% 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,3250CA92.6% 778.7 / 45.0 249E171D (GAG→GATnudKGDP‑mannose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (16/3);  new base A (115/115);  total (131/118)
Fisher's exact test for biased strand distribution p-value = 3.94e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

GTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCC  >  CP009273/2576178‑2576471
                                                                                                                                                   |                                                                                                                                                  
gTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTAtatc                                                                                                                                                  <  1:2806330/150‑1 (MQ=255)
gTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTAtatc                                                                                                                                                  >  1:3050551/1‑150 (MQ=255)
  aTCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTTTTATTTTTCAGTCAAAAAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGc                                                                                                                                                >  1:459150/1‑150 (MQ=255)
    cTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCg                                                                                                                                              <  2:791364/150‑1 (MQ=255)
       cGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGtt                                                                                                                                           >  2:1527558/1‑150 (MQ=255)
           gcgcAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGa                                                                                                                                       <  2:2077251/150‑1 (MQ=255)
              cAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGatca                                                                                                                                    >  2:2243556/1‑150 (MQ=255)
               aGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGatcat                                                                                                                                   >  2:4068008/1‑150 (MQ=255)
                ggTAGCGCGGCTCAATCTGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGatcatc                                                                                                                                  <  1:3670447/150‑1 (MQ=255)
                ggTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGatcatc                                                                                                                                  <  1:704822/150‑1 (MQ=255)
                    gcgcGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCa                                                                                                                              <  2:2511781/150‑1 (MQ=255)
                      gcgGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAgc                                                                                                                            <  1:1554784/150‑1 (MQ=255)
                                  gATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAAc                                                                                                                <  2:3182927/150‑1 (MQ=255)
                                  gATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACAGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAAc                                                                                                                >  2:3301533/1‑150 (MQ=255)
                                   aTTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACg                                                                                                               <  1:3879889/150‑1 (MQ=255)
                                      tGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCa                                                                                                            <  1:3384213/150‑1 (MQ=255)
                                      tGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCa                                                                                                            <  2:1609073/150‑1 (MQ=255)
                                        tCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGc                                                                                                          >  2:2014521/1‑150 (MQ=255)
                                        tCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGc                                                                                                          >  2:542014/1‑150 (MQ=255)
                                          ggATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTc                                                                                                        <  2:2114029/150‑1 (MQ=255)
                                           gATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCg                                                                                                       >  1:2722676/1‑150 (MQ=255)
                                              tGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGc                                                                                                    >  2:4344518/1‑150 (MQ=255)
                                                aGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAACTGTGACGTTTGCAAAGAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCAc                                                                                                  <  1:2609130/150‑1 (MQ=255)
                                                aGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGACCGGCAGCTCGAGCAc                                                                                                  >  2:2755971/1‑150 (MQ=255)
                                                aGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCAc                                                                                                  <  1:1028743/150‑1 (MQ=255)
                                                aGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCAc                                                                                                  <  2:3118147/150‑1 (MQ=255)
                                                  ccAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACtt                                                                                                <  1:947886/150‑1 (MQ=255)
                                                       aataaACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAata                                                                                           >  2:1787461/1‑150 (MQ=255)
                                                          aaaCAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATGGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatct                                                                                        <  1:2128297/150‑1 (MQ=255)
                                                            acaaaGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACAACGTCTAACCGTCACGTATATAGCCGGTTTTGATAATCTCAAGCGCCTGGCTGAACGGCAGCTAGAGCACTTCAATatcttc                                                                                      >  2:1703742/1‑150 (MQ=255)
                                                             caacGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttca                                                                                     >  1:795392/1‑150 (MQ=255)
                                                              aacGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcat                                                                                    >  2:315161/1‑150 (MQ=255)
                                                               acGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcatc                                                                                   >  1:4302308/1‑150 (MQ=255)
                                                               acGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcatc                                                                                   <  2:3792806/150‑1 (MQ=255)
                                                               acGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcatc                                                                                   <  2:4086091/150‑1 (MQ=255)
                                                                  aaaTATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcatcttc                                                                                <  1:3878966/150‑1 (MQ=255)
                                                                   aaTATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCg                                                                               >  1:3746897/1‑150 (MQ=255)
                                                                   aaTATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCg                                                                               >  2:645431/1‑150 (MQ=255)
                                                                    atatTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGa                                                                              >  1:2613383/1‑150 (MQ=255)
                                                                      attattGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGaca                                                                            <  1:137669/150‑1 (MQ=255)
                                                                         attGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCCCTTCAATATCTTCATCTTCGACACCg                                                                         >  2:3081817/1‑150 (MQ=255)
                                                                         attGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCg                                                                         >  1:3111474/1‑150 (MQ=255)
                                                                         attGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCACCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCg                                                                         <  2:923180/150‑1 (MQ=255)
                                                                           tGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGcc                                                                       >  2:1918827/1‑150 (MQ=255)
                                                                           tGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGcc                                                                       >  1:2167375/1‑150 (MQ=255)
                                                                           tGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGcc                                                                       >  2:2327835/1‑150 (MQ=255)
                                                                            gTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCa                                                                      >  2:2173266/1‑150 (MQ=255)
                                                                             ttATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCAc                                                                     >  1:3211928/1‑150 (MQ=255)
                                                                              tATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCAcc                                                                    <  2:245704/150‑1 (MQ=255)
                                                                                tttttCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCgg                                                                  >  2:2734707/1‑150 (MQ=255)
                                                                                   ttCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGt                                                               >  1:3609931/1‑150 (MQ=255)
                                                                                     cAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTa                                                             <  2:1931838/150‑1 (MQ=255)
                                                                                            ttAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCt                                                      <  2:204429/150‑1 (MQ=255)
                                                                                            ttAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCACCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCt                                                      >  2:1179311/1‑150 (MQ=255)
                                                                                               aaTGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGAt                                                   >  1:351427/1‑150 (MQ=255)
                                                                                                aTGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGAtt                                                  <  1:1527555/150‑1 (MQ=255)
                                                                                                aTGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGAtt                                                  <  1:3293163/150‑1 (MQ=255)
                                                                                                aTGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGAtt                                                  <  1:1098303/150‑1 (MQ=255)
                                                                                                  gtgACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTAGACACCGCCACCGGCGTTAGCGCGCTGATTGt                                                >  1:4076613/1‑150 (MQ=255)
                                                                                                    gACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCa                                              <  1:1025181/150‑1 (MQ=255)
                                                                                                    gACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCAGATTGTCa                                              >  1:3146313/1‑150 (MQ=255)
                                                                                                     aCGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCAc                                             >  2:4152693/1‑150 (MQ=255)
                                                                                                     aCGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCAc                                             <  2:4089641/150‑1 (MQ=255)
                                                                                                       gTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTg                                           >  1:4231613/1‑150 (MQ=255)
                                                                                                       gTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTg                                           >  1:3932989/1‑150 (MQ=255)
                                                                                                        tttGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGt                                          <  1:2835485/150‑1 (MQ=255)
                                                                                                            cAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCCGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTc                                      >  2:827656/1‑150 (MQ=255)
                                                                                                            cAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTc                                      <  1:2246696/150‑1 (MQ=255)
                                                                                                             aaaTAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCg                                     <  2:413256/150‑1 (MQ=255)
                                                                                                              aaTAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCgg                                    >  2:359998/1‑150 (MQ=255)
                                                                                                               aTAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGc                                   >  1:38016/1‑150 (MQ=255)
                                                                                                                tAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGca                                  <  1:3464254/150‑2 (MQ=255)
                                                                                                                tAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGca                                  <  1:1940285/150‑2 (MQ=255)
                                                                                                                  gTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGAt                                <  2:4227931/150‑1 (MQ=255)
                                                                                                                  gTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGAt                                >  1:2890922/1‑150 (MQ=255)
                                                                                                                    tAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaa                              >  1:3376536/1‑150 (MQ=255)
                                                                                                                    tAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaa                              <  2:594189/150‑1 (MQ=255)
                                                                                                                    tAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaa                              >  1:1180354/1‑150 (MQ=255)
                                                                                                                      aGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaaaa                            >  1:2483910/1‑150 (MQ=255)
                                                                                                                      aGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaaaa                            >  2:526504/1‑150 (MQ=255)
                                                                                                                       gCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCCGCTCGAGCACGTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCGGATTGTCACAGTATTCGGCGATaaaaa                           >  1:3732950/1‑150 (MQ=255)
                                                                                                                        cAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaaaaaa                          <  1:825265/150‑1 (MQ=255)
                                                                                                                         aaTAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAg                         >  1:1735408/1‑150 (MQ=255)
                                                                                                                            aaCACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTgg                      <  1:169188/150‑1 (MQ=255)
                                                                                                                             acacCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGa                     >  1:1781528/1‑150 (MQ=255)
                                                                                                                              cacCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGAt                    >  2:3791177/1‑150 (MQ=255)
                                                                                                                              cacCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGAt                    >  2:1603155/1‑150 (MQ=255)
                                                                                                                              cacCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGAt                    <  2:3321454/150‑1 (MQ=255)
                                                                                                                               acCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGAtt                   <  1:2040324/150‑1 (MQ=255)
                                                                                                                               acCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGAtt                   <  1:2040131/150‑1 (MQ=255)
                                                                                                                                  gTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTTTCACTGTATTCGGCGATAAAAAAGTGGATTAGc                >  1:1914356/1‑150 (MQ=255)
                                                                                                                                  gTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGc                >  2:2797800/1‑150 (MQ=255)
                                                                                                                                  gTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGc                <  2:400730/150‑1 (MQ=255)
                                                                                                                                   tCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCGTCATCTTCGACCCCGCCCCCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCt               >  1:664498/1‑150 (MQ=255)
                                                                                                                                    cTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTc              <  2:1361558/150‑1 (MQ=255)
                                                                                                                                           tCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATATTCGACACCGCCACCGGCGTTATCGCGCTGATTGTAACAGTATTCGGCGATAAAAAATTTGTTTAGCTCAGTcaca       >  1:1169788/1‑150 (MQ=255)
                                                                                                                                            cACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCacac      >  2:2329429/1‑150 (MQ=255)
                                                                                                                                               gTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGc   >  2:3936715/1‑150 (MQ=255)
                                                                                                                                               gTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTAATATTCGACACCGCCACCGGCGTTATCGCGCTGATTGTCACTGTATTAGGCGATAAAAAAGTGGATTATCTCAGTCACACCGc   >  1:1492119/1‑150 (MQ=255)
                                                                                                                                                tataTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGcc  >  2:3450612/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCC  >  CP009273/2576178‑2576471

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